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LVI. The Thermal Conductivity of Potassium Chrome Alum at 
Temperatures below One Degree Absolute. 


By C. G, B. GARRETT, 
Royal Society Mond Laboratory, Cambridge *. 


[Received April 12, 1950] 


ABSTRACT. 


An account is given of measurements of the thermal conductivity of 
potassium chrome alum made between 0-14° and 0-30° K., from which 
it is concluded that the mean free path of the phonons in the crystal 
is about 4mm. In the method employed, a long single crystal of the 
salt, isolated from the outside world, is demagnetized to a uniform low 
temperature: a temperature gradient along the crystal is then set up, 
and the time constant of the subsequent heat flow measured. The paper 
concludes with a consideration of the lowest temperature to which a 
sample of material could in a reasonable time be cooled by direct thermal 
contact with a cold mass of paramagnetic salt. 


§ 1. [yTRODUCTION. 


TEMPERATURES below 1° K. may be obtained by the process of magnetic 
cooling, in which a paramagnetic salt is used as the refrigerant. In 
order to extend the study of these very low temperatures to other materials 
it is desirable to have information about the time required for cooling by 
direct contact with the cooled salt, and one of the controlling factors in 
this process will be the thermal conductivity of the paramagnetic crystals 
themselves. In the following paper an account is given of measurements 
of the thermal conductivity of a typical paramagnetic salt, potassium. 
chrome alum, in the range of temperature between 0-14° and 0-30° K. 
Such measurements are also of interest in connection with the theory of 
heat conduction in insulating materials. ; 

Precise measurements of the thermal conductivity of potassium chrome 
alum have been made in the liquid helium temperature range by Bijl 
(1949), but below 1° only the order of magnitude of the conductivity is 
known (Kiirti, Rollin and Simon 1936). The latter determination was 
made by demagnetizing a long single crystal of the salt from an inhomo- 
geneous field, so that the ends were cooled to somewhat different 


Oe ccc ee ee ee eae er aa = Tae anaes 
* Communicated by Sir Lawrence Bragg, F.R.S. 
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temperatures: the subsequent approach to temperature equality was 
watched by means of local magnetic susceptibility measurements made 
with mutual inductances placed around the two ends, and the conductivity 
estimated from semi-qualitative considerations. This method has the 
intrinsic advantage that it does not require separate thermometers or a 
heating coil attached to the crystal, thereby avoiding trouble from 
contact thermal resistances and heat inflow along metallic leads. In the 
present investigation a similar method has been used, in which however 
a device is employed for setting up a temperature gradient whenever 
required, without the necessity for a fresh demagnetization. In addition 
a more exact analysis of the heat flow has been made, and care has been 
taken to keep temperature differences as small as possible, with a view 
to investigating the temperature variation of the conductivity over the 
observable range. 


$2. PRINCIPLE OF THE METHOD. 


The salt, in the form of a long single crystal, is demagnetized from a 
uniform magnetic field, so that its final temperature is uniform. A 
non-uniform magnetic field, obtained from a suitable system of coils, 
is applied to the crystal in such a way that the field varies linearly along 
the length of the crystal with a gradient of the order of 100 gauss/cm. 
A temperature gradient is set up and heat flows from the warmer to the 
cooler end of the crystal. After allowing time sufficient for the 
temperature to become once more approximately uniform, the non- 
uniform field is removed, and a temperature gradient is now set up in 
the opposite sense to the first. The thermodynamics of this process will 
be discussed in §3. The temperatures near the two ends of the crystal may 
be measured by systems of mutual inductance coils suitably placed, 
and the rate of approach of those temperatures determined*. 

From the rate of approach of temperature it is possible to calculate 
the diffusivity of the material, K/pC, where K is the thermal conductivity, 
p the density and C the specific heat per unit mass. The value of C 
appropriate to the mean temperature of the crystal may be taken from 
previous measurements, and hence the thermal conductivity may 
be derived. Since, in the temperature range above 0:2°, C~A/T?, and 
since the conductivity itself is proportional approximately to T?, the 
diffusivity, and hence the equilibrium time of the flow process, are 
extremely rapidly varying functions of temperature. It is therefore 
possible, for a given size of crystal sample, to make measurements only 
over a limited temperature range. 


* In principle it would have been possible to do this during the course of the 
first heat flow, after applying the non-uniform field: in practice, however, 
fluctuations in such a field, to which the measuring system is very far from 
astatic, render accurate measurement impossible at this stage. 
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§ 3. THERMODYNAMICS OF THE TEMPERATURE GRADIENT INITIATION. 


The principle of the device for establishing the temperature gradient 
may best be understood in terms of the Mollier diagram for the salt 
(fig. 1), in which the entropy S is plotted against the enthalpy E, defined 
so that dE=T dS—MdH. For the purpose of constructing such a diagram 
it may be assumed that the susceptibility follows a Curie law y=A/T 
and that the specific heat is of the form C—A/T?, as is certainly a 
reasonable approximation throughout the temperature range considered. 
It may then be seen that ‘lines of constant field” are parabole, each 
one given by allotting to H a suitable value in the equation 


oe (ce eee weienese 2) 


“Lines of constant temperature”, on the other hand, are straight, 
being given by the equation 
Hie ii aes parte pie al at ohi ell) 
in which the appropriate values are given to T. These lines must 
necessarily stop at the parabola corresponding to He=0: 
2X2 
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Initially the crystal is assumed to be in a condition represented in 
fig. 1 by the point 0. The initial magnetization in the non-uniform field 
is carried out isentropically, the two ends of the crystal reaching point 
such as a and b, in different fields and at different temperatures. In 
the subsequent heat flow the total enthalpy remains unchanged, so that 
the ends will reach points c and d on the same isotherm, determined by 
the condition that a’c’=b'd’. (These considerations should strictly be 
applied to all pairs of corresponding points in the two halves of the 
crystal.) On removing the non-uniform field isentropically the ends 
move to points e and f, representing temperatures respectively higher 
and lower than that of the initial state o. Having achieved such a gradient 
in the presence of zero magnetic field, the subsequent rearrangement of 
temperature may be studied experimentally and the diffusivity deduced. 

After the second heat flow the final state of the system is represented 
by the point g, fixed by the condition that e’g’=f’g’. It may be verified 
that the temperature at g is somewhat higher than that at O, corresponding 
to a second-order increase in the entropy consequent upon the irreversible 
_ processes which have occurred. 

To calculate the order of magnitude of the effect and the exact 
distribution of temperature set up, we may proceed analytically, using 
the arguments outlined qualitatively above. If x represents the axial 
distance between any point in the crystal and one end, the total length 
being 2/, and if H(z) represent the field at that point, then it may be 
shown that the final temperature T’ at that point is given in terms of 
the initial uniform temperature T by the expression 


T= T{1— 55 [ Ho} —3|. A (a yp de |. ees) 


If the field H(~) varies uniformly along the length of the crystal following 
the equation H(~#)=H)+ ha/l, it being assumed that h<Hp, then 


T’T/1— AB ket). 5. 0 nn) 


The temperature gradient produced is therefore uniform. For potassium 
chrome alum A/A=1-36 . 10~® gauss“? (Garrett 1948). If H)=300 gauss 
and h=100 gauss, then the extreme difference of temperature produced 
at 0-2° should be, according to this formula, 0-008°. 


§ 4. EXPERIMENTAL DETAILS. 


The specimen was a single crystal of chrome alum, 65 mm. in length and 
of octagonal cross-section, the distance between opposite faces being about 
15mm. The crystal was cut in such a way that the axis of the sample 
coincided with one of the cubic axes of the crystal. It was suspended 
by fine nylon threads from a framework of german silver, which itself 
fits inside the conventional german silver salt-tube of a magnetic cooling 
cryostat. The measuring coils were wound on paper formers over the 
outside of the salt-tube, arranged so as to provide an estimate of the 
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averaged susceptibility over the last 15 mm. at each end of the crystal. 
Measurements of the two mutual inductances could be made alternately 
by comparison with a standard inductometer using a 40 cycle A.C. bridge 
network, detection of the balance point being made with the help of an 
amplifier and vibration galvanometer. In this way a temperature 
difference on the Curie scale of 0-01° at 0-2° could be determined with 
an accuracy of about 5 per cent. 

The initial cooling was attained with a water-cooled armoured solenoid, 
which gives a maximum field of 10 Kg. After demagnetization, the 
non-uniform magnetic field required for setting up the temperature 
gradient was produced by (i) a solenoid of internal diameter 120 mm. 
and field-factor 36 gauss/amp.; and (ii) a set of coils fitting inside the 
solenoid and around the liquid nitrogen dewar, designed to produce a 
linearly varying field of 10 gauss/amp./em. A current of 8 amp. passed 
through these two coils in series for an appropriate interval of time was 
found to set up a temperature gradient of a suitable magnitude. 

After switching off this non-uniform field, measurements of mutual 
inductance were made for the two coil systems alternately, the time of 
each determination being noted. Both sets of measuring coils were 
calibrated at temperatures in the helium range, so that mutual inductance 
readings could be used in each case to calculate T*. Unfortunately the 
corrections to T* for the shape of the crystal are unknown, and this 
uncertainty, although still moderately small at the lowest temperature 
investigated, will be a serious drawback in any attempt to extend the 
method to lower temperatures. As the field due to the primary of each 
of the mutual inductances is effective over roughly a spherical volume 
15 mm. in diameter near the end of the crystal, it has been supposed that 
the measurements will give Tx)... directly. 

It may be shown (see Appendix) that, to a good approximation, the 
difference in temperature between the ends of the crystal falls off _ 
exponentially with time according to the law 


AVorexp (a4), 2 OO ee Fae (5) 


in which, as before, 21 is the overall length and « the diffusivity of the 
crystal. In practice there will be, superimposed on these changes of 
temperature, a steady upward drift due to ‘ heat leak ’”’, which probably 
represents a surface recondensation of gaseous helium. Smoothed curves 
representing the terminal temperatures are drawn; from this graph 
another curve is plotted, showing the difference in temperature plotted 
logarithmically against the time. The linearity of this graph checks the 
validity of the exponential law (equation (5)) and the diffusivity of the 
crystal at a mean temperature may be estimated from the gradient of 
the fitted straight line. The relation between T* and the true temperature 
T is known from Bleaney’s measurements (Cooke 1949), as is also the 
specific heat at each temperature. The thermal conductivity itself may 
then be derived from the experimental values of the diffusivity. 
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§ 5. RESULTS. 


Measurements were made during three separate helium runs in 
October 1948. The results are given in the accompanying table. 


_——- ~ 


Ay? /Lar? K 10-°K 

Derren ee i meh a8 erg units Ts 
8 0-163 0-145 574 0-693 1-58 5-2 
6 0-167 0-149 478 0-682 1-86 5-6 
20 0-182 0-166 372 0-635 2-22 4-9 
20 0-232 0-218 112 0-475 5-52 5-3 
20 0-247 0-235 74 0-427 7-56 5:8 
20 0-252 0-240 74 0-405 715 5:2 
20 0-260 0-248 67 0-395 7-68 5-1 
20 0-276 0-266 42 0-358 11-2 5-9 
8 0-285 0-274 46 0-340 9-6 4-7 


| 


It will be seen that, as expected, the thermal equilibrium time does 
in fact depend very sharply on temperature. This, incidentally, explains 
why it has not so far been possible to extend the measurements to 
temperatures below 0-14° (where the state of transport of heat is no 
longer large in comparison with the heat leak) or above 0-30° (where an 
automatic recording device would have been necessary). The results 
are plotted against temperature in fig. 2, which shows that the conductivity 
is proportional to T? in this temperature range. It may be seen from the 
table that the average value of K/T® is 


537010" ergsisees, Sem: degree 7 te 2has eG) 


with a mean deviation of 0-12.105. Owing to the possibility of 
systematic errors the actual uncertainty is rather greater than this. 

Fig. 2 also shows the results of measurements made by Bijl (1949) 
in the liquid helium range, where the conductivity was found to follow 
a T?? law. Bijl found that the absolute magnitude of the conductivity 
appeared to depend on the rapidity with which the crystal was cooled 
from liquid nitrogen to helium temperatures. No attempt has been 
made to alter this rate of cooling during the present measurements, but 
it seems that our conditions correspond more closely with Bijl’s ‘‘ slow 
cooling ’’, and it will be seen from fig. 2 that a comparison of the two sets 
of measurements tends to confirm this. Kirti, Rollin and Simon (1936) 
estimated the thermal conductivity at 0-18° to be of the order of 10-4 watt 
units, and this also is in agreement with our results. 


§ 6. Discussion. 


In an insulating crystal the principal mechanism of thermal conduction 
is that associated with Debye waves or phonons. The phonon con- 
ductivity is given in order of magnitude by the expression 


Ky~1sC,, . . . . . . . . . (7) 
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Cy being the lattice specific heat per unit volume, s the velocity of sound 
in the crystal and J; the phonon mean free path. In this temperature 
range J, should be a length specifying the separation between lattice 
boundaries, and so should be independent of temperature; s is also 
independent of temperature, and, since the lattice specific heat C, should 
be proportional to T°, it would be expected that the phonon conductivity 
would also be proportional to 'T3. 


Fig. 2. 


10 


K 
(ergs seci! 
cm! deg) 


10° 


Temperature 


In § 5 it was shown that a T*® conductivity law is indeed followed for 
chrome alum at temperatures of the order of 0-2°K., although in the 
liquid helium range K ooT?:*. Taking C;~10?T® ergs/° (Casimir 1940) and 
s~105 cm./sec. we find that 

Loxton, ae een eae Be et aang (9) 
a length considerably smaller than the macroscopic dimensions of the 


erystal (6-5 cm. 1-5 cm.). Probably this is due to a mosaic structure in 
the crystal, as suggested by de Haas and Biermasz (1938) to account 
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for the “saturation” effect observed by them in experiments on the 
size dependence of the thermal conductivity of potassium chloride rods. 
Such a domain structure may also be associated with hysteresis effects 
in microwave experiments observed by Bleaney and Penrose (1948), and 
with the observation by Kraus and Nutting (1941) of a shattering of 
the crystal structure noticed on cooling a very similar material (ammonium 
chrome alum) through the temperature range 80—100° K. 

Akheiser and Pomeranchuk (1944) have made the interesting suggestion 
that at the lowest temperatures the conductivity of paramagnetic 
crystals should again become large, due to the alternative mechanism for 
the transfer of thermal energy provided by the spin system (“ exciton 
conductivity’). This is not in disagreement with the results presented 
in this paper, for calculation shows that at 0-14° K. the exciton conduc- 
tivity should be some 5000 times smaller than the phonon conductivity. 
At somewhat lower temperatures the thermal equilibrium times become 
exceedingly long, but, if the predictions of Akheiser and Pomeranchuk are _ 
indeed correct, the diffusivity at 0-003° K. should be equal to that at 
0-14° K., while at still lower temperatures it ought to increase as 1/T?. 
In view of the technical interest of such a possible increase in conductivity 
it is obviously of considerable importance that experiments be extended 
to this lower range of temperature. 

From the information now available, it is possible to arrive at a 
tentative conclusion concerning the lowest temperature to which it is | 
possible in a reasonable time to cool a sample of material by thermal 
contact with a cooled mass of chrome alum. It may be seen from 
equations (5) and (7) that the thermal equilibrium time within a crystal 
is proportional to T~*/7J,-1, where T is the temperature, / the macroscopic 
size of the sample and /,the phonon mean path. In an actual experiment 
(see table of results) an equilibrium time of 500 seconds was found when 
I~65 mm. and T=0-15°, 1, being subsequently deduced as 0-5 mm. 
From these figures it is possible to calculate what would be the equilibrium 
time at that temperature under the conditions of an actual contact 
cooling experiment ; and hence to deduce how low the temperature 
could be before the equilibrium time would become intolerably long. 
In the experiments of Mendoza (1948) in which powdered salt was 
compressed around a copper insert under considerable pressure, an 
acetone-based lacquer being used as bonding agent, it was possible to 
reduce the distance between any grain of salt and the nearest copper 
surface to about 2mm.; but on the other hand the phonon mean free 
path is probably given now by the size of the crystal grains, a distance 
very much smaller than the phonon mean free path of 0-5 mm. found for 
alarge crystal. Asa basis for comparison put /=1 mm. and 1,=5:10-3 mm. 
Using the above analysis it is found that at 0-05° K. the equilibrium time 
for such a sample should be of the order of 1000 seconds, which is about 
the longest time tolerable in a contact-cooling process. It may therefore 
be concluded that a temperature of 0:05° K. probably represents the limit 
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to which it would be practicable to cool a sample by direct thermal 
contact with a cooled salt. In this calculation no account has been taken 
of any surface thermal resistance between the salt and the metal. This 
assumption is based on the conclusion of Mendoza (1948) that below 
0-2° K. the surface resistance is not the controlling factor. 
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APPENDIX. 


The exact solution of the thermal conduction problem, taking into 
account a diffusivity varying as T5, would involve elaborate numerical 
computation. It is sufficient however to suppose the diffusivity constant 
at some mean value, and the problem is then soluble by standard methods 
(Carslaw and Jaeger 1947). Since radiation from the sides is of course 
negligible, the crystal may be treated as an infinite slab between two 
parallel planes, there being no heat flow across either boundary. If the 
initial temperature distribution be f(z), then the temperature at some 
later time ¢ at any point x is given by 

2 » 21 ’ 

T= xi), f(x’) da’ + 72 "cos = aT) COs 7 aids aA) 
where « represents the diffusivity of the material. Suppose a linear 
initial temperature distribution of the form f(«)=T)+s(v—l).  Sub- 
stituting in equation (9), all terms under the summation sign for which 
n is even vanish, and there remains the expression 


8sl THX reat 
T=T,— 2) COS (3) exp (- <) 


8s] 3arx i Ont 
— rae ie exp 4/2 


ae PR cs fens eee | L()) 


It is possible, however, to neglect all except the first of the exponential 
terms, provided that we are not interested in the distribution at times 
very shortly after the initiation of the heat flow. Since the coefficient 
of ¢ in the exponent in equation (10) is the same at all points along the 
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length of the crystal, it is immaterial whether temperature measurements 
are made over an infinitesimal region at the very end of the crystal, or 
whether instead an average is observed of the temperature near the end 
of the crystal. In either case 


AL ~exp (=amtj4i), <= ns) eee 


so that a graph of log T against t should be a straight line, from the 
gradient of which x may be found. 
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ERRATA. 


In the author’s paper on Quantal Aspects of Scientific Information, which 
appeared in the March issue (p. 294), the relation between physical scales and 
proper scales is unfortunately obscured by an error in the example of §4(c). 


The last two sentences on p. 294 should read : “ We can now give a probability 
p of just } to a statement in the form: a occupies Vi equal intervals. The 


figure of } implies . . . . that its terminal is as likely to fall within the Vith equal 
as to fall outside it.” 


As pointed out in §4(c), the number of conceptually-separate occupance- 
relations, or metron-content, (§4(b)) is not the same thing as the number of 
physically-separate occupance-relations between a magnitude and the scale on 
which it is observed (§4(a)). Logically, 7 metrons enable us to add i elements to 
a formal representation of a measurement. The number of occupance-relations 
on the scale of physical magnitude graduated as in §4(a) may be quite different. 
In the example given, (§4(c)) i=2?/(4dz?). But the physical scale-unit of x as 
defined in §4(a) is 4x, and if the probable error is constant, only x/Ax equal 
intervals can be observationally separated on the scale of x, though i concep- 
tually significant intervals can be defined. The probability is 1/2, not that the 
ith conceptual interval be occupied, as stated, but that the Vith physical scale- 
unit be occupied. The number of occupance-relations on the physical scale 


is, infact, of greater practical importance, though it lacks of course some of 
the invariant features of metron-content. 


D. M. Mackay. 
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LVII. The Crystallography and Thermodynamics of Order—Disorder 
Transitions in Certain Mixed Ketones. 


By VeRA DANIEL, 
British Electrical Research Association*. 


[Received March 20, 1950.] 
[Plates XX. & XXI.] 


SUMMARY. 


Previous work by Oldham and Ubbelohde (1939, 1940) on ketones of 
equal chainlengths and different positions of the ketogroup has shown 
that some of these compounds have practically identical crystal structures 
and physical properties; these will be called equivalent ketones. In 
mixtures of equivalent ketones all constituent ketones have identical 
carbon chains, and each ketone contains one ketogroup. Thus, mixtures 
of equivalent ketones differ from each other only in so far as the con- 
figuration of the ketogroups affects the bonding of their structure. In 
the solid state, the system formed by mixtures of different proportions of 
equivalent ketones is analogous to a crystalline compound, the bonding 
of which may be altered at will, within certain limits. 

In the present paper the thermodynamic stability of different con- 
figurations of ketogroups in mixtures of equal amounts of equivalent 
ketones is treated theoretically. It is assumed the ketogroups contribute 
to the bonding only by virtue of their dipole moment. 

Different mixtures of six equivalent ketones of chainlength eighteen 
were examined experimentally. They always form one of two structures, 
one orthorhombic (or monoclinic) and the other hexagonal, the hexagonal 
one being relatively disordered, as far as the packing of the chains is 
concerned. Two kinds of order—disorder transitions occur, one of which 
corresponds to a segregation of the individual ketones into separate 
phases, while the other one is analogous to the well known transitions in 
paraffins. The conditions of stability of the different configurations of 
ketogroups agree with the theory. In particular, the binding energy 
of a plane array of ketogroups can be derived from transition temperatures, 
and agrees with the electrostatic interaction energy derived for this 

array by Froéhlich (1946). 


INTRODUCTION. 


THE relationship between the physical properties of substances and 
their chemical and crystallographic structure is difficult to investigate 
because their structure cannot be modified at will. A physical property 


* Communicated by Sir Lawrence Bragg. 
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such as, say, dielectric constant or thermal conductivity, depends on a 
number of factors. The understanding and control of such a property 
would be easier if these factors could be altered separately and in small 
steps. However, this is usually not possible because a change from one 
chemical compound to another generally alters a number of structural 
features simultaneously, and in rather large steps. 

Systems ‘which permit the alteration of structural features in small 
steps can arise in two ways. Firstly, there exist series of compounds 
such that structural features within the series vary in a regular way and 
in more or less small stages. Such are the series of homologues in organic 
chemistry, for instance the paraffins, or the groups of analogous inorganic 
compounds such as the alkali halides. Here the modification of the 
structure is effected by chemical means and may be restricted to a more 
or less distinct structural feature, such as the length of a paraffin chain. 
Secondly, systems with variable structure may be produced by mixing 
two or more chemical compounds. Here changes can be made indefinitely 
small, but they are generally not restricted to any one particular structural 
feature. 

In the present investigation both the above methods are applied 
simultaneously, in that the materials investigated are mixtures of closely 
related compounds, where both the original compounds and their mixtures 
have very similar crystal structures. These mixtures form a system 
which can in some ways be varied in small and quantitatively controlled 
steps. It may exist in forms of different degrees of order, and the 
conditions of stability of these forms can be derived thermodynamically, 
on the basis of certain simple assumptions. 


THE STRUCTURE AND PROPERTIES OF LONG-CHAIN 
KETONES. 


The compounds used for the present work are long-chain ketones 
of the general formula C,H,,,,,CO C,H,,,,, with p and q greater than 3. 
They differ from normal paraffins only in that two hydrogens are replaced 
by an oxygen. The C=0 group hasa large dipole moment, and a ketone 
of this type may be considered as a paraffin to which a dipole has been 
rigidly attached, approximately at right angles to the chain. 

The structure of normal long-chain paraffins has been fairly fully 
investigated in the past. Miiller (1928, 1930) determined the crystal 
structure of the paraffin C,)H¢ (see fig. 1). This structure is typical 
for the long-chain paraffins and many of their derivatives. Its main 
features are that each molecule forms a comparatively rigid flat zig-zag 
chain and that the chains are parallel to each other in layers. A projection 
of a layer of chains on the basal plane shows a face-centred arrangement 
(fig. 1(a)). There exist different modifications of this structure. The 


chains may be normal or inclined to the basal plane, and the fit of the 
chain-ends may vary. ; 
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Long-chain paraffins show a variety of polymorphic transitions (Piper, 
Chibnall and Collaborators 1931). Many of these transitions are connected 
with changes of the angle between the chain axes and the basal plane of 
the structure, while in other cases the chains are normal to the basal plane 
throughout the transition. The latter type of transition was investigated 
by Miiller (1932) with the help of X-rays. Miiller found that the transition 
consists in a very anisotropic expansion of the basal plane ab (fig. 1(@)) such 
that the chains in it approach a hexagonal arrangement. That is the 
cross section of the chains becomes statistically circular rather than flat. 
This change of structure is gradual in some paraffins and discontinuous 
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The structure of the paraffin Cy9H¢o. (Simplified from Miiller 1928.) 


in others. It represents a typical order—disorder transition, as is also 
confirmed by the calorimetric data of Garner, Bibber and King (1931) and 
Ubbelohde (1938). Dielectric measurements by Jackson (1935a, b) Miiller 
(1936, 1939) Pelmore (1939) and Sillars (1939) on paraffins with dipolar 
groups indicate that in some paraffin lattices the chains have considerable 
freedom to rotate and/or twist. 

As regards ketones as distinct from paraffins, a number of ketones 
CpHopi1 CO CgHogy1 can be produced for any one long-chain paraffin 
of a given chainlength p+q-- I. Oldham and Ubbelohde (1939) measured 
the melting points of all ketones of chainlengths seventeen and eighteen, and 
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in addition measured the latent heat of melting and entropy of melting of 
all ketones of chainlengths seventeen. Table I. gives their results in a 
slightly simplified form. It shows that for all ketones of the same 
chainlengths, with p and q greater than three, all the measured quantities 
are very similar. The only exception is C,H,,COC,H,, which differs from 
the other ketones in that its ketogroup is symmetrically placed in the 
centre of the chain. All the unsymmetrical ketones of equal chainlengths 
and with ketogroups away from the chain-ends have closely similar latent 
heats and entropies of melting and melting points. 


TABLE I. 


Ketones of total chainlength p+q+1=17. 


Tx AH AS 
P OF k cal./mole cal./mole/° 
it 47-8 14-1 44-0 
2 47-6 13-9 43-2 
3 41-6 12-7 40-4 
+ 41-5 eT 40-3 
5 40-9 12-7 40-5 
6 40-6 12-9 40-9 
7 41-8 13-1 41-7 
8 50-0 is: 42-3 


Ketones of total chainlength p+q+1=18 


agi 

1 51-5 

2 51-0 

3 45-8 Ty ....melting point, 
+ 45-2 4H ....latent heat of 
5 44:8 melting, 

6 45-0 AS ....entropy of 

i! 44-6 melting. 

8 45-8 


Oldham and Ubbelohde (1939) examined the ketones of chainlengths 
seventeen by X-rays and microscopically without finding significant 
differences between them. A similar examination was carried out by 
the author for ketones of chainlengths eighteen*. 


ee Se PE Ce SE ce ee mee aaa ereE ey Mite se yee) ST 

* The author is greatly indebted to Dr. Oldham of the Royal Institution for 
the loan of small samples of five ketones of chainlengths eighteen with p—3, 4 
5, 7 and 8. A sixth ketone, with p=6, was synthesized according to 
Dr. Oldham’s advice. 
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X-ray powder photographs were taken of six ketones of that chainlength 
and C,H,,COC,,H., was investigated in somewhat greater detail, using 
single crystals; but no complete structure determination was made. 
All the ketones have unit cells of orthorhombic shape, and crystallize 
with their chain axes normal to the basal plane. (This does not, of 
course, exclude a monoclinic space-group, and that of C,H,,;,COC,,H,, 
is probably C§, or Ct.) No significant difference was found between 
the unit cell dimensions of the six ketones and inspection of the powder 
photographs indicates that in all the ketones the chains are arranged 
very similarly. The X-ray evidence will later be discussed in more 
detail. 

The structure of an unsymmetrical ketone of the kind discussed here 
has some interesting features which have been considered by Fréhlich 
(1946)*. As fig. 2 shows, the ketogroups belonging to a layer of chains 
may be situated at one of two levels. Within either of these two levels 


Fig. 2. 


Schematic representation of a layer in a crystalline ketone. The positions 
which may be occupied by oxygen atoms are marked by crosses. 


0 ee et Ve Ds el a eS eee 

* Fydhlich’s paper is mainly concerned with a theory of order—disorder 
transitions in pure ketones. Miiller’s (1936, 1938) dielectric measurements for 
the ketones C,H,,COC,H,, and C,,H,;COC,,H,; appeared to indicate the 
inception of such a transition near the melting point : they show the dielectric 
constant to change gradually, over a temperature range, from the value for the 
solid to that for the liquid. This gives the impression that a transition would 
take place at a temperature above the melting point if the material did not 
melt before this temperature is reached. However, dielectric measurements on 
very carefully purified samples of CsH,,CO CsH,,, CeHis COC,,H,; and 
C;H,,COC,,H,; indicate that in pure solid ketones the dielectric constant is 
almost independant of temperature, and that there is no gradual transition 
from the value for the solid to that for the liquid. The change from solid to 
liquid in a pure ketone is discontinuous at the melting point. and the gradualness 
of the transition found by Miiller appears to be due to small traces of impurities. 
Thus there is no evidence that pure ketones show order—disorder transitions. 
The dielectric measurements mentioned will be described elsewhere. 
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the ketogroups form a plandr arrangement, such that the dipoles 
represented by the ketogroup are roughly parallel to each other. 
Frohlich shows that this arrangement is very favourable from the point 
of view of electrostatic interaction. If a level is fully populated, so that 
each chain contributes a dipole to it, the electrostatic interaction of the 
dipoles contributes quite a large term to the binding energy of the 
structure. 

The arrangement of the paraffin chains in the structure leaves it 


undecided whether ketogroups shall be concentrated in “dipolar ~ 


planes” or distributed at random between the two possible levels. 
However, the two possibilities imply different energy conditions. The 
random arrangement has a higher entropy. On the other hand, the 
electrostatic interaction is stronger if the ketogroups are concentrated 
at one level, because the electrostatic interaction energy between dipoles 
varies inversely with the third power of their distance. If the dipoles 
are distributed between two levels the numerical value of the total dipolar 
interaction energy falls to about a quarter of that for a fully populated 
dipolar plane. If the dipoles are concentrated in dipolar planes the 
structure will have a lower entropy and a lower energy than if they are 
distributed at random between two levels, for the dipolar interaction 
reduces the total energy H of the structure by restricting some of the 
degrees of freedom. Thus the free energy G=H—TS of the form with 
dipolar planes has a lower H and a lower 8 than the random form. 
Whether or not it has a lower free energy depends on the numerical 
value of the energy of the dipolar plane. This energy is, as it were, 
balanced against the entropy term. If it is sufficiently large, the form 
with dipolar planes will be stable. The thermodynamic conditions of 
stability in question will be discussed in more detail later. 

It is a matter for experiment to decide which of the possible forms 
of an unsymmetrical ketone is in fact formed. From the point of view 
of X-ray crystallography the two forms, with and without dipolar planes, 
differ very little from each other. A preliminary investigation of single 
crystals of the ketone C,H,;COC,,H,; showed that it would need a 
detailed determination of the crystal structure to decide which of the 
two possible forms is present. However, the two arrangements differ 
radically in their ability to form solid solutions with other, similar, 
ketones, as will be seen from the following. 


MIxTURES OF EQUIVALENT KETONES, 


As Oldham’s and Ubbelohde’s work (1939) shows, unsymmetrical 
ketones of equal chainlengths but different position of the ketogroup 
have practically identical crystal structure and thermal properties unless 
their ketogroups are near the chain-ends. These ketones, of the general 
formula C,H,,,,COC,H,,,, with p>3q>3 and p+q, will be called 
equivalent ketones. Because of the similarity of the “ equivalent ” 
ketones one would expect them to form mixed crystals which were very 
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nearly ideal solid solutions. However, this cannot be the case if the 
ketones form dipolar planes with a high energy of electrostatic interaction, 
because in a solid solution these planes must be brokenup. Consequently, 
if the energy of a dipolar plane is large, two equivalent ketones may be 
mutually immiscible, in spite of their similar structure. On the other 
hand, if no dipolar planes are formed in the pure peers ketones, 
they should form solid solutions with each other. 

Oldham and Ubbelohde (1940) in a study of structures with flaws, 
measured the melting points of a number of binary mixtures of equivalent 
ketones. Their results show that in thermodynamic équilibrium these 
ketones do not mix in the solid to any appreciable extent. This indicates 
that the pure equivalent ketones crystallize with their ketogroups segre- 
gated in dipolar planes. 

The two-phase character of binary mixtures of equivalent ketones 
shows that the energy of dipolar planes must be rather large for the 
following reason. In pure ketones there exist two levels among which 
the ketogroups might be distributed at random, while the configuration 
of the paraffin chains is retained. (See fig. 2.) In a mixture of ‘two 
equivalent ketones with different positions of the ketogroup there are 
four such levels, and the formation of dipolar planes means a greater 
loss of randomness than in a pure ketone. This means that in a mixture 
of two equivalent ketones the energy of the dipolar planes has to outweigh 
a larger entropy term if it is to enforce the segregation of dipoles in 
dipolar planes. The fact that the segregation does occur in binary 
mixtures implies that the energy of the dipolar planes is sufficient even 
for that. However, what will happen on mixing three, four, or more 
equivalent ketones? In these cases there will be still more levels 
available among which the ketogroups can be distributed at random, 
and there must be a number of ketones for which the energy of the dipolar 
plane will be insufficient to enforce the segregation of the dipoles and a 
solid solution will be formed. 


Tur THERMODYNAMICS OF MIXTURES OF EQUIVALENT 
KETONES. 


The qualitative reasoning in the preceding paragraph may be made 
quantitative, subject to some simple assumptions. It may however be 
useful first to restate the problem in somewhat different terms. 

There exist many substances which crystallize at low temperatures 
in an ordered structure which changes to a disordered structure as the 
temperature rises. In all these cases an energy term—making for 
order—militates against an entropy term which makes for disorder, 
The entropy term contains the temperature as a factor while the energy 
term does not and so the entropy term must outweigh the energy term 
ait some sufficiently: high temperature. - , 

The case of two substances which in the solid state tend to segregate 
is analogous to: an order-disorder transition, Here again an energy 
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term makes for segregation while an entropy term makes for solid solution. 
Two such substances are said to form a solid solubility gap which closes 
at some temperature, so that above that temperature a solid solution is . 
always formed. 

Mixtures of equivalent ketones have the characteristics both of a single 
substance undergoing an order-—disorder transition and of a number of 
substances forming a multiple solid-solubility gap. The fact that we 
have a number of equivalent substances makes it possible to increase 
the entropy term—which makes for disorder and solid solution—in 
discrete steps. The disordered form may be made stable by increasing 
the number of ketones instead of the temperature. 

In the following we consider mixtures of equal amounts of n equivalent 
ketones and deduce their conditions of thermodynamic equilibrium. 
To this end we make the following assumptions : 

1. The ketones considered have equal total energies and entropies 
for equal temperatures and volumes when in equilibrium. 

2. Liquid mixtures of the ketones are ideal solutions. 

The first assumption depends for its justification on the equality of 
the latent heats and entropies of melting of ketones of chainlength 
seventeen (Oldham and Ubbelohde 1939) and the similarity of the 
structure of the ketones. The second assumption implies that in liquid 
ketones the ketogroups are distributed at random. 

For solid mixtures of » ketones, Gibbs’ phase rule implies that the 
mixture can consist of not more than ” phases in equilibrium (except at 
a so-called n—2 fold point). With regard to the nature of the possible 
phases we make the assumption : 

3. Any ketone in the solid mixture may crystallize either with or 
without dipolar planes. That is, the ketogroups may either all be 
concentrated in one level or else distributed at random between all the 
2n possible levels. This assumption implies that the solid mixture may 
be either n-phase or single-phase, since if the form with dipolar planes 
is stable for one ketone, it must be stable for all the ketones. This, 
because equal amounts of the m ketones are present and thus the specific 
volume is equal for all of them and because they are equivalent as defined 
in the first assumption. If the form with dipolar planes is stable the 
mixture will be n-phase, each phase consisting predominantly of one 
ketone. If the form without dipolar planes is stable the ketones will 
form a single-phase solid solution. We shall make one more simplifying 
assumption with regard to these two possibilities : 

4. Ifthe mixture is single phase, it is a completely random solid solution. 
If it is n-phase, the mass of each ketone is concentrated in its “‘ own ” 
phase, the remainder of its mass which is distributed among the other 
n—1 phases being negligible. 

In principle, conditions are likely to exist in which each ketone contains 
a fair quantity of the other ketones in solid solution, while still forming 
a.separate phase. However, such ‘conditions are. considered to be 
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restricted to narrow ranges of temperature so that they may be neglected. 
This is equivalent to assuming the transition from the form with dipolar 
planes to the form without dipolar planes to be a transition of the first 
order. 

In a mixture of equal amounts of n equivalent ketones all the ketones 
behave alike and so it is sufficient to write down the thermodynamic 
functions for one of them. 

For a pure ketone we put the Gibbs free energy per mole in the solid 
state 


Cp il 6 ate We Pca re ee ee AL) 
and in the liquid state 
(GS el (Sole ce se sk een ar eee) 


the notation used is that of Oldham and Ubbelohde (1940). G is defined 
for constant pressure so that H contains a volume term. 

As the liquid mixture of 2 equivalent ketones is assumed to be an 
ideal solution, the partial total energy of each equivalent ketone is the 
same as for a pure liquid ketone, but the partial entropy is increased 
by a term equal to R log n because the volume available for each ketone 
is increased n fold (see for instance Fowler and Guggenheim, Statistical 
Thermodynamics, p. 161). Thus for the liquid mixture of » ketones 


(Gie=(fE)— 11 (S.)4- Relop win ee te fee) 


The solid mixture may be either n-phase or single phase. If it is 
n-phase the assumptions made imply that all except a negligible part 
of each ketone is in the same state as in a pure ketone. 


Thus Ge =H. ees hE eh PH C4) 


If the solid mixture is single-phase, the partial entropy of each ketone 
will be increased by a mixture term, which is again R log n if the ketones 
form ideal mixed crystals (see Fowler and Guggenheim, Statistical 
Thermodynamics, p. 185). In addition there will be an entropy increase 
of R log 2, due to the fact that in a pure unsymmetrical ketone all dipoles 
are on one level, while in the mixture the dipoles of any one unsymmetrical 
ketone are on two levels. Apart from this configurational increase of 
entropy R log 2n there may be present a further increase of entropy 
because the removal of the dipolar planes is likely to make the whole 
structure looser. This dynamic part of the entropy is difficult to 
estimate. We shall treat it as an unkown constant s, independent of n. 
This is most likely an over-simplification, as the dynamic part of the 
excess entropy probably depends on the number of levels on which dipoles 
may occur, and thus on n, but it seems reasonable to neglect this 
dependence, in first approximation. The total energy H of the mixture 
will be increased by the energy of the dipolar plane hy, and if there is a 
further loosening of structure beyond the removal of the dipolar planes, 
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it may be increased by a further term h,. Both these terms are 
independent of », in first approximation. This makes the free energy 
of the single-phase mixture 


Gg Hip thy TiS Blog 202)s a) a) 


The assumptions implicit in this equation mean that the single phase 
mixture may have another structure than the pure ketone—hence h, 
and s—but that the thermal characteristics of this structure are 
independent of n. 

The equations which have been deduced permit calculation of the 
melting points of the n-phase and single phase mixtures. For the 
melting point T, of the n-phase mixture G,=(G), so that 


(Hy) —H)=T,((89)— 8) + RT, log n 


where (H,)—H,)= AH is the latent heat of melting of a pure ketone and 
(S,)—S)=AS its entropy of melting. Using the equation 4H=T,4S 
valid for the melting point T, of the pure ketone, the melting point of 
the n-phase mixture is given by 


1 i 
AH (= — qe) =R log, rere Fit 


Equation (6) is simply the Hildebrand equation for the melting points 
of inmiscible substances, written forthe special case of equal amounts 
of substances of equal 4H and AS. 

The melting points of the single phase mixtures follow from (G)=G, 
and should obey the equation 


1 1 hp+h 
AH (7 — 7) — a ==R log 2-8"... See 
This implies that T,, the melting point of the single phase mixtures, 
should be independent of n. 

Comparison of equations (6) and (7) shows that with increasing » the 
melting point of the n-phase mixtures decreases while that of the single 
phase mixtures remains constant. This implies that even if the multi- 
phase mixture is thermodynamically stable up to the melting point 
for low values of 7, it must: become unstable as compared with the single- 
phase form if » is made large enough. The condition in which the two 
forms are in equilibrium. at a temperature T can be derived by. putting 
G,=G, which leads to the equation 


hy t+h,=T(R log 2n+s), . . . . . . . (8) 


The lowest value for » at which the single phase form is stable in the 


solid is obtained if T is the melting point of the single phase form, that 
Is wie bes 
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MeLTING Points or Mixep Kerongs. 

Equations (6) to (8) derived in the last paragraph may be compared with 
experiment by measuring the melting point of mixtures of equivalent 
ketones. By the courtesy of Dr. Oldham, the author had available 
small samples of all six equivalent ketones with chainlengths eighteen. 
Mixtures were prepared of equal amounts of these, and their melting 
points were measured in the usual way, the mixtures being held in glass 
capillaries. These measurements led to the observation that mixtures 
of 2, 3 and 4 ketones can exist in a metastable state. If the mixtures 
are allowed to solidify slowly, and kept at room temperature for some 
time, they reach a stable condition to which corresponds a certain 
“melting point. If, however, they are quenched in cold water and then 
quickly immersed in warm water, they melt at a lower temperature, 
a metastable melting point. If the liquid mixture is kept for some time 
at a temperature between the metastable and the stable melting points, 
it ultimately solidifies again in the stable form. These observations show 
that, on quenching, the mixtures in question solidify in a mestastable 
form which changes in the solid state into the stable form if given time. 
The persistence of the metastable form varies markedly with the number 
n of ketones in the mixture. Mixtures of four equivalent ketones 
persist in the metastable state for several days, while in mixtures of 
two ketones the metastable melting point could only be observed if 
the quenching in water and ice was followed very quickly by immersion 
in warm water. No metastable melting points were observed in single 
ketones. 

There exist in general several combinations of equal amounts of ” 
ketones for any number of n. With six equivalent ketones of chainlength 
eighteen, the number of combinations is, respectively, 15, 20, 15, 6 and 1 
for mixtures of 2, 3, 4, 5 and 6 ketones. It would have needed pro- 
hibitive amount of material to make all possible mixtures. Three 
combinations were made of four ketones, and it was found that the values 
measured for the stable and the metastable melting points did not 
depend noticeably on the individual ketones making up the four fold 
mixtures. Variations found. were within a range of 1°C. For other 
mixtures only one or two combinations were used for each value of n. 

The results of the measurements are shown in fig. 3. To facilitate 
comparison with the theoretical equations 6-8, the reciprocals of the 
melting points are plotted against the logarithm of n. The theory 
postulates that 1/T should vary linearly with log n if the mixtures are 
n-phase, while for single-phase mixtures 1/T should be independent 
of n. It may be seen that for n=1to4 the stable melting points lie 
on a straight line (1). The metastable melting points depend much less 
on n and lie on another curve (II), together with the stable melting points 
of the mixtures of five and six ketones. The experimental results can 
be interpreted in terms of the theory, if it is assumed that curve I 
corresponds to n-phase mixtures, while curve II corresponds to single- 
phase mixtures. This means that up to n=4-6 where the two curves 
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intersect, the n-phase state is thermodynamically stable while the single 
phase state is stable for higher values of 1 (at the melting point). The 
existence of metastable single-phase mixtures can be easily understood 
because the formation of a multi-phase mixture from the melt involves 
extensive sorting out of the mixed ketones. This segregation is a process 
of diffusion and needs time. Thus, if the mixture is quenched it has to 
solidify in the random single-phase form. This interpretation explains 
why mixtures of four ketones persist in the metastable state for several 
days, while mixtures of two ketones become stable in a few seconds 
(the temperature being equal in both cases). The segregation of four 
ketones should need more time than the segregation of two. 


+ stable melting points 
O metastable melting points 


0-7 0-8 1-0 
Melting points of mixtures of n equivalent ketones. 


According to equation (6) the slope of curve I should give 4H, the 
latent heat of melting of a pure ketone. The value found from the 
experimental data is MH=16-6k cal/mole. Oldham and Ubbelohde 
(1940) give the latent heat of melting of one of the equivalent ketones 
used in the present work, namely C;H,,COC,,H,;, as AH=14/6k cal/mole. 
The agreement is satisfactory, as equation (6) embodies some assumptions 
which are unlikely to be strictly true. The value of Rlogn for the 
entropy of mixing of the liquid mixture of ketones is an upper limit, 
and this term will be less if assumption 2 is not obeyed and the liquid 
mixture of ketones is not an ideal solution. If that is the case the value 
of JH deduced from equation (6) is spuriously large. Similarly, it is 
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presumably an exaggeration to assume that the ketones are entirely 
immiscible in the solid multiphase form. If they are somewhat soluble 
in each other, equation (6) will again givea spuriously high value of 4H. 
As both the factéts mentioned tend to make AH too large, the agreement 
found shows that the assumptions made are reasonably correct. 

Equation (7) postulates that the melting points of the single-phase 
mixtures should be independent of ». Experiment shows that these 
melting points do not depend very strongly on n. However, they are 
not constant ;. mixtures of fewer ketones have higher melting points 
than mixtures of larger numbers of ketones. This is not surprising, 
for the constancy of h,, hy and s in the equations could hardly be 
expected to be more than a crude first approximation. Some dipolar 
interaction must occur even in random mixtures and this will be more 
marked the more densely populated the levels containing dipoles. Dipolar 
interaction should iné¢rease the melting points of single-phase mixtures 
of few ketones, where the levels are relatively densely populated, in 
comparison with mixtures of many ketones. 


Tue CrystaL StrucTuRES oF Mixep KETONES. 


The evidence of the melting points indicates the existence of two 
forms of the mixed ketones, a multi-phase form and a single-phase 
form. It remains now to investigate the structures of these forms, with 
a view to interpreting the energies and entropies involved in the 
order—disorder transition. 

X-ray powder photographs were taken of a number of mixtures of 
the six equivalent ketones of chainlengths eighteen, both in the stable 
and the metastable state. As the most important Debye—Scherrer 
lines of the ketones occur at Bragg angles of 2 to 20°, an experimental 
arrangement had to be used which gave sharp lines for these angles. 
A Unicam single-crystal gonimeter was used, in conjunction with a slit 
collinator, the slits being about 0-2mm. wide. The specimens were 
enclosed in cellophane tubes. For the measurement of spacings flat 
film was used, the specimen-film distance being 10 cm., and for a survey 
of the diffraction pattern a cylindrical filmholder of 3 cm. radius. 

The results of the X-ray examinations were as follows: In all cases 
when according to the melting point data an n-phase mixture was 
expected, its X-ray diffraction pattern was consistent with a multi-phase 
mixture of the component ketones (see Pl. XX.(6)). In all cases when a 
single-phase mixture was expected, it had an hexagonal structure, 
different from that of the pure ketones (see Pl. XX. (qa)). The structures 
of all single-phase mixtures were identical, within the limits of error. 
In the following these results will be discussed in greater detail. 

The six equivalent ketones have closely similar X-ray patterns in 
that they have identical unit cells and the intensities of most Debye- 
Scherrer lines are very similar. The reason for this is that all the 
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ketones have identical paraffin chains which are apparently arranged 
very similarly in the different ketones. As the ketones contain only 
one oxygen atom for a paraffin chain of eighteen carbon atoms, the 
intensities of the Debye-Scherrer lines are mostly determined by the © 
paraffin chains. Only if the structure amplitude of the paraffin chain 
for an X-ray reflection is small does the position of the ketogroup become 
important for the intensity of that reflection. This happens to be the 
case for the 00/ reflections with J up to about 17, and as the 001 reflections 
up to J=5 are clearly distinguishable on powder photographs, it is easily 
possible to identify a ketone from a powder photograph. Otherwise, 
however, the X-ray diffraction patterns of the ketones are essentially 
due to their paraffin chains. 

The multi-phase mixtures have orthorhombic unit cells, the dimensions 
of which are closely similar to those of the pure ketones (see Table II. a), 
The signficance of the small discrepancies found is not easy to assess, 
because, as has been pointed out by McArthur (1944), even with exceedingly 
pure ketones and paraffins unit cell dimensions are not reproducible 
within the geometric limits of accuracy. The powder patterns are 
too complex to be analysed in detail, but as far as visual examination 
can show, the diffraction patterns of the multi-phase mixtures are very 
similar to those of the pure ketones with regard to the intensities of the 
reflections other than the first five 001 reflections. This means that the 
arrangement of the paraffin chains in the multi-phase mixtures is very 
similar to that in the pure ketones. This agrees with the prediction 
that these mixtures should be multi-phase mixtures of the pure ketones. 

As for the postulate of the theory that the individual ketones in the 
multi-phase mixtures should form separate aggregates with dipolar planes, 
the X-ray data are insufficient to give a conclusive answer. The only 
reflections which are different for the individual ketones, and could be 
identified on powder photographs, are the first five 00/ reflections. The 
intensities of these reflections are only little dependent on whether the 
individual ketones form separate aggregates. However, even if more 
data were available, the investigation of this problem would meet the 
following difficulty. 

In general, multi-phase aggregates of crystals may be expected to 
scatter X-rays incoherently, because the regularity of the repeat of the 
crystal structure will normally be interrupted at the boundaries between 
aggregates. However, this need not be the case with a multi-phase 
mixture of equivalent ketones, where the individual ketones have 
identically arranged paraffin chains. Aggregates of the different ketones, 
each containing dipolar planes, may be packed together so as to scatter 
coherently, as if the whole conglomerate were one single crystal. If 
this were the case, only a very detailed X-ray investigation could decide 
whether dipolar planes are formed or not. Such an investigation would 


be very interesting from the point of view o X-ray optics, but has so 
far not been undertaken. 
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The above considerations show that the X-ray evidence on the multi- 
phase mixtures is not quite conclusive. However, the evidence of the 
unit cell dimensions and the arrangement of the paraffin chains confirms 
the multi-phase character of these mixtures, and the intensities observed 
are consistent with it. This appears to establish a strong case for the 
conclusion drawn from the melting point data, that the stable mixtures 
of up to four ketones consist of aggregates of the individual ketones, 
each forming dipolar planes. 

As for the single-phase mixtures, X-ray. photographs were obtained 
of stable mixtures of five and six ketones and of metastable mixtures of 
four ketones ; metastable mixtures of three and two ketones were not 
photographed, because of their short life. All these mixtures gave 
X-ray patterns corresponding to a structure with an hexagonal arrange- 
ment of chains in the basal planeab. The structure of the pure ketones 
is orthorhombic (or alternatively monoclinic with a monoclinic angle 
of very nearly 90°) and the chains in the basal plane a b are arranged in 
a distorted hexagon (see fig. 1(a)). For the single-phase mixtures the 
distortion of this hexagon vanishes so that a=1/3(b). 

The unit cell of the single-phase form may be referred to hexagonal 
indices, when the reflections 110 and 200 of the orthorhombic notation 
coalesce into one reflection, with indices 100 in the hexagonal notation. This 
reflection which is shown on Pl. XX.(a), and another one, corresponding 
to the coalescing of the orthorhombic 310 and 020 into the hexagonal 
210 were the only hk0 reflections which could be distinguished on the 
powder photographs. Even the second line, the hexagonal 210, was 
so weak that monochromatic radiation had to be used for its detection. 
This indicates strong thermal vibrations at right angles to. the chains. 
However, the reflections present establish a hexagonal packing of the 
chains in which each is surrounded by six nearest neighbours—the 
closest packing for cylindrical rods. 

The 00/ reflections of the single-phase mixtures correspond to a ¢ axis 
equal to that of the pure ketones and multi-phase mixtures (see Table IT.a). 
This implies that the paraffin chains are at right angles to the basal plane, 
as in the pure ketones. The intensities of the 001 reflections of any 
single phase mixture of four ketones are closely similar to those of the 
multi-phase mixture of the same four ketones. As has been pointed 
out above, the intensities of these reflections cannot be used to determine 
whether the ketogroups are distributed at random or arranged in dipolar 
planes. It may be noted in this context that the high orders of the 001 
reflections with /=18—20 remain visible on the pattern of the single-phase 
mixtures. These reflections correspond approximately to the so-called 
backbone spacing of the paraffin chain, 1-27A. The presence of these 
reflections indicates that thermal vibrations are less violent in the direction 
of the chains than normal to them. 

The hexagonal structure of the single-phase mixtures is more symmetrical 
than that of the pure ketones and the volume of its unit cell is about 
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3 per cent larger than for a pure ketone. The hexagonal form thus 
represents a disordered packing of paraffin chains. It is analogous to 
the hexagonal high temperature form which occurs in some paraffins 
(Miiller 1932). Miiller explained this structure as being due to a twisting 
and turning of the chains which makes their cross section statistically 
circular instead of flat. It is interesting that the single-phase form of the 
mixed ketones is hexagonal over a wide temperature range, while the 
paraffin of the same chainlength melts without becoming hexagonal. 
This may be due to the presence of ketogroups at different levels. The 
dipole moment of these groups implies that they repel or attract each 
other, according to their configuration, and this is likely to cause some 
twisting of the chains even in their state of lowest potential energy. 

On the whole, the X-ray data of the mixtures of equivalent ketones 
are in good agreement with the theory. They show that the single-phase 
form is more symmetrical and thus less ordered than the multiphase 
form, being hexagonal, while the multiphase form is orthorhombic or 
monoclinic. The multiphase form has a structure which is consistent 
with a physical mixture of the pure ketones. However, it has to be borne 
in mind that the X-ray data could not establish the presence or absence 
of dipolar planes. It would be desirable to refine the X-ray investigations 
of these mixtures, using single crystals. . 


OBSERVATIONS UNDER THE MICROSCOPE. 


Although the mixtures of equivalent. ketones were not available as 
single crystals, for the purpose of the X-ray work, they could be made 
to grow crystallites distinct enough for investigation under the microscope. 
The technique used was as follows*:—a small quantity of mixture, of © 
the order 5-10 milligrams, is put on a glass slide and covered with a 
coverslip. The mixture is then molten and the weight of the coverslip 
makes it spread in a thin film. If the specimen is now cooled slowly, 
only one or a few crystals form in the direction of the thickness. The 
crystals tend to be oriented with their ¢ axes normal to the glass plates, 
the more so the slower the cooling. In any specimen it is generally 
possible to find some crystals filling most of the thickness and with 
their ¢ axes normal to the plates. 

Observation in convergent, polarized light, shows that all single phase 
mixtures are optically uniaxial, the ¢ axis being the unique axis. This 
is in agreement with the hexagonal structure deduced from the X-ray 
diffractions pattern. The multi-phase mixtures are optically biaxial, 

‘with a small axial angle, and no significant difference can be found 
between these mixtures and the pure ketones. This agrees with their 
X-ray diffraction pattern. 


a ee 


* The author wishes to thank Mr. Scott Harley of the Burmah Oil Company, 
Ltd., for pointing out this method. 
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In parallel polarized light metastable mixtures of four ketones may 
be observed to change into the stable form. The hexagonal form presents 
a picture in which dark areas—crystals with their unique ¢ axes normal 
to the plate—alternate with coloured crystals of other orientations (see 
Pl. XXI.(a)). When the hexagonal form changes into the orthorhombic 
form this picture changes. A pattern appears within the dark areas, 
which consists of bright birefringent needles (see Pi. XXI.(b)). These 
needles seem to represent minute cracks, more or less in crystallographic 
directions. The contraction manifested by these cracks is in agreement 
with the fact that the volume of the hexagonal form is 3 per cent larger 
than that of the monoclinic form, 


A Low TEMPERATURE TRANSITION IN THE SINGLE-PHASE 
MIXTURES. 


When single-phase mixtures of ketones were observed under the 
microscope, it was found that at about 15° C. they change into a form 
very similar to the multi-phase form (see Pl. XXI.(c)). This change was 
observed in detail ‘on a microscope stage which could be cooled, using 
solid CO, (details of this microscope stage will be described elsewhere). 
It was found that on cooling from room temperature birefringent cracks 
began to appear at about 14-5-15-5°C. and became more intense as 
the temperature fell, until their intensity reached saturation round 
13-5-14° C. These cracks make a pattern very similar to that formed 
when a single phase mixture of four ketones changes into its multi-phase 
form. The only difference is that now the cracks are longer and more 
distinct. On heating, the cracks vanish again, but at a slightly higher 
temperature ; they begin to close at about 14-5-15-5° C. and the last 
cracks disappear at about 16-17°C. The transformation is rapid, it 
follows the change of temperature of the microscope stage without 
observable delay. 

The low temperature transition occurs equally in all single phase 
mixtures of four, five and six ketones, at the same temperatures. The 
crystal structure of this form was examined by X-rays using a cooling 
attachment to a “‘ Unicam ”’ single-crystal goniometer. This attachment 
will be described elsewhere. Pl. XX. (c) shows the X-ray powder pattern 
of the low temperature form. It is indistinguishable from the pattern 
of the multi-phase form, also with regard to the higher order lines not 
shown on the Plate. When the multi-phase form is cooled below 14° GC. 
it remains unchanged, both with regard to its X-ray structure and its 
appearance under the microscope. 

As for the physical significance of the low temperature transition, it 
might be thought at first sight that it corresponds to a segregation of the 
ketones into the multi-phase form, for equation (8) shows that at sufficiently 
low temperatures it is always the multi-phase form which is stable. If 
a mixture of five or six ketones is, in equilibrium, single-phase at its 
melting point, it must, in equilibrium, break up into the n-phase form 
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below a temperature given by equation (8). However, the low temperature 
transition cannot correspond to this effect, for two reasons. Firstly, it 
occurs at the same temperature, whether the mixture concerned consists: 
of four, five or six ketones. This contradicts equation (8). Secondly, the 
transition is very rapid. This is in contradiction with the observation 
that it takes several days for a mixture of four ketones to change into 
the stable form, at 18-20°C. At 15°C. the time for the change should 
be still longer. The rapidity of the low temperature transition implies. 
that it cannot be due to a lengthy process of diffusion. 

As the low temperature transition is not due to a segregation of 
different ketones, it must be due to a contraction of the hexagonal form, 
the ketogroups remaining distributed at random. The existence of such 
a transition is not very surprising because it is known that pure paraffins 
show transitions of this kind from an almost hexagonal high temperature 
form to a monoclinic or orthorhombic low temperature form (Miiller 
1932). A mixture of ketones without dipolar planes should not be 
very different from a paraffin if we assume that the chief contribution 
of the ketogroup to the binding of the structure is their dipolar interaction. 
Indeed, the single phase mixtures of five and six equivalent ketones have 
melting points of 27-5 and 27-0 respectively, as compared with 28° C. 
for the paraffin C,,H,,. However, the comparison is complicated by 
the fact that the paraffin C,,H;, has a structure in which the chains are 
inclined to the basal plane, while they are normal to it for the mixtures 
of ketones. 

Table II.b shows the dimensions of the unit cell for the low temperature 
form of the single-phase mixtures. These agree with the cell dimensions 
of the multi-phase form, within the limits of error. As far as the 
intensities are concerned the powder patterns of the single-phase low 
temperature form and the multi-phase form are indistinguishable. This 
implies that the packing of the chains is very similar. 

The information provided by the X-ray data makes it possible to 
derive an expression for the free energy of the low temperature single-phase 
form. The free energy of the hexagonal form (equation (5)) contains an 
entropy of mixing R log 2n. This should be the same for the low 
temperature form, if the ketogroups remain distributed at random 
between the possible levels. However, there should be no further difference 
between the entropy of the low temperature form and that of the pure 
ketone, because the volume and arrangement of the paraffin chains 
is the same for both. This should imply an equal mobility of the paraffin 
chains in the first approximation. Thus the entropy term for the low 
temperature form is S)+R log 2n. As for the total energy, it should 
differ from that of a pure ketone only in that dipolar planes are not 
present. ‘Thus the total energy of the low temperature form should be 


H,+A/p, and the free energy 


G,=Hj-thp—T(Sp+R log 2m), - - ss (9) 


+ 
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The transition from the low temperature form to the hexagonal form 
may in principle be a transition of the first or second order. The 
experimental evidence indicates that it occurs over a temperature interval 
of 1-2°C. and is thus rather sharp. (The difference between the 
temperature of the transition or heating and cooling appears to be a 
hysteresis effect). Thus it seems justifiable to treat the transition as 
one of the first order. If this is done the transition temperature, Ts, 
can be derived by equating the free energies of the low temperature and 
hexagonal forms, which gives 


red WE Ta A ad oe Da i | 


THe Enereay or A DreoLaR PLANE. 

Comparison of equation (10) with with equation (8) shows that the two 
equations contain three unknown variables, namely h,, hp and s. 
Attempts have been made to find a third equation containing these 
variables, so that they could be determined, but so far without success. 
However, equations (8) and (10) allow and approximate calculation of hp, 
provided some rough estimate may be found for s, because 


Nye lely log 20-- (i — 13 )8. .or | ul ee a 


Here n=4-6, T, is the transition temperature between multi- and single- 
phase mixtures for n=4-6, namely 301° K, and Ty, is the transition 
temperature between hexagonal and low temperature forms, namely 
288° K. As T,—T; is small, hp is little sensitive to the value of s. 

It is easy to find an upper limit for s, because it must certainly be less 
than the entropy of melting for a pure ketone. This has been measured 
by Oldham and Ubbelohde (1939) for the series of equivalent ketones 
of chainlength seventeen, and found to be 41-7 cal/mole degree on the 
average (see Table I.). This value for s gives hj)<1310+550 cal/mole. 
However, the entropy of melting is certain to be higher than the 
entropy of a transition in the solid, and it is possible to find a closer 
approximation for s from the data of Ubbelohde (1938) or Garner, 
van Bibber and King (1931) who measured specific heats and latent heats 
of transitions for different paraffins. Ubbelohde’s data, which have been 
evaluated by Frohlich (1949) give a change of entropy, connected with the 
transition, of about 6 cal/mole degree. Garner’s data give an entropy of 
transition of about 21 cal/mole degree. Assuming a value intermediate 
between the two alternatives gives 


hp=1500 cal/mole 
which should be a fair first approximation, 


The experimental value for hy may be compared with the energy of 


electrostatic interaction for a dipolar plane, which has been calculated 
by Froéhlich (1946) to be 
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where , is the dipole moment of the C=0 group, r the distance between 
dipoles in the dipolar plane and <,, the dielectric constant of the medium 
-between them. ‘The latter is rather an indefinite quantity because the 
dipoles are not far enough apart to permit the material between them 
to be considered as a continuum. This means that «,, might have 
any value between unity and the non-polar part of the dielectric 
constant of a ketone, namely 2—2-5. Evaluation of U/2 with e, 
indeterminate gives 
U 2080 


2 Ew 


cal/mole, 


and this may be compared with the experimental value hyp=U/2=1500 
cal/mole. The two quantities become equal if «,=1-4, which seems a 
very plausible value. Thus the experimentally determined energy of 
the dipolar planes can be accounted for on the basis of dipolar interaction. 


CONCLUSIONS. 


Ketones of the formula C,,H,,,,,COC Ho 41 with p23, q23 and pq 
were considered as equivalent ketones throughout the present paper. 
‘This corresponds to the assumption that all the bonding forces of these 
molecules are equal except for the presence of a ketogroup at different 
parts of the paraffin chain. The ketogroup was assumed to influence 
the bonding only by means ofits dipole moment. With these assumptions, 
the treatment of mixtures of equal amounts of equivalent ketones could 
be based exclusively on statistical mechanics and classical electrostatics. 
‘This model was on the whole surprisingly successful. It was found 
that, as far as the packing of the paraffin chains is concerned, the 
mixtures always have one of two structures, of which one is disordered 
compared with the other one. The conditions of stability of these 
structures could be derived from the theoretical model, and were found 
to depend on the configuration of the ketogroups. Two kinds of order— 
disorder transitions were observed, one corresponding to the closing 
of a multiple solid-solubility gap and the other analogous to an order— 
disorder transition ‘in a paraffin. 

In considering the results obtained it has to be borne in mind that 
the crystallographic data available are somewhat crude and need refining 
by means of single-crystal work. A more detailed crystallographic 
investigation would presumably reveal minor differences between mixtures, 
the structures of which were classed as identical in the present paper. 
Such an investigation might prove interesting from the point of view 
of X-ray optics. Apart from this, an X-ray study of the transition from 
‘a metastable single-phase to a stable multi-phase mixture might provide 
interesting data on internal diffusion. 

The data obtained in the present work on the conditions of the order— 
disorder transitions permitted the calculation of the electrostatic energy 
.of a dipolar plane, in first approximation, but not the calculation of the 
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other entropies and energies involved in the transitions. It would be 
interesting to determine those, for instance by the measurement of specific 
heats. Such measurement would also provide a check on the validity 
of the assumptions made. 

Apart from their intrinsic interest the mixtures of oauivelent ketones. 
may provide a tool for studying the dependence of physical properties 
on structural features in long-chain materials. The different equivalent 
ketones on the one hand and the various mixtures of ketones on the 
other hand represent a number of closely related structures, in which 
structural features vary in different respects. This may prove useful, 
for instance, for the study of molecular vibrations. 


ACKNOWLEDGMENTS. 


The author wishes to thank Mr. C. G. Garton for valuable suggestions, 
Dr. J. W. H. Oldham for the loan of samples of ketones and advice on 
chemical synthesis, and Professor H. Frohlich and Dr. H. N. V: Temperley 
for helpful discussions on theoretical aspects of this paper. The author 
is indebted to the Director of the British Electrical and Allied Industries 
Research Association for his encouragement and for permission to publish 
this paper. 

REFERENCES. 


Frouticn, H., 1946, Proc. Roy. Soc. A, 185, 399 ; 1949, Theory of Dielectrics 
(Oxford), p. 148. 

Fow ier, R., and GuecEnHerm, E. A., 1939, Statistical Thermodynamics 
(Cambridge). 

GARNER, W. E., vaN Bipper, K., and Kine, A. M., 1931, J. Chem. Soc., 1533. 

Jackson, W., 1935 a, Proc. Roy. Soc. A, 150, 197; 1935 b, Ibid., 153, 158. 

McArruvr, I., 1944, Proc. Leeds Phil. Soc. 4, 170. 

Mtuusr, A., 1928, Proc. Roy. Soc. A, 120, 437; 1930, [bid, 127, 417; 1932, 
Ibid., 188, 514 ; 1937, [bid., 158, 403 ; 1938, Lbid., 166, 316. 

OLpHaAM, J. W. Lae and UBBELOHDE, A.R. , 1939, Trans. Faraday Soc., 35, 328 ; 
1940, Proc. Roy. Soc. A, 176, 50. 

Prrer, S. Ho CuIBNaLL, A. C., Hopkins, 8. J., Pottarp, A., Surrn, J. A. B., 
and WILLIAMS, 1Dig lite 1931, Biochem. J., 25, 2087. 

Prtmorg#, D. R., 1939, Proc. Roy. Soc. A, 172, ‘502. 

SILLARS, R. W., 1939, ’Proc, Roy. Soc. A, 469, 66. 

UBBELOHDE, A. R., 1938, Trans. Faraday Soc., 34, 228. 


| VERA DANIEL Phil. Mag. Ser, 7, Vol. 41, Pl. XX. 


(a) 


C,+C,+C,+C, 
metastable, 18° C. 


(6) 
€,+C;,+C,+C, 
stable, 18° C. 


(c) 
C,+ C; + C,+ Ce 
metastable, 11° C. 
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(a) 
C,+C,+C,+C, 
metastable, r.T. 


(b) 
C,+C,+C,+C, 
stable, 7.T. 


(¢) 
C,+C, +C,+C, 


metastable, 12° C. 
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§ 1. INTRODUCTION. 


It was shown recently by the author in a letter to the Editor of the Physical 
Review (Thirring 1950 a), that the quantization of higher order equations 
leads to results which correspond to the new subtraction method called 
regularization. But by applying the usual formalism of quantized fields to 
higher order equations one ought in general to interpret the field as 
belonging to particles which possess states of several rest masses. One is 
therefore confronted with the question whether one should accept a 
realistic or a formalistic interpretation. This question arises in the 
regularization method as well as in the field theories of Bhabha, Podolsky, 
etc. It is possible to circumvent this problem by using transcendental 
operators which possess only one zero-point in momentum space, as was 
done by Born and Green, but one encounters here the obstacle that a 
transcendental operator becomes hopelessly difficult, if one tries to intro- 
duce generalized derivatives in a space with an electromagnetic field. 
Another way of avoiding the problem of different rest masses is to take 
field equations of the form 


(TA p2)"$(a)=0. 


In this paper we shall discuss the next simplest case after the Klein— 
Gordon equation, namely 
(C1?—p?)(L?—p?)o(2)=9. 
The time independent Green-functuon of this equation is the exponential 
potential 
exp(—ur) 1 8 exp (—pr) 

© is aly 7) ae ae 
As this potential has the advantage over the Yukawa potential that it 
makes no r-* singularity in the tensor force, one is inclined to take this 
equation for the meson rather than the Klein—Gordon equation. By 
taking a reasonable meson mass and a simple scalar coupling we can 
adjust our coupling constant to well known nuclear data, and after having 
so fixed every constant we may proceed to calculate the self energy of the 
nucleon due to this meson field. It turns out, that the second order self 
energy is less than 8 per cent of the rest energy of the nucleus, but that 
in higher order approximations there are still divergent terms. The 
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reason is that one can make the meson self energy due to virtual nucleon 
creation finite only by taking another nucleon equation but not by chang- 
ing the meson equation. The goal of this paper is, of course, not to give a 
final solution of the divergence problem, but to show that by taking higher 
order equations some renormalization terms, which used to be infinite, 
are now finite and have an actual small numerical value. In this paper 
we use the new formalism of quantized fields, developed by Schwinger and 
the calculation method of Dyson—Feynman. 


§ 2. FoRMALISM. 


First we shall define our notation, which is very similar to the one used 
by Feynman (1949). We choose natural units A=c=1 and take real 
world coordinates with the metric tensor %y=1, 941:=—922=933= — 1. 


———s 
The scalar product of four-vectors we write ab =a,b)—ab, an arrow 
indicates a spatial vector. The four-dimensional volume is written with 


out index Uk adloyd k= dhe dhe, Hh, dhs, By 


= wemean 3, —x3—, — _ 
Ox, Of, | On, Ore one 


the d’Alembertian is 


ad a 
Ces pat aS Ne ea 
L'= 02, an, 1" 


The surface-element of a space-like surface is 
do,—=—dx,dagd23, dxydx.dx3, dxydx,dx,, dx, dx, day. 
The nucleons are represented by a quantized Dirac field 


0 “= 
((V—m)b(x)=0, Vay, =BE +B a grad. 
we 


The y’s obey the commutation rules y,y,+y,y,=29,,.The scalar product 
of a four-vector with the matrix vector we write a=a,y,,. 


The adjoint spinor is given by p= p*ys, wl 7 Mi REL lee 0. In momentum 
space the equations are : 
#=U, exp (—tpz), (p—m)U,=0, U,(p-+-m)=0, pP=m 
For the vacuum expectation value we have [Dyson (1949), Feynmen (1949)] 
P<p(x), Pu") >p=he(w—x' (iV +-m)Ap(e—2") 


r» 20 ¢ ptm 
Oe mT) Dae 


| 


te(a— 


me xP {—tp(a—2"')} dp 
X~ 


[mol 


Ed) — 
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P is the operator which arranges factors chronologically. We shall only 
consider free nucleons with sharp energy-momentum, and get therefore 
for the one particle expectation value 


P<H(x), p(x") >,=«(e—w')U,U, exp {ip(e—a’)}. 


Coming now to the meson field we first carry out quantization of the 
vacuum field and enter then immediately into the interaction representa- 
tion for the field interacting with the nucleons. From here we can then 
deduce easily the field equations in the Heisenberg representation. 
Our vacuum equation is 


(eS ea His) =O ahah wis Metheny iene 151) 


The index # indicates the component of the field function in isotopic spin 
space. Hquation (1) is derived from a Lagrangian L 


1k, & ad, 84, 
L=— 5 (sa baa 0x, On, a utr) Soe iy Bes (2) 


The commutation relations for a field with ne order Lagrangian are 
given by Weiss (1938), Podolsky and Kikuchi (1944), Chang (1946), 
Chang (1948), De Wet wee Le Couteur (1948) 


[5 — : aes Ey oap \AV— iw egt see ava ( a) 


OL {a’ } 0 OL {a’} Pas: 
| Eo Dad,(a’y xp) dx; OOP, lx" )/Oark dar , |av ST eae) 


In our case, if we replace the space integration by the integration over 
an arbitrary space-like surface, we get 


[[S&. C4,(0") |dow=— iB, pare an 
7) Ona a") 
[ [4 2 gible — gz C780) doi... (40) 
This gives for the commutator function 
[Pi(x), bx(’)]=O¢D(w—a’)  . we 5) 
PPO) dot, [ap Ple) do,=0. Patera 


a is here any space-like surface through the origin. D(«) has to be a 
solution of (1), therefore the associated function 


D(z)=4e(z) Diz), . - - . - we C7) 


obeys (1) everywhere but at the origin. To investigate its behaviour 
there we consider the integral 


b del HF Deo), eS ee) 
k 
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where & is any region round the origin. ye Gauss’ theorem (8) 
becomes 


ta} a aes 
[Sige (C20) + | deeDe) 


o is the surface of k. 
If we make the region infinitesimally thin and choose the surface to be 


space-like we get 
0 
8 =i done B—-2us Diet 
) . Vag fs 


D(a) is therefore the invariant Green-function of (1), the same connection 
with opposite sign holds between the Green-function and the function 
in the case of the second order field. The invariant Green-function of 
(1) is 

hae. Ilo exp =the yg Oe 
D, is the Green-function of the Klein-Gordon equation. Though it is 
evident from the integral-representation one can easily prove that (9) 
is the right Green-function of (1) 


5 (8) = (8), 0) + (29D, (0) 


=A5(t)e— bs + ay ee a 2 eee 


(C2 —p2)2D(x)= 


For the commutator function we obtain 


—~ lee 
D(x)=2€(x%)D(x)= — om Bene ia 2h cee 
D,(«) is the D-function for the second order field with rest mass » and as 
is necessary, (11) is an odd invariant solution of (1). For the calculation 
we need Feynman’s function which is 


5, 


P<Hy(o\by(a)>o= Dayle) eS (SPI) 


2 Gxt) eS 
(5a) 
One notices the similarity between (5a) and Feynman’s cut-off function. 
Passing from the vacuum fields to interacting fields we keep our field 
equations, but the state vector is now no longer a constant but varies 
according to the fundamental equation of the interaction representation 


6¥(c . 
Sa BW) ¥(C). Tee ay eae) 


For the interaction energy #(x) we take a simple scalar coupling 


(x)= —p/(2) \orpy(%) blac) rahi CD OMROG De Deere ol (13) 
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The factor is introduced to make the coupling constant g, dimension- 
less, 7, are the isotopic spin operators. The field equations in the 
Heisenberg representation are easily available. If we indicate the field 
quantities in the Heisenberg representation for a moment by Clarendon 
type, the general connection between the derivatives of the field quantitites 
in the Heisenberg representation and the interaction representation is. 
given by Schwinger (1948) 

we, =U) FO) p60) 410-4) | [H(x), F(x’)] do} U(c). (14) 


Ox, ox, 


U is the unitary operator which performs the transformation between 
the Heisenberg and the interaction representations. Using (6), (13), 
(14) we obtain 


Op; 1 94; ee it a2 


a) ) (15) 
mene tee (J) ete 
Ox, L GP; U Ox, LJ dU, 
CPP = UAL PLP GU + Ue 2)geblw) 7H ()U, 
therefore the field equations for the meson field in the Heisenberg repre- 
sentation are : <*- 
(Pwr =guv(2prap. - 2 2 + (16) 
To prove the third equation of (15) we need 
—ld 5 
2 ee Aldara. 2) do,=0. 
iL CPD(2) do,= <= = | 2D, (x) do, =0. 
For the spinor field we get [Schwinger (1948)] if we use 
SB), (yy@ble'))p} dop=8,9 (A, B}=AB+BA, 
and the identity . 


[Plac)yb(ae), abla’) =x) fb(x), Wa")}— (a), P(x’) (x) 


the following field-equations in the Heisenberg representation 


GV—m+gprV/(2)7p)P=0. . - - 2 = + (17) 
The field quantities in the Heisenberg representation are easily expressed 
in terms of the field quantities in the interaction representation [Schwinger 
(1949 b), Dyson (1949)]. For the interaction energy in the Heisenberg 
representation we get 
0o ,\n 
Hg, ae aa) 


n=0 n!} 


[7 dey. [dey PQ (aa)AB(e2) «-- WUeq))- (18) 


First we consider the expectation value of the second term of the series 
(17) for a state where two nucleons and one meson are present, which is 
nothing but twice* the second order interaction between two nucleons. 
eee ee er 


* Half of it is subtracted if we transform Hp. 
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We get, if we integrate over the space volume, 
2E i= JH2(%) 4 % 
= — 12g pg fd xo da ,P <rp(aq) 7a (a%o)(@1)7 ph (%1) Pr(Lo) bx(%1) > 2, 0+ 


For simplicity we neglect, in evaluating the expectation value, exchange 
terms which give rise to mixed densities. 


: 1 —lé re 
Em = 20 soma Ge Ga | __ Lele tle) + pilaa)et(er) 


cp {—tk(a,—x = 
x eh, iavdkes na cata LS) 
The upper index on the nucleon densities pj(w)= <ah(a)7,4h(x) > indicates 
to which nucleon the density belongs. To obtain the static energy we 


assume the densities to be time independent, carry out the time inte- 
gration and get 


2 => — a | 
Buy =— de | p(eo)pP( ey) exp {—wley—ay)} dity de. E> (20) 


Or if we assume the densities to be sufficiently concentrated and approxi- 
mate them by a 6-function 


2 
Biye=—p Ze 1} exp{—p(e1—@)}. Pee rare Oa 


x, and x, are the positions of the nucleons. The negative sign shows that 
we get attractive forces. The same result yields a classical meson theory 
if we start from equations (16) and (13) and replace the nucleon by a 
point source. We are now able to determine the coupling constant g. 
We may either choose symmetric theory g=9,=92=93, J1=0, or neutral 
theory g=94, 91=92=93=0. In both cases pg?/47 is the strength of the 
potential. To fit the proton—proton scattering data for an exponential 
potential we need a strength of about 100 MeV. if we choose a meson mass 
of 282 m (Rosenfeld 1948). This gives a coupling constant g?/47=0-71. 
It is clear that the expontial potential wants a stronger coupling than the 
Yukawa potential with its singularity. Now everything is ready for 
working out the self energy. 


§ 3. CALCULATION OF THE SECOND ORDER SELF ENERGY. 


The second order self energy is given by the expectation value of the 


second term in (17) for a state where one nucleon and no meson is present. 
The expectation value 


<jp(y)p (wp) P(x )Wlwy)h (eo) h(@1) >, 0 
is of course 


<4 (Xo xf (aig xs (224 xb (1) > 1 <(a9) (21) >o- 
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The one-particle part is Schwinger (1949 a) 


<p (94h (9) WP ary) > = <P(®o bey) >1 <P (@o)H(@1) >o 
co <4p(aq)(21) > 1 <p(ap)eb(ay)>o 
which makes the self energy 


Eyer = J <Ho(%) >1,0 4% 


= — 1279.94) diy dx, <thgT hs > 1 P<$((%)$x(%1) >o 
x [P <p, (2o)4(21) >1P <p, (x1) p(ao) >o 
+ P<, (27 )Wal%o)>1P <Py(Mo)ba(1)>ol ~~ (22) 


Inserting (5a) and the expectation values for the spinor field and taking 
into account that the vacuum expectation value of two spinors gives 
Kronecker operators in spin and isotopic spin spaces we get 


Eyap= —i4y?g? | day dx, dp’ dx €(%p—2)U,p exp {ip(%)»—2,)} 


+ p+m —i exp {—tk(%—%}) 
X (to 81) BH Tne aa 


Ai Tie 1 
(23) 


It is easy to see that the two different terms in the one particle 
expectation value give rise to the same integral (23) with different sign 
of k. The sign can be restored by a change of integration variable, we 
may therefore replace the second term by a factor two in the first one. 
Further, the square of every isotopic matrix 7 is the unit matrix, all kinds 
of mesons give therefore the same contribution to the self energy, which is 
equal for neutron and proton. We will now drop the index on the 
coupling constant and calculate the contribution to the self energy of one 
kind of meson. The integration over space and nucleon momentum in 
(23) yields then 
p—k+m 


Pag = 
= —j U0, a os Ue Hw + (4 
i Bayi | Um (Pa (P— ph) 9? a (24) 


We have integrated over the 2) unit volume, because we shall renormalize 
our spinor function to the unit volume. To carry out the operation in 
spin space which is required in (24) we observe that it is compatible with 
the Dirac equation to choose the vector U,(p) in spin space so that 
U,U;=N(p+) 5. As desired this corresponds to an unpolarized wave 
since Uo,,U=0. The normalization constant N we may take to normalize 
either to one mesonic charge per unit volume U,U,=1, N=1/4m or to 
one electric charge per unit volume 02,01, N=4p,/1. For conveni- 
ence we take the former possibility. 
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After an easy spur computation the integral (24) is then evaluated 
by a straightforward calculation using Feynman’s method (Feynman 1949) 
and we obtain 


—9u? josie 494 a mn (Y i aa 


Kye 2(2a)2m? 0 (1 —x)* + ap4/m? Pre (477)? m 0 (~—%)(%—2p) 
(25) 


The zeros of the denominator in (25) are worked out to be 


Dn — op? 2 
#1 »=exp (+iy), i Peal S00] a= ~ = =0-0118, y=0-153. 


? 


l—a 
(26) 
Introducing the further abbreviations 
@y9—1L=—atirv/(2a—a«?)=P exp (418), P=0-1535, PB—=1-647, 
(27) 
the elementary integration yields 
aes : {sin 2y . In P+(B—y) cos 2y}=6-3 
SCS) nu ahnge emt WO SERS aa 
Hy iL 
———_———- dx = —— {sin y ln P+ (B—y) cos y}=7-7. 
I, (*—2)(x—2°) sin ue te +(B—y) y} 
Inserting into (25) we get 
—q2 214- 2 
Be! (E) wn OEE Gh 0:0 16-1 eee eet io) 
7 \m) 47 7 


The negative sign corresponds to the fact that we have attractive forces. 
For pseudoscalar coupling the self energy is even less, namely 


2 


B.uecpsy= tt! 0-00454. 


The second order self energy gives therefore only a small contribution to 
the observable mass of the nucleon. 

It might be interesting to calculate the classical self energy. If we 
replace the nucleon density in (16) for the static case 


(4 sn? )*6 (a) == ple) Fae te nO) 


by a point source 


p(a)=guv/(2)8(a) 


and use the above mentioned solution of (16a) we get for the self energy 


1 /(2) ex Ve 
Bete (et) = — 5 [ ple x)p(a) He = 5 { mev/(2) (2 5(a ye ) ee 
_ =P, , ,, dedk  —g pp ' —9 
(27)! { 2778(ko) (2 =e Ta om — year O768m. . (30) 
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Comparing the integral representation of the classical self energy given 
in (28) with (24) we see that the self energy of the quantized field in (24) 
corresponds to a spread out of the nucleon density. This causes the 
reinforcement of the denominator in (24) which cuts off the contributions 
of the high frequencies and makes the quantum theoretical self energy 
a bit smaller than the classical one. 


§4. HIGHER ORDER APPROXIMATIONS. 


In the fourth order contribution to the self energy there are terms which 
are even smaller than the second order self energy and which correspond 


to double emission of mesons. The Feynman graphs of these terms 
would be 


Fig. 1. Fig. 2. 


but there are still divergent terms arising from the creation of nucleons 
as in the following graph: 

Calculating the contribution to the self energy due to this graph we 
encounter the following integral : 

(p—k-+m)tr(p’+m)(p'—k—™) 2 ay, 
| a9 ww 0) 

Though the denominator is of degree 12 and the numerator is of degree 11, 
this integral is divergent as it contains the self energy of the meson 
which has the divergent integral 


pipe ie EM) 
J pee Em) 


dp’. 
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On this integral the form of the meson equation has no influence and it is 
infinite in our theory as well as in the usual one. The same fact mars 
the other higher order terms. 

In conclusion we may state, that it is possible to make finite terms 
which are infinite in the usual theory, by taking higher order equations, 
but it does not yet seem to be possible to select from the manifold of 
linear equations of arbitrary order just those equations which might be 
realized in nature. 


I should like to take this opportunity of thanking Professor Schrédinger 
for valuable discussions in connection with this paper. 


Note added in proof.—In the meantime the same equation was proposed 
independently by Bhabha (Phys. Rev., 77, 665, 1950). For a closer 
discussion of the difficulties of this theory see Pais and Uhlenbeck (to 
appear in Phys. Rev., 77, lst July, 1950). 


Fig. 3. 
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SUMMARY. 

The hypothesis that the change in the c/a ratio in magnesium solid 
solutions is due to the changes in the Fermi energy, arising from the 
distortion of the Brillouin zone, is developed into a more general and 
quantitative theory. The case of monovalent solutes is considered, and. 
agreement with experiment is shown to be not unsatisfactory. The 
experimental results enable deductions to be made regarding the density 
of electronic states in pure magnesium. 


§ 1. INTRODUCTION. 


Tue changes which occur in the lattice parameters a and ¢ of magnesium 
solid solutions, as the concentration of the solute is varied, have been 
measured with great accuracy by Hume-Rothery and Raynor (1940) and 
Raynor (1942). A simple qualitative theory has already been proposed to 
account for the variation of the c/a ratio in those solid solutions for which 
the valency of the solute is greater than that of magnesium (Jones 1949). 
The experimental results show unambiguously, in all cases, that this c/a 
ratio is determined, not by the atomic concentration of the solute, but by 
the valence electron concentration in the alloy. 

It is the purpose of this paper to express the theory in a more general and. 
quantitative form than hitherto, and to extend it to include those cases 
where the valency of the solute is less than the valency of magnesium, 
e.g. to the cases of lithium and silver. 

The changes brought about in the overall form of the lattice by the 
addition of other elements in solid solutions may be regarded as the sum of 
two effects, (i) a change of atomic volume, and (ii) the development of a 
shearing strain. The change in the atomic volume must depend upon 
many factors, e.g. the difference between the size of the solute ion and the 
size of the magnesium ion, and also upon the difference in the potential 
fields of the solute and solvent ions in which the valence electrons move. 
The complexity of the causes of the change of volume can clearly be seen 
in the experimental observations which show no obvious regularities 
between the dilatation at a given concentration and the solute metal. 


* Communicated by the Author. 
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Since the structure of pure magnesium is almost exactly that of hexagon- 
ally close packed spheres (the c/a ratio for Mg is 1-6237 and the value for 
ideal packing 1-6330), and since the metal is comparatively elastically 
isotropic, there is no reason to believe that differences of ionic size, or of 
spherically symmetrical potential fields, will play a significant part in 
producing the shearing strain. It might, therefore, be expected @ priori 
that the shearing strain would result mainly from the effect of the Brillouin- 
zone structure on the valence electrons of the metal. This supposition is 
confirmed by the fact, already mentioned, that the shearing strain is 
substantially independent of the particular solute metal and depends, to a 
large degree, only on the ratio of the number of electrons to atoms in the 
metal. 


§ 2. A CONVENIENT FORMULATION OF THE EXPERIMENTAL RESULTS. 


It may be assumed that the relations between the lattice parameters 
and the concentrations of solute atoms will not be significantly different at 
the absolute zero of temperature from the observed relations at room 
temperature. Thus it will be sufficient to express the condition for 
equilibrium by minimizing the energy rather than the free energy. 

The energy per atom W, expressed relative to the energy of the separated 
atoms, is a function of the atomic volume », of the shearing strain, and, in 
the first place, of the solute concentration x. It is convenient to introduce 
the variable s=(c/a)?* so that a=(4v/1/3)1%s-12 and c—(40/1/3)"8s. The 
conditions for equilibrium can then be written : 


(i) eo) AY ise ok (ee ato WE aay 


In order to obtain the variation of c/a with x, as determined experimentally, 
we should have to eliminate v between these two equations. We shall 
make the assumption that the changes of volume arise from forces which 
act equally in all directions and can thus be simulated by the effect of a 
uniform pressure or tension. Such a uniform pressure or tension, in an 
anisotropic elastic crystal, produces a shearing strain given by 


813 +S33—811—8}9 


where | a 
2 (811 +819) +48 13835 


and the s;; are the elastic coefficients of the crystal. Thus, if we subtract 
this strain from the observed value, we get the strain induced by the solute 
atoms at constant volume. This quantity which we can determine from 
the experimental results can now be compared with the relation between 
s and x given by the second of equations (1) alone. The reason for this 
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procedure is that the dependence of W on v is too complicated to be amen- 
able to calculation, whilst the dependence of s on x at constant v, we 


believe, arises from a particular effect which can be comparatively easily 
isolated. 


Fig. 1. 
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(a) The variation of the c/a ratio with electron concentration in Mg solid solutions. 
(b) Detail of (a) for very small concentrations. 


If Sa/a and Sce/c are the observed strains, we get, for the change of s at 


constant volume, 
8—8\'_, foe 8a\_ ag, oc oa 
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For pure magnesium k=0-0226. In fig. 1 {(s—sp)/so}, is plotted, not 
against atomic concentration, but against the increase of the number of 
electrons per atom in the alloy which we denote by. If /” is the number 
of valence electrons per atom of the solute, n=(”—2)a. The values of 
AN assumed for the solute metals are as follows : 


Li Ag Al Ga In Tl Sn Pb 
N= k 1 3 3 3 3 4 4 


§3. A THEORY OF THE ORIGIN OF THE SHEARING STRAIN. 

Since fig. 1 shows clearly that variations of s at constant volume depend 
only on n, we shall use, from this point onwards, n as the variable in place 
of x; and also since all changes, now to be considered, take place at constant 
volume, there is no further need to make explicit reference to v, and there- 
fore this quantity will no longer be shown as a variable. The condition 
for equilibrium becomes therefore (@W(s, )/ds),=0, which is the required 
equation to determine s as a function of n. For the derivative of s with 
respect to n, we have therefore 


ds @W //@W 
FPS = saan | ele . . . . . . (5) 


Since (W/v) is equal to the strain energy, apart from a constant, and since 
the strain energy in terms of s can easily be shown to be $(s—89)?/s9°, 
where 

P=$(C11+€12+2C33— 43), 


and c¢,; are the elastic constants of the crystal, we find 


Ledom So\ 07W 
1 EN (e) et ee 


It will now be assumed that the only part of W which contributes to the 
second derivative on the right-hand side of (6) is the Fermi energy. This 
Fermi energy per atom, which we denote by w,, may be considered as 
made up of several parts as follows : 


Wp=- Wp paw paw, ge ee ee tee beer) 


where wy is the energy corresponding to the completely filled Brillouin 
zone of two electrons per atom. In magnesium or its alloys some electrons 
must overlap into the second zone, at various points, leaving vacancies in 
the first zone. The overlapping electrons and the holes may lie in distinct 
regions of momentum space. The energy of one such region containing 
n; overlapping electrons is denoted by w; and the energy of n, electrons 
which would just fill one particular region of holes is written w; Thus 
we see that w, is given by (7), where the summations are over all distinct 
regions of holes (j) and overlapping electrons (i). 
w,is a function of s and n; only, but ,is itself a function of s and n, and 
therefore . 


(= _ (ew; dw;\ (dn; 
08. a me), + Ge), (Ge). , 7 ; ; (8) 


on the Lattice Spacings in Magnesium Solid Solutions 667 


Now (dw,/0n;),=E,(s, n;); where E; is the energy at the top of the filled 
levels in the region (7). For equilibrium, we have 


Ky=h,= fe: estrone i = Pe ah ey 1 =E,=6, . ° . . (9) 
where € is the energy at the surface of the filled states, and thus 
on; On; 
zB, (5), SPU, (S2). Oe... at com 
since Te LN Dis. ee a x), ET) 


Thus, since w, is independent of n, we have, using (8) and (10), 


wy  ,, Ow, (On; aw, (On, 
ade 2 ast oe), 2 aaaes (ae) ‘ayant ae 


If N(¢) denotes the total density of states at the energy ¢, and N,(Z) the 
contribution from the region of momentum space (7), then clearly 


(FS) =e eee ty toe, eer 3) 


with exactly similar relations, apart from a negative sign, holding for the 
regions of the holes. | 
Thus (12) becomes 


Oty al C1 OF 0K; \ 
Fa WH EN Fe), +58, ( = ), ee ei 
Finally, if we define coefficients ¢; by the equations 
Sy (OE; 
b= 2() : ee eee) 
(6) and (14) give 
i ds C 1 
a7 (=) xq BANE 24N}) Ee te 06) 


§ 4. COMPARISON WITH EXPERIMENTAL ReEsULts. 


In order to apply equation (16) to the magnesium solid solutions, we 
have to make detailed assumptions with regard to where the holes and the 
overlapping electrons occur in relation to the Brillouin zone. Fig. 2 (a) 
shows the form of the zone which contains, for all values of s, exactly two 
electrons per atom. We shall assume that the holes in the first zone occur 
at points along the line D,D, and regard the centre of this line D as the 
point of highest energy *. We shall also assume that electrons overlap at 
the points A, and that for n>0 overlap may also occur at the points B. 

All points like those indicated by D in fig. 2 (a) are equivalent and lie 
at equal distances from the origin for all values of s. The neighbour- 
hoods of D form altogether six identical regions in which the surfaces of 


nr 

* Until calculations of the regions of highest energy within the Brillouin zone 
have been made any assumption such as this is necessarily arbitrary since the 
geometry of the zone is not sufficient to indicate the relative energies. It is of 
some interest that the coefficient dp depends upon the assumed position, of the 
point of highest energy and that the experimental data relating to the variation 
of c/a are most conveniently explained by the present choice. 
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constant energy are simply connected closed surfaces. There are two 
separate regions corresponding to the points A, and one to the points B. 
The distances from the origin, k,, kp and kp, are given by : 


ky = (1/3/40)? 380, 


ba=(V/3/40)"9 5, 


3 
by=(v/3/40)!® 5 (14+ zo) 3 


In order to calculate the coefficients ¢;, we notice, first, that 0E,/ds is taken 
at constant ,; and, therefore, the changes in E; must be closely equal to 
the changes in E; when n;=0, 7. e. of the energy of the centre point of the 
region, and secondly, in equilibrium the value of E; is always equal to ¢. 


Fig. 2. 
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(a) Brillouin zone for Mg. (6) Density of states for overlapping 


electrons N, and for holes Np. 


For these reasons, we may assume that d,, for example, will be given 
without serious error by d log E,/d logs. It is to be noticed that the ¢’s 
are unlikely to depend at all sensitively on n, and, of course, according to 
the above assumption, they are completely independent of n. For the 
application of (16), this is their most important property since it means that 
the right-hand side of this equation only varies with n significantly through 
the densities of states N;. Numerical estimates of the ¢’s may be obtained 
by assuming that E,, E, and Ep vary as the squares of the distances from 
the origin to the points A, B and D respectively. In this way we find : 


18 13 
—— Is =— 2, = See 
Pa bp bp 3-453 4] . 
Thus (16) may be written 


1 ds Enrl 78 
a (2) es ose Oe 
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where N,, Nz and Ny denote the densities of states of single regions about 
their respective points. 

Since, even for n=—0, there is overlap at A and holes at D, we can expand 
N, and Np about their values at n=0. Thus, for instance, 


dN di 
N,(2)=NY+ (Zz), (=), 


BEN at ee Fhe ee Si, oe yes (18) 


where N is the total density of states in pure magnesium. Ny cannot, of 
course, be expanded in this way, since it is zero up to a certain (very small) 
value of » which we denote by vn). For n>n), we may assume that 


Np=«/(f—Eg), 


where « is a constant. Also we have n—n,.=(C—E,)N, and therefore 


N,= Guerin ang 7 ak earaais = (1-O) 


TN (n—n 


Finally, substituting for N,, Np and Np into equation (17), and integrating 
with respect to n, we find 


io) =A(n—n,)??+Bn+Cn?, ... . . . (20) 
0 v 
aN rd. =) a : 
Shee \enae if n>, 
where Aes @ NS ‘ ee + (21) 
0 ; if N<Np, 
2 (39.0) wo) (§ 2 
B x (GiN8 A (= : (22) 
1 /39:dNp dN,\/¢ 
chat as ma TA og EAN ald aa 
and | C N2 & dt wt) (= Che 


The experimental results, for solutes with >2, indicate (cf. fig. 1) that 
B must be zero or exceedingly small. Hence we deduce that NY and Nf 
must be very nearly equal. 

It is clear from equation (23) that C must be negative, for dN)/d¢ is 
necessarily negative and dN,/d¢ positive (see ¢.9. fig. 2(6)). Thus 
monovalent solutes, for which n<0, must decrease the c/a ratio. This 
conclusion is in agreement with the experimental results for lithium and | 
silver. 

It is possible to make a direct estimate of the number C by assuming that 
N, and Ny, have approximately parabolic form as shown in fig. 2 (6). If 
c, is the energy width of the electrons overlapping at A, and «, the energy 
width of the distribution of holes, then 

dN, Na INyn Np 
pices oa Te, 
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For magnesium (=7-3 eV. as given by the soft X-ray emission spectra 
(Skinner 1939), p=6-23x101dynes/em.2, and v=23-08x 10-4 cm.? 
Thus the pure number (¢/uv)=0-813. Hence, neglecting the difference 
between 39/41 and unity, we obtain 
o_ 07813 (Na | No 
Rr OMNINGG hes Gere 
and, since N=6N,+2N,, and N,=Np, we get, if we assume that N may 
be taken as 3/C, which is the free electron value and probably slightly too. 


big, 
0:813¢ /1 li 
Oncaea cacy) 


If we take e,=e,—1 eV. which is roughly the value of «, suggested by the 
X-ray data, we find C=—0-247. The value of C which has been used to 
plot the curves of fig. 1 is —0-213, and is thus quite compatible with the 
above estimate. The curves of fig. 1 are given by 


(=) = —0-213n?+0-150(n—0-005)??,  — n>0-005, 
0 v 


= —0:213n? : n<0:005, 


where it has been assumed that overlap sets in at the points B when n 
exceeds 0-5 per cent. 

The agreement is not good in the case of silver, but the curve for 
negative 7 goes exactly through the only known point for the magnesium 
lithium solid solution, viz. that at the phase boundary limit of 17 atomic 
per cent lithium *. At this concentration the observed (s—sg)/sp is equal 
to —:0067. 

It must be admitted that the interpretation of the experimental results, 
which we have adopted, depends largely on the observations at very small 
n. Apart from this region, a fairly good straight line could be drawn 
through the whole set of points both for positive and negativen. The 
interpretation would then be that the constants A and C were relatively 
small, and that B was positive, implying that Nj>>N,. To distinguish 
with certainty between these possibilities more data relating to the strain 
produced by monovalent solutes would be necessary. 
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SUMMARY 


Plastic anisotropy in rolled sheet is examined by tension and cupping 
tests. The earing positions are correlated with the strain-ratios measured 
in tension tests at various orientations to the direction of rolling. The 
materials used are copper giving four ears at 45°, brass giving four ears. 
at 50°, and brass giving six ears at 0° and 60°. A theory of plastic 
anisotropy due to Hill is found to be in good agreement with the 
experimental data for materials producing four ears. The theory is 
extended to describe more complex states of anisotropy. 


§1. MeTHoD OF CoRRELATION, 


Tue deformation imparted during the rolling of sheet metal produces an 
elongation of the crystal grains and a tendency towards a preferred 
orientation. The resulting anisotropy can be reduced or accentuated at 
will by a suitable programme of rolling reductions and interstage heat- 
treatment. Directionality in the final product is revealed by earing in 
a cupping test on a circular blank, or by the variation of mechanical 
properties in tension tests on specimens cut from the sheet at different 
angles to the direction of rolling. It is evident that the plastic behaviour 
of a metal with crystallographic anisotropy is directly dependent on the 
plastic properties of a single crystal of the metal and on the degree of 
preferred orientation. If the mechanism of distortion of a grain in an 
aggregate were completely known, it would be possible, in principle, to 
predict the behaviour of the aggregate when loaded in any way. This 
cannot yet be done since it is not known, for instance, how the operative 
glide system in an individual grain is influenced by the constraint of 
neighbouring grains. For example, in a face-centred cubic metal there 
are twelve operable slip-directions at ordinary temperatures, whereas a 
general uniform strain can be produced by a combination of shears in 
only five independent directions. Sufficient conditions for the deter- 
mination of a unique active system are unknown. For this reason alone, 
preliminary attempts to relate the earing positions to the rolling texture 
seem to be premature ft. 


* Communicated by the Authors. 

+ For a general review of past work see T. Ll. Richards, “ Progress in Metal 
Physics”, 1 (1949), pp. 281-305. (Butterworth’s Scientific Publications, 
London.) 
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A more practical alternative is to measure the macroscopic properties 
of the anisotropic sheet in simple tests and, from the data so obtained, 
attempt to predict the behaviour of the sheet when deformed by any 
complex stress system. The question immediately arises as to which 
properties of the sheet are relevant to its behaviour under combined 
stress, in particular when deep-drawn into a cup. It is obvious that 
only mechanical properties are relevant, and that of these the breaking 
stress and elongation can be ignored when we are not concerned with 
fracture. Again, when we are interested only in large plastic deformations, 
the elastic moduli are irrelevant except in so far as they affect the residual 
stresses after unloading. There remains only the plastic stress-strain 
properties, which incorporate (i) the dependence of work-hardening on 
the amount of deformation, and (ii) the relations between the ratios of 
the components of the stress and the plastic strain-increment. In 
general, both factors govern the distortion of the sheet when taken along 
some complex loading-path. In the cupping test, however, the 
conditions are fairly simple, each element of the rim being subjected 
to a circumferential strain of approximately the same amount if the 
material is not too anisotropic. Moreover, ear-formation presumably 
depends mainly (though not entirely) on conditions near the rim, and 
here the state of stress is essentially a pure circumferential compression. 
The blank-holder pressure is usually negligible compared with the yield 
stress, and hinders the thickening of the edge only slightly because of 
the elasticity of the apparatus. Hence factor (ii) appears likely to 
predominate in controlling the positions of the ears. The relevant 
fundamental test is therefore the measurement of strain-ratios in uniaxial 
compression at various angles in the plane of the sheet. For convenience 
it is obviously preferable to perform tension tests instead, although this 
involves the assumption that the strain-ratios are unaffected by micro- 
scopic internal stresses which exert different influences in tension and 
compression. This seems entirely reasonable since the final annealing 
should remove those stresses responsible for the Bauschinger effect and 
allied phenomena. The measurement of strain-ratios in tension tests is 
the basis of the present method of predicting earing positions. In essence, 
the method has been used previously by Baldwin, Howald and Ross (1945), 
but the details differ, as will be explained later. 


§2. THEORETICAL CONSIDERATIONS. 

According to a theory of plastic anisotropy proposed by Hill (1948), 
the relations between the ratios of the components of the stress and 
strain-increment are 

at de, A de, dy... 

(G-fHjc,=Hoe, 7 (F+-H)o=Hos Ne? 7 et 
where the x-axis is the direction of rolling and the y-axis is transverse 
to this in the plane of the sheet. The applied stress is considered to be 
in the plane of the sheet, with components g,,, a,, and 7,,. V5; Giitiaand 
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N are parameters specifying the current state of anisotropy, which will 
normally vary during a process of plastic distortion. Similar linear 
relations have been independently derived, though from different physical 
assumptions, by Jackson, Smith and Lankford (1949) and by Dorn (1949). 
Dorn’s relations contain five parameters, and cannot be derived from a 
plastic potential (see equation (10) below). 

(i) Strain-ratios in a tension test.—For a uniaxial tension o applied at an 
angle « to the direction of rolling, the (x, y) components of stress are 


= 2 , 
Tz=TCOS* 4, Dy 


=o sin? a, Toy G SIN  COB.C. 1 2 (4) 
The strain-increment in the width direction of the tension specimen, 


with current width w, is 


dw : ya : 
aa =de, sin? «+de, cos? «—2dy,,sinacosa, . . + (3) 


while the strain-increment through the thickness ¢ is 
dt 
re = — (de, +de,) . . : . . . . ° (4) 
since there is no change of volume. Combining these four equations 
we obtain 
dwjw  H+(2N—F—G—4H) sin? « cos? « 
dir F sin? «+G cos? « : 
If the state of anisotropy is not altered by the deformation during the test, 
the ratios of the parameters remain constant and (5) can be integrated 
to give 
In'(w ow)  H+(2N—F—G— 4H) sin? « cos? « ins 
In (t/t) F sin? «-+G cos? « "5 
where w, and ft, are the initial width and thickness. The ratio r(«) 
between the logarithmic strains in the width and thickness directions 
then has a constant value throughout any one test. In a completely 
isotropic material F=G=H=N/3 and r(«)= 1. If the sheet is isotropic 
only in its plane, N=F+2H=—G-+ 2H, 7(«)= H/F. In a non-isotropic 
material the ratios H/G and H/F ean be found from tension tests along 
and perpendicular to the direction of rolling. Thus : 


Ey a =r(/2). eer pee eet (7) 


(5) 


(a), say, . (6) 


A third ratio, N/H say, can be obtained from a test in some other direction. 
The theoretical relation (6) is then completely known, and its validity can 
be checked by comparing it with the experimental results over the whole 
range of orientations. 

Some general properties of (6) are worth noting, Stationary values of 
r(x) occur when «=0°, 90°, and also for a certain intermediate value f 
except when (2N—F—G—4H)/H lies between (F--G)/G and (G—F)/F. 


— 
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The intermediate stationary value is a maximum when (2N—F—G—4H)/H 
is positive, and a minimum when it is negative. In particular, when the 
properties are the same in both the rolling and transverse directions, so 
that F=G, the r(«) relation is symmetrical about B=45°; its value at 
the latter is a maximum if N>F+2H and a minimum if N<F+2H. 

In an anisotropic specimen a tensile stress generally produces a shear 
strain as well as an extension; that is, a principal strain axis does not 
usually coincide with the direction of the applied stress. It is easily 
shown from the above equations that the shear-strain increment, referred 
to axes along and perpendicular to the applied stress, is equal to 


dw [(N—F—2H) sin? «—(N—G—2H,) cos? «] sin « cos « (8) 
ee (H-+-(2N—F—G—4) sin? « cos? a] “ges 


The presence of a shear strain has been demonstrated experimentally by 
Hazlett, Robinson and Dorn (1949). In the present tests the angle of 
shear, as computed from (8) with the deduced values of the parameters, 
was never more than about 5°, even at the largest extensions. The 
corresponding eccentricity of loading was negligible in view of other 
experimental errors and the unavoidable scatter in the results. It can 
be seen from (8) that the axes of stress and strain-increment coincide 
when «=0°, 90°, or an intermediate value y satisfying 


N—G—2H 


tan2 


y coincides with the stationary point f of the r(«) relation only when F=G. 

(ii) Angles of necking.—A sufficiently thin flat strip necks at an oblique 
angle when pulled in tension. Ina perfectly isotropic specimen the angle 

is about 55° to the axis. According to theory (Hill 1948), the line of 
the neck develops in the direction of zero rate of extension in the uniformly 
deforming strip. In a non-isotropic specimen the two directions of zero 
extension-rate make different angles with the axis, except when «—0°, 
90°, or y. If, then, the neck does in fact coincide with a direction of 
zero extension, a further check on the preceding theory is provided by 
measurements of the necking angles. For example, for‘a test in the 
rolling direction the theoretical inclination of the neck to the axis is 
tan-! 4/(1+G/H). 

Unfortunately the fractures in the copper and brass specimens used in 
the present experiments were too ragged for sufficiently precise measure- 
ments of the necking-angles. It appeared, however, that the necks 
usually formed at greater angles to the specimen axis than would be 
expected theoretically. This discrepancy is probably due to a slight 
necking in the width direction, which occurred towards the end of the 
considerable range of uniform strain and extended over a large part of 
the gauge length. The theory seems to give closer agreement in a metal 
where the range of uniform strain is comparatively small (Hill 1948). 


ae 
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(iii) Haring positions.—The state of stress and strain in a circular blank 
which is being deep-drawn into a cup can, in principle, be calculated from 
equations (1) along with some assumed law of work-hardening under 
combined stress. The positions and gradual development of the ears 
would be included in the results of the analysis. However, although the 
calculation can be made for an isotropic blank (Hill 1949) the analysis 
becomes extremely involved when the blank is anisotropic. Consequently, 
some simpler method must be adopted for predicting the earing positions. 

In view of the approximate symmetry of the deformation about the 
tip of an ear or the bottom of a hollow, it seems a reasonable assumption 
(Hill 1948) that the ears and hollows develop at the points on the rim where 
the principal axes of the stress and the strain-increment coincide. Since 
the stress in the rim is a circumferential compression, it follows from the 
previous theory that the ears and hollows should begin to form at points 
on the rim whose angular positions are 0°, 90° and 90°—y, relative to the 
direction of rolling. There are therefore four ears and four hollows (at 
most), and the equations (1) must be modified to represent the anisotropy 
of a material giving more than four ears (see below). Now, when the 
yield criterion and stress-strain equations are connected on the plastic 
potential basis, Hill has shown that the initial yield stress has stationary 
values in the directions for which the axes of stress and strain coincide. 
It is natural to assume, therefore, that the ears form in the positions 
where the (circumferential )yield stress is initially a minimum (numerically) 
and the hollows where it isa maximum. Then, if N is greater than both 
F+2H and G+2H the ears are in the 90°—y positions, and if N is less 
than both F+2H and G+2H the ears are in the 0° and 90° positions. 
If N is intermediate to F+2H and G+2H only two ears are formed ; 
these are in the 0° positions if F>G and in the 90° positions if F<G. 

In particular, if F=G, there are four ears in the 45° positions ifN>F-+2H, 
and in the 0° and 90° positions if N<F+2H. These are also the positions 
where the strain-ratio 7 is a maximum in the special case F=G. This, 
indeed, is what would be expected if the circumferential strain in the 
blank did not depend on the angular coordinate, since the ears correspond 
to the parts of the blank which suffer most radial extension. When 
applied generally, this consideration leads to the hypothesis that the 
ears develop at positions where the strain-ratio r is a maximum for a 
uniaxial stress in the circumferential direction. However, this alternative 
hypothesis is unlikely to be correct since in actuality the circumferential 
strain varies considerably around the rim of a non-isotropic blank. In 
fact, the experimental evidence reviewed later is more in accord with 
the assumption that the ears form at positions where the strain-ratio r is 
a maximum for a uniaxial stress in the radial direction. This third 
hypothesis was suggested on empirical grounds by Baldwin Howald and 
Ross, but it is difficult to see a reasonable theoretical basis for it. 

It is worth remarking that the absence of ears in a cupping test does not 
prove that. the sheet is completely isotropic, but only isotropic in its 
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plane. According to the above theory this implies N=F+2H=G-+2H, 
but H is not necessarily equal to F and G. To prove complete isotropy 
it is essential to perform, in addition, a tensile test in the plane of the 
sheet to see whether the strain-ratio is unity or not. 

With suitable rolling and heat treatment some metals can be made to 
give six ears in a cupping test, e.g. brass, or even eight ears, ¢. g. certain 
aluminium alloys (Chevigny 1949). To represent these more complex 
states of anistropy it is necessary to generalize the stress-strain 
equations (1). This may be achieved (Hill 1950) by writing 


dey bey) Mags go 9h ae a 


Of/dc,  afldc,  af/or,, ; 
where f, the plastic potential, is taken to be a polynomial of degree n 
in the components of stress*. f is homogeneous if the strain-ratios do 
not depend on the absolute magnitude of the stress. The corresponding 
number of ears is 2n (or less, in special circumstances). The plastic 
potential for equations (1) is the quadratic 


(G+H)o?—2Ho,0,+(H+F)o?+2N77,,. SRR me ys) of (i tet 


To represent a material giving six ears, f must be a cubic (at least) and 
for a material giving eight ears f must be a quartic (at least). 


§3. MATERIALS AND TREATMENTS. 


The materials used were copper and brass. In pure copper there are 
two common recrystallization textures, viz. (100) [001] giving four ears 
at 0° and 90°, and (113) [211] giving four ears at 45°. Only the latter 
was studied here since the former has been thoroughly investigated by 
Baldwin, Howald and Ross. In 70/30 brass the usual earing positions 
are at 45° or at 0° and 60°. Both these were examined. 

The copper was phosphorus deoxidized, with an average grain size of 
006mm. It was received as a hot-rolled and pickled strip of width 7 in. 
and thickness in. ‘The strip was machined to 0-210 in. to remove surface 
pitting, and rolled to a nominal thickness of 0-030 in. (85 per cent 
reduction). It was then annealed at 650° C. for half an hour. 

The brass was received as a hot-rolled strip of width 7 in. and thickness 
gin. To obtain a texture giving four ears the strip was rolled 50 per cent, 
annealed at 425°C., rolled 84 per cent, and annealed at 650°C. The 
final nominal thickness was 0-030 in. To obtain a texture giving six 
ears the strip was rolled 50 per cent, annealed at 650° C., rolled 84 per cent, 
and annealed at 650° C. The final thickness was again 0-030 in. Burghoff 
and Bohlen (1942) have determined the (111) pole figures for both treat- 
ments. 


a 

* In the differentiation 7, and 7, are treated as distinct, the terms containing 
the shear-stress component being regarded as split into two equal parts. Only 
even powers of the shear-stress are allowable since the axes of anisotropy are 
the axes of reference. 
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In the cupping tests 4 in. diameter blanks were drawn with a punch of 
diameter 2 in. over a die of radius 0-187 in. The nominal blank-holder 
clearance was 0-003 in. Thelubricant was a mixture of tallow and 
graphite. Cups drawn from the copper strip had four large ears at 45° 
to the rolling direction. The four-eared brass cups had large ears at 
approximately 50° to the direction of rolling, and the six-eared brass 
cups had small ears at 0° and 58° (+2°). 

The tensile test pieces were produced in a blanking press and milled 
to size, the edges being polished on emery paper. Specimens (in 
duplicate) were cut at 0°, 223°, 45°, 67$° and 90°. The gauge length 
was 2 in. and the width 0-5in. The tests were carried out ina Hounsfield 
tensometer. Extensions in the elastic and initial plastic range were 
measured by a Lindley extensometer. Width and thickness measure- 
ments (average of three along the gauge length) were taken with a 
hand-micrometer reading to 0-0001 in. The readings were continued up 
to the stage at which necking began. 


§ 4, EXPERIMENTAL RESULTS. 


With all the materials tested it was found that the strain-ratio 7 did 
not vary with strain, within experimental scatter. The constancy of r 
was maintained until necking intervened ; this occurred at extensions of 
the order of 30-40 per cent. The inference to be drawn is that the 
initial anisotropy in the rolled strip was not appreciably altered by the 
uniform tensile strains. 

Values of r(x), derived from the slopes of the best straight lines through 
the experimental points (means of at least two tests), are shown in fig. 1. 
Variations in the material, errors of measurement, and scatter in the 
results, combine to produce an uncertainty of order +0-1 in the derived 
values of r. Ratios of F, G, H, N were obtained by fitting (6) to the 
results for the copper and brass giving four ears. There is some 
uncertainty in the values of the parameters, amounting perhaps to 
+0-1, owing to the arbitrariness in choosing the best fit. 

The curve in fig. 1(a) for the copper corresponds to F=G=H=N/4. 
The theory in § 2 (iii) predicts four ears at 45° when F=G and N>F-+2H, 
in agreement with experiment. Comparison may be made with the 
data of Baldwin, Howald and Ross, for copper giving 45° ears (see the 
lower curve in fig. 18 of their paper). The parameters for their material, 
which had received a different treatment, are F=G—=H=N/5. 

The parameters for the brass giving four ears, corresponding to the 
curve in fig. 1(b), are 

Mae HN le Ga 1-2 35-5. 


The angle y in (9) then has the value 40°, approximately, so that 
90°—y=50°. Since N is greater than both F+2H and G+ 2H, the theory 
predicts four ears at 50°. This is what was observed. It is, of course, 
hardly possible to be precise about the actual angular position of the tip 
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of an ear, especially since there are always slight differences between each 
of the four ears. In view of the uncertainty in the derived values of the 
parameters, a discrepancy of less than 5° between theory and experiment 
would probably not be significant. 

As a further check on the validity of equation (1), a compression test 
normal to the sheet was carried out on a cylindrical brass disc of diameter 
0-25 in. and thickness 0-029 in. The section became elliptical, with the 
major axis in the direction of rolling. After a 50 per cent compression 
the axes were 0-38 in. and 0-325 in. (approximately, since the edge became 
ragged)*. The corresponding ratio of the logarithmic strains is 1-6. 


Fig. 1. 
r@ r@ 
O22 4S 67% 90q0 O22. 45: if O78 90,6 
(a) Copper with four ears (6) Brass with four ears 
at 45°. at 50° 


r(c0) 
2 


(e) Ss 4 4 4 —t 
224 45 674 90.0 224 45 674 S0q° 
(c) Brass with six ears - (d) Theoretical variation of 


at 0° and 58°. yield stress with orienta- 
tion in brass with six ears. 


Regarding the test as equivalent to the application of equal biaxial 
tensions in the plane of the sheet, equation (1) gives de,/de,=G/F=1-6, 
in perfect agreement. 
In fig. 1(c), for the brass giving six ears at 0° and 58°, the curve 
corresponds to the relation 
oe: 6 cos® «+6 sin® «+38 sin* « cos? «—10 sin? « cost « 
3 cost «-+-6 sint x-+sin2 « cos? « 2 
This was obtained by assuming a homogeneous cubic polynomial for the 
plastic potential defined in equation (10), and choosing the ratios of the 
six numerical coefficients so that the following five conditions were 


* Calcium oleate was the lubricant (coefficient of friction ~0-015). With 


these specimen dimensions the ellipticity would not be significantly influenced 
by friction. 
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satisfied : (i) r=2 in the rolling direction ; (ii) r=2 in the 45° direction ; 
(iii) y=1 in the transverse direction ; (iv) the axes of stress and strain 
to coincide in a tension test at 30° to the direction of rolling; (v) the 
ratio of the strain in the rolling direction to the transverse strain to be 
1-5 in a compression test normal to the sheet. The latter is regarded as 
equivalent to equal biaxial tension in the plane of the sheet. The value 
1-5 is near to that observed in a compression test on a cylindrical disc of 
_ diameter 0-25 in. and height 0-03 in., lubricated with calcium oleate. 
After a 50 per cent compression the axes of the ellipse were 0-375 in. and 
‘0-325 in., approximately. 
The plastic potential so determined has the form 


f= 303— 6070, —60,07+ 403+ (40,42lo,)72, . . . (12) 


According to the hypothesis of §2 (iii) the fourth condition ensures that 
either an ear or a hollow occupies the 60° position. The other possible 
positions of ears and hollows are calculated to be 0°, 90°, and cot~14/11=17° 
approximately. It will be seen that the calculated curve passes close to 
the two experimental points not deliberately fitted, viz. at 224° and 67}°. 
This method of comparison was preferred to the obvious one of fitting the 
curve to all five points, and then comparing predicted and observed 
earing positions, because there is considerable latitude in fitting a cubic 
curve to so few points when they are uncertain to +0-1. 

It remains to distinguish between the calculated positions of the ears 
and hollows. If the yield criterion appropriate to (12) is f=constant 
(depending on the strain-history), the ears are predicted at 0° and 60° if 
it is supposed that they form at points where the circumferential yield 
stress is a (relative) maximum. This is opposite to what was found above 
for copper and brass sheet giving four ears, where the ears form at points 
where the circumferential yield stress is, theoretically, an (absolute) 
minimum. We are unable to suggest an explanation of the difference ; 
it may be an accidental product of the theory, depending in a complex 
way on the entire distribution of the stresses and strains in the initial 
stages of the cupping test. The calculated variation of the yield stress 
with the orientation of the tension specimen is shown in fig. 1 (d). 

It is interesting to notice that the earing positions are close to the angles 
for which the strain-ratio r has maximum values for a radial tension. 
This is the hypothesis of Baldwin, Howald and Ross, referred to pre- 
viously. The present experiments do not indicate whether this is the 
real factor determining earing, or whether it is merely that the maximum 
values of r happen to fall rather close to the actual earing positions on 
account of some more fundamental property (which may, indeed, be 
peculiar to brass and copper). 

The variation of the strain-ratio r with orientation has been measured 
for an aluminium alloy (61S-T6) by Hazlett, Robinson and Dorn (1949). 
The composition was 0-25 per cent Cu, 0-6 per cent Si, 1:0 per cent Mg, 
and 0-25 per cent Cr. The alloy was solution heat-treated, quenched, 
and aged. Values of the parameters giving best fit of equation (6) to 
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their data are F:G:H:N=1:1-15: 0-65: 2-25. Since N is only 
slightly different from F+2H or G+2H the shear strains (8) should be 
negligible, and this was found to be so. The earing positions were not 
measured. The same authors found that the strain-ratios for a 
magnesium alloy with only slight anisotropy varied during a tension test. 
Measurements of strain-ratios in an aluminium alloy, for tests both in 
_ and out of the plane of the plate, have been reported by Klingler and 
Sachs (1948). 

Strain-ratios for low earhou steel have been measured by Lankford, 
Snyder and Bauscher (1950), who make the interesting suggestion that 
plastic anisotropy may sometimes be advantageous to good performance in 
an unsymmetrical pressing operation. 


§ 5. STRESS-STRAIN CURVES. 


For the correlation of stress-strain curves, obtained in tension tests at 
various angles to the direction of rolling, it is necessary to relate some 
measure of the amount of work-hardening to the previous strain-history. 
Assuming that the amount of hardening can be specified by a scalar 
quantity (which appears probable when the state of anisotropy does not 
vary), the natural measure to take is the value. of the expression f, 
representing the yield criterion and plastic potential. Jackson, Smith 
and Lankford (1949) claim to have found evidence that the degree of 
hardening in anisotropic sheet steel is a function only of the amount of 
plastic work per unit volume (a hypothesis which is a reasonable approxi- 
mation for some isotropic metals.) The present tests do not support 
this hypothesis in the cases of copper and_ brass. 

Thus, the stress-strain curves at 0°, 224°, 45°, 674°, and 90° for the 
brass giving four ears did not vary from a mean by more than 0:8 tn/in.? 
at any strain up to 50 per cent extension. In view of the variation 
between nominally similar specimens, there was scarcely any significant 
difference between the stress-strain curves in the various directions. 
Hence, if the hardening depended only on the amount of plastic work, 
the tensile yield stress should not vary significantly with direction in 
a given mechanical state. This clearly follows from the hypothesis 
that the same mechanical state can be produced by performing the 
same amount of work in tension tests in different directions; since 
the stress-strain curves coincide, the final yield stress is the same in all 
directions. This conclusion, however, is not consistent with the assump- 
tion that the yield criterion is f=constant, since this implies variations 
of the yield stress with direction. For example, according to (11), the 
ratio of the yield stresses in the 0° and 90° directions should be 
/ {(F-+-H)/(G-+H)} or about 0-9; the ratio of the yield stresses in the 
40° and 90° directions is 0-8, approximately. 

A similar conclusion applies also to the brass earing at 0° and 60°, 
where the stress-strain curves showed even less variation and the theoretical 
yield stresses more. 


Properties of Rolled Sheet in Tension and Cupping Tests 681 


The stress-strain curves for the copper, however, showed significant 
differences (fig. 2), especially between the 0° and 90° curves. Since the 
theoretical yield stresses are roughly equal at 0° and 90°, the same 
conclusion is again reached, namely that the hardening (as measured by the 
value of f) is not a unique function of the total plastic work per unit volume. 
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Fig. 2. 
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Stress-strain curves for copper at various inclinations to the 
rolling direction (experimental points omitted for clarity). 
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LXI. The Fluetuation and Fading of Radio Echoes from Meteor Trails. 
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[Plates XXII. & XXIIT.] 


ABSTRACT. 

This paper is concerned with the amplitude fluctuations observed in the 
radio waves reflected from meteor trails. Subsequent to the initial 
diffraction phenomena the radio echoes exhibit fluctuations of period 
between 0-01 and 0-1 seconds, when observed on radio wave frequencies of 
36 to 72 Mc./s. These fluctuations have been investigated simultaneously 
on two radio wave frequencies and it is concluded that the phenomenon is 
due to wind gradients in the high atmosphere which cause the trail of 
meteor ionizatiom to break up. In this way several reflecting regions are 
formed which give an interference pattern moving over the observing 
station. 


§ 1. INTRODUCTION. 


OBSERVATIONS of the radio echoes from meteor trails indicate that the 
amplitude fluctuations during the life of the echo are complex in character. 
For example, during the Giacobinid shower of October 1946, Lovell, 
Banwell and Clegg (1947) investigated the life histories of several echoes 
and found that many showed rapid fluctuations in amplitude superimposed 
on the general decay of the echo. Similar effects have been found 
by McKinley and Millman (1949). These observations were made 
using either cine-photography or photography of intensity modulated 
displays, and it was evident that further progress in the study of this. 
phenomenon could only be made by investigating the amplitude of the 
echo over much shorter time intervals. This paper describes the 
information which has now been obtained on the fluctuations using 
a technique which enables the amplitude of successive radio pulses to be 
measured with a time separation of 6-6 milliseconds. 

During the complete life history of a characteristic radio echo from a 
meteor trail it is possible to distinguish four types of fluctuation. 


(a) Diffraction Phenomena. 


Approximately 10 per cent of meteor echoes show a rapid initial rise in 
amplitude, followed by oscillations of increasing frequency. A typical echo. 
of this type is shown in fig. 1 (Pl. XXII.). The oscillations, which are usually 


* Communicated by Dr. A. C. B. Lovell. 
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complete with 0-5 seconds of the first recorded pulse correspond to the 
passage of the head of the meteor through successive Fresnel diffraction 
zones,.and they can be used to measure the meteor’s velocity (Ellyett and 


Davies 1948, Davies and Ellyett 1949). 


(b) Short Period Flutter. 


Although the majority of echoes have durations of less than 0-2 sec., 
some 25 per cent persist for much longer periods, and show violent 
fluctuations of the type illustrated in fig. 2 (Pl. X-XIII.). These oscillations, 
which have periods of the order 0-01 to 0-1 seconds, are occasionally observed 
directly following the Fresnel fluctuations, but more usually occur with 
echoes where the diffraction phenomena are not observed. 


(c) Long Period Fluctuations. 


In a small percentage of very long duration echoes, a much slower 
amplitude variation, of the order of 1-0 seconds period, has been detected. 


(d) Noise Fluctuations. 


In addition to the foregoing types of amplitude oscillation all echoes 
invariably show small fluctuations between adjacent pulses, which are due 
to random noise in the receiver (Ellyett 1949). 


The phenomena described in (a) and (d) are now well understood, and 
will not be discussed further in this paper. The short period flutter (5) 
and long period fluctuations (c), however, present a number of interesting 
physical problems, and have been investigated in detail. 


§ 2. TECHNIQUE. 


Most of the results have been obtained by a dual wave technique, 
giving continuous pulse by pulse amplitude records on 72 and 36 Mc./s. 
simultaneously, for the whole duration of an echo. Additional information 
has been obtained (on 72 Mc./s. only), from earlier techniques, using 
automatic photography of successive pulse amplitudes from sections of the 
echo. 


(a) Dual wave technique. 

This technique was developed by Ellyett * and a block diagram of the 
equipment is shown in fig. 3. ‘Two transmitting and receiving systems 
were used, operating at frequencies of 72 and 36 Mc./s. Upon receipt of 
the first pulse exceeding twice the signal to noise ratio reflected from a 
meteor trail, a recorder unit (Davies and Ellyett 1949), attached to the 
72 Mc./s. receiver output, actuated the dual wave unit. Both receiver 
outputs were fed continuously via separate amplifying channels in this 
unit, to a double beam cathode ray tube. A camera with continuously 
moving film was switched on simultaneously, so that successive pulses 


ne ee eee NS 
* Dr. C. D. Ellyett returned to New Zealand before the final completion of the 
dual wave technique. 
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were separated on the film, giving a pulse by pulse picture of the echo _ 
amplitude throughout its life. The camera motor required approximately 
0-05 sec. to reach a speed sufficient to separate the individual pulses, so 
this period was lost at the beginning of each echo. Once the recorder had 
actuated the system a hold-on relay came into operation, maintaining the 
unit in action, until switched off manually on conclusion of the echo. 

The two transmitters were locked together to the same pulse repetition — 
frequency, and their peak powers were adjusted so that echoes appearing 
on both wavelengths simultaneously had approximately the same signal to 
noise ratios. The aerials were horizontal folded dipoles with reflectors, 
at heights of half a wavelength above the ground, and looking in a northerly 
direction. Common transmit-receive systems enabled the same aerials 
to be used for both transmitting and receiving. The pulse repetition 
frequency used in these experiments was 150 per sec. 


Fig. 3. 
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Block diagram of the dual wave apparatus. 


(b) Single wave technique. 

An automatic recorder unit designed for observation of the Fresnel 
zone pattern has been described previously (Davies and Ellyett 1949). 
This unit photographs some 100 pulses on one photographic frame, and 
takes 1 sec., to reset. For long duration echoes detailed amplitude 
information can therefore be obtained from short sections of the echo 
separated by intervals of 1 second. 


§ 3. THE SHORT Prertop FLUrrer. 


’ Photographie records of short period flutter haye been obtained during 
most of the major showers of 1948 and 1949, although the dual, wave 
technique. was used only during the 1949 Geminid shower. Examination of 
the results shows that the individual fluctuations which comprise the 
flutter vary continuously in amplitude and frequency during eae lifetime 
of an echo. The results have been analysed to find, 
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(a) the relationship between the echo duration and the mean time 
period of the flutter ; 


(6) the relationship between the mean flutter periods observed on 
frequencies of 36 Mc./s. and 72 Me./s.; and 


(c) the relationship between the flutter period and the type of meteor. 


(a) Relation between Echo duration and Mean Period of Flutter. 


A histogram relating the mean period of the flutter to the echo duration, 
for echoes observed on a frequency of 36 Mc./s., is shown in fig. 4. It is 
clear that the mean period is practically the same for each duration group. 
The processes causing the phenomenon of high speed flutter must therefore 
be independent of the duration of the echo. 


Fig. 4. 


0:25 0:5 1:0 2:0 40 3:0 16-0 32:0 Secs. 
Echo duration. (36 Mc./s.). 
Mean flutter period as a function of echo duration. 


(b) Relationships between the Mean Periods of Fluctuation on 
frequencies of 38 Mc./s. and 72 Mc./s. 
The relationship between the flutter periods on 36 Mc./s. and 72 Mc./s. 
was investigated with the aid of the dual wave apparatus. 
Table I. gives the mean flutter periods for the two frequencies, and the 
ratios between these periods. 


TABLE I. 
oo 2 Eee eee eee 
Mean flutter Ratio of the mean flutter period 
period at 36 Mc./s. to the mean flutter 
(secs.) period at 72 Mc./s. 
Sr nnn ne 
72 Mc./s. echoes — 0-045 1-0 
36 Mc./s. echoes while accom- 
panied by a 72 Mc./s. echo —- 0-082 1-84 
All 36 Mc./s. echoes 0-091 2-02 


pee ee ee 
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It is evident that the period of the oscillations on a frequency of 36 Mc./s. 
is approximately twice that of those observed on 72 Mc./s. In general the 
echo observed on a wavelength of 4 m. is of much shorter duration than the 
echo observed on 8 m..from the same trail (Lovell 1948) but there is no 
significant difference between the mean periods on 36 Mc./s. before and 
after the end of the corresponding 72 Mc./s. echo. 

Since the flutter period is independent of the duration of the echo, the 
probability P(T) of a fluctuation of period T occurring is proportional to 
T.N(T), where N(T) is the total number of fluctuations of this period 
observed. The relationship between T.N(T) and T for the echoes 
photographed on the 72 Mc./s. apparatus is shown in fig. 5. This is 
compared with the corresponding curve for 36 Mc./s. echoes in fig. 6. 


Fig. 5. 


50 - Numberof Pulses: 
— > Seconds 


Flutter period (T). 


P(T)dT: Probability of a fluctuation of period between T and T+dT 
occurring at 72 Mc./s. ({T=0-066 secs). 


If the mean flutter periods T, and T, observed on the 36 Mc./s. and 
72 Mc./s. equipments are related to the wavelengths by a power law, the 
exponent 7 will be given by 

ya logis T,/T, 
301 a 
This exponent was found for 36 echoes observed simultaneously on both 


wavelengths. Measurements were restricted to times when both echoes 
were present, and the distribution of values obtained for » is shown in 
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Flutter period (T). 
(Continuous curve). P(T)dT: probability of a fluctuation of period between T 
and T+4dT occurring at 36 Mc./s. (dT=0-066 secs.) 
(Dotted curve). Corresponding curve for 72 Mc./s. scaled up by a factor of 2. 
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fig. 7. The mean value was found to be 0-99. It therefore seems reason- 
able to conclude that the periods are proportional to the wavelength of the 
radiation employed. 

Thirty-six echoes exceeding 0-5 sec., in duration at 72 Mc./s., whose 
amplitudes remain at a reasonable value for a sufficiently long period, 
have been analysed pulse by pulse. In general, the 72 Mc./s. fluctuations 
are of shorter periods than those at 36 Mc./s., and there are sections in most 
of the echoes where the ratio of the periods is exactly two'to one. Fig. 8 
shows an echo in whicha high degree of correlation persists for an appreciable 
part of the echo life. In this example amplitude maxima at the higher 
frequency correspond approximately to the maxima and minima at 
36 Mc./s. The echo lasted for 1-9 sec. at 72 Mc./s. and 5:8 secs. at 36 Mc./s., 
only the first 2-0 secs. of the 36 Mc./s. echoes are shown. The individual 
pulse amplitudes have been smoothed by a sliding group of two. 
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Duration tn) 
Time: 05.21 hrs. Dec. 13th 1949. 


Pulse repetition frequency 150 ¢/s. 
Duration of 72 Mc./s. echo 1-9 secs. 
Duration of 36 Me./s. echo 5-8 secs. 


Individual pulse amplitudes smoothed by a sliding group of 2. 


Fig. 9 is an echo of shorter duration showing higher speed fluctuations. 
The pulse amplitudes are unsmoothed. Fig. 10 is an example of an echo 
in which the fluctuations are less clearly defined. There is a two to one 
ratio in flutter periods for certain sections of the echo, but in other parts 
the fluctuations are confused and often of small amplitude. In general 
the fluctuations at 72 Mc./s. are faster than those at 36 Mc./s. The first 
3-9 sec., of the echo are shown on both frequencies, although the 36 Me. /s. 
echo persists for a further 14 sec. 


(c) Relation between flutter period and types of meteor. 
In order to discover whether the flutter period is related:to the type of 
meteor by which the trail was originally produced, or if it depends only 
upon the subsequent behaviour of the completed column of ionization, the 
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mean flutter periods for three major showers, and for sporadic meteors 
have been measured on 72 Mc./s. only, using the single wave technique. 
The results are given in Table IT. 


ig. 9: 
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(36 Mel) 
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Duration (Secs) 
Echo 
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(72 Mc/s) 
0 fao:s 
Duration Gecs.) 
Time 22.24 hrs. Dec. 13th 1949. 
Pulse repetition frequency 150 c/s. 
Duration, of 72 Mc./s. echo 0-6 secs. 
Duration of 36 Me./s. echo 2-0 secs. 
TaBueE Il. 
Mean Period of Flutter (72 Mc./s.). 
ce ee Se Ra aie 
No. of Mean period 
Shower observations of flutter 
of period (sec.) 
Perseids (Aug. 1948) 80 0-048 
Geminids (Dec. 1948) 362 0-032 
Quadrantids (Jan. 1949) 74 0-035 
Perseids (Aug. 1949) 353 0-039 
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It is evident that the méan period is independent of the type of meteor 
which produces the initial column of ionization. 


§ 4. Lona PERIopD OSCILLATIONS. 


Some experimental evidence for the existence of a longer period ampli- 
tude fluctuation has been found on 36 Me./s. only. Out of a total of 
30 echoes obtained by Ellyett * during the 1949 Perseid meteor shower, 
two, with durations of 81-0 sec. and 50-0 sec., show long period fluctuations 
emerging after about 30-0 sec., the period of the fluctuations gradually 
increases throughout the remainder of the echo life, and: the mean period 
is of the order of 1-0 sec. 

Similar behaviour has not yet been found using the dual wave equip- 
ment, owing to the shorter durations of the echoes on 72 Mc./s., but 
evidence of waves in one 7:0 sec. echo on 72 Mc./s. has been given by 
Lovell, Banwell and Clegg (1947), using cine-photography. Further 
information on the behaviour of the 72 Mc./s. echoes was obtained from 
a total of 90 cine records of the summer daytime and Perseid meteors of 
1947, which show 10 cases of possible wave occurrence. A detailed 
analysis of the periods of the oscillations cannot be made from these 
records, owing to the possibility of a beating effect with the camera 
shutter mechanism. 

As there are only two positive cases of the occurrence of long period 
amplitude fluctuations at 36 Mc./s., and none at 72 Mc./s., it is not yet 
possible to decide whether this is a separate phenomenon, or merely the 
tail of the distribution curve shown in fig. 4. 


§ 5. THEORETICAL EXPLANATIONS OF THE High SPEED AMPLITUDE 
FLUCTUATIONS. 


Experimental evidence from several sources (for summary see Ellyett 
1950) shows that there are both steady winds, with velocities of the order 
of 50-100 m./sec., and winds of a turbulent character, present in the region 
70-120 km. above the earth. 

Several suggestions have been made that the phenomenon of high speed 
flutter might be explained by the presence of such wind gradients in the 
upper atmosphere. (Herlofson 1948, Ellyett 1949, McKinley and Millman 
1949). The results obtained in this paper will now be considered in terms 
of such wind motions. When a meteor trail is first formed, the ionization 
along the track is approximate linear, and it is known that the reflection of 
radio waves at 72 Mc./s. from the trail is specular (Lovell, Banwell and 
Clegg 1947, Hey and Stewart 1949). If the trail is subsequently deformed 
by winds, then other sections could present right angle reflecting points to 
the observer. An interference pattern due to the interaction of waves 


* Unpublished. 
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reflected from various sections of the trail will then be formed on the sur- 
face of the earth. If the reflecting centres move relatively to each other, 
this interference pattern will move over the ground, producing a series of 
maxima and minima in the echo amplitude at any point. 

If two sections of the trail predominate, and reflect approximately equal 
amounts of energy in the direction of the receiver, large amplitude 
fluctuations will be produced. The frequency of the fluctuations will be 
proportional to the frequency of the radio waves. The maximum possible 
flutter frequency will occur when the reflecting points are moving in 
opposite directions along the line of sight of the observer, with the maxi- 
mum wind velocities present in the upper atmosphere. Any other flutter 
frequency less than this maximum is possible. 

If V is the component of the velocity of the sources relative to each 
other in the direction of the observer, A the wavelength, and T the period 
of the flutter, then 

m\ 

iT: 

V is independent of the wavelength, so that the flutter period T should 
be proportional to A. 

The distribution curve for all flutter periods at 72 Mc./s., scaled up by a 
factor of two is shown together with the corresponding distribution for 
36 Mc./s.in fig. 2.(Pl. X XIIT.). Thereisareasonable correlation between the 
twocurves. The most probable value of V is approximately 25 metres/sec., 
while the highest speed fluctuations correspond to a value of V of the 
order of 100 metres/sec. The slowest fluctuations observed have periods 
of approximately 0-5 sec. at 36 Mc./s. and 0-25 sec: at 72 Mc./s. 

If more than two reflecting centres, of differing amplitudes, and moving 
with different velocities, are present, the fluctuations will become more 
complex, and correlations on widely different frequencies could not be 
expected. The mean flutter period, however, would still be proportional 
to the wavelength. 

Additional information supporting the wind theory has been obtained 
by the visual observation of an auxiliary cathode ray tube. In several 
cases an echo was observed at 36 Mc. /s. a few seconds before the appearance 
of the 72 Mc./s. echo. This suggests strongly that the track is undergoing 
distortion, and that a favourable electron distribution in the trail has 
temporarily caused the 72 Mc./s. echo to appear. 


Vw 
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ABSTRACT. 


The paper contains a discussion of the possible effect of high altitude 
winds on the ionized columns created by meteors. It is shown that the 
various types of range drifts observed in the long duration radio echoes 
from meteor trails can be explained in terms of the influence of winds on 
the ionized column. A comparison of the characteristics of the fading 
of the radio waves returned from the ionosphere, with the characteristics 
of the fluctuations of the radio echoes from meteor trails, leads to the 
conclusion that some of the ionospheric phenomena may be due to meteor 
ionization in the E region. 


$1. INTRODUCTION. 


THERE is now much evidence from the study of the radio wave reflections 
from meteor trails to show that the trail undergoes complex disturbances 
after its formation, and several authors have suggested that winds at 
altitudes of 80-120 km. are an important disturbing factor (Herlofson 
1947, Hey and Stewart 1947, McKinley and Millman 1949). The purpose 
of the present paper is to examine the available evidence and to estimate 
the extent to which winds can account for the observed phenomena. The 
principal phenomena involved are those of range drifts and amplitude 
fluctuations found in the long duration echoes from meteor trails. The 
paper includes a discussion of the correlation of these effects with recent 
work on the fading of radio waves from the ionosphere (Ratcliffe 1948). 


§ 2. WinD Processes at 70-100 km. HEIeur. 


Evidence from the movement of polar aurore, of noctilucent clouds, 
and of meteor dust trails seen by reflected sunlight, gives conclusive proof 
of wind motion at heights of 70-100 km. (Stérmer 1939). The mean drift 
velocity is of the order of 50-100 m./sec. Similar values for the wind 
velocities have been obtained recently by Mitra (1949) from direct 
ionospheric radio measurements. 

In addition to steady winds, turbulent motion also occurs. Olivier (1947), 
analysing the movements of 1492 visual and telescopic meteor trails, finds 
proof of winds with different velocities and directions along a given trail, 
and both vertical and horizontal components within quite small volumes. 


* Communicated by Dr. A. C. B. Lovell. 
+ Now at Canterbury University College, New Zealand. 
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; Based on these observations, a trail of meteor ionization, initially 
linear, may be subject to one or more of the following types of motion :— 
(i) Drift or tilt of the trail as a complete unit ; 
(ii) Gradual curvature of large sections of the trail ; 
(iii) Relative movement of short sections of the trail, due to turbulence. 


These effects will be considered in turn. 


§ 3. INTERPRETATION OF OBSERVATIONS IN TERMS OF WIND Motion. 
(i) Range Drift—Drift or tilt of the trail as a complete unit. 


It is shown by the work of Lovell, Banwell and Clegg (1947), and by 
Hey and Stewart (1947) that the reflection from a meteor trail is markedly 
aspect sensitive at 72 Mc./s. In fig. 1, if M is the line of motion of the 


Fig. 1. 


meteor, most of the echo energy received at O will come from the region 
close to the foot A of the perpendicular from O to the trail. The elliptical 
curve represents the coverage of the aerial over a horizontal plane at a 
height of 100 km. The meteor trails which give echoes are contained 
completely within this area *. 

If the trail drifts as a unit nearly parallel to itself to a position M,, 
under the influence of winds, the reflecting region will move to A’. The 
velocity of such movement should be of the order of 50 m./sec., and could 


*The polar diagram has a volume distribution for meteor collection 
(Clegg 1948), but for present illustrative purposes the section of this volume is 
sufficient. 
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be either towards or away from the receiver. The change of range in 
drifts of this order is very slight, being only 1-25 km. in 25sec. Drifts. 
of this magnitude have been found by McKinley and Millman (1949) who 
report that the majority of long-enduring discrete range echoes exhibit a 
drift in slant range at the rate of approximately 1-0 km. in 25 sec. Their 
result is in good numerical agreement with the hypothesis of a general 
wind drift. The further observation by McKinley and Millman that the 
drift is predominantly away from the observer still awaits explanation. 


(ii). Range Drift—Gradual curvature of large sections of the trail. 


In addition to this slow drift, much greater range movements have been 
seen in a-proportion of echoes. The available experimental evidence on 
these gross range drifts observed on 72-73 Mc./s. is given in Table I. 


TaBe I. 
Observations on Range Drift. 


Equipment Range drifts 


Freq. ‘Aetint tt -"Beana Total No. | No.of | Order and | Velocity 
Author 


(Me. /s.) Bea ciel bs of echoes | travelling sign of of 
ue eon examined | echoes | movement | movement | 


Prentice, 72 Vertical | Main beam 1836 55:0 5 to 10 km, Widely 
Lovell, Yagi Pale (3%) Mostly spread from 
and Appreciable decreasing | afew 
Banwell side lobes range km./sec. to 


1947 60 km./sec. | 


Hey, and Vertical 15° 20% Mostly 0 to 25 
Stewart Yagi decreasing km. /sec. 
1947 range 


Ellyett * y Array of . 6-5° > 1000 (0-1%) | + 7:0 km. + 0:39 
5 Yagis : km/sec. 
at 15° 
elevation. 


* Unpublished. From an analysis of more than 1000 cine records obtained during the 
major showers of 1946-48 using the technique described by Lovell, Banwell and Clegg (1947). _ 


The first point of interest emerging from the table is the agreement 
between the three independent sets of observations with respect: to the 
magnitude and the velocity of apparent movements, which are of a very 
high order ; speeds of 25 km./sec. being by no means exceptional. The 
suggestion was made by Hey and Stewart (1947) that since these drift 
velocities were far in excess of any possible direct air movement, the 
distortion of the trail must have produced a change in the point of reflec- 
tion, leading to the very high apparent rate of movement. Type (ii) of 
the classes of air motion suggested is identical with this mechanism. 
Referring to fig. 2 the initial reflecting region of the trail at A is shifted 
to A’ by progressive distortion. The actual displacement of the trail is 
exaggerated in the figure, the region A’ being visualized as lying very close 
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to its original position along M. The one definite range drift obtained 
from an analysis of the cine records (Table I.) is the only echo available 
with sufficient information to allow an estimate to be made of the wind 
gradient. As the range R=167km. and the drift HA’=-+7 km. in 
18-0 sec., the angular displacement EOA is 16°. This yields a wind 
velocity gradient of 9 m./sec./km. This value, being less than one-fifth 
of the mean wind velocity, appears not unreasonable. 

On the mechanism of fig. 2, since reflection effectively comes from some 
new point along the line M, the range always increases, and must do so in 
all cases where the normal to a linear trail lies within the beam area for 
initial reflection. 


Fig. 2. 


Although a percentage of outward drifts do occur, the major group 
appears to be directed towards the observer. To account for this group a 
slightly different mechanism must be postulated. Consider a meteor 
track, given by the heavy section AA in fig. 3, passing obliquely through 
the aerial beam, so that the projected perpendicular reflecting point lies 
outside the beam angle. In such a case no echo is received on 72 Mc./s. 
unless the track undergoes distortion and presents a suitable directed 
portion A’ within the beam. As the bending increases A’ may decrease in 
range towards A”. ; 

An interesting feature of the results given in Table I. is the marked 
difference between the numbers of echoes showing range drift with 
vertically and horizontally directed beams. This can also be explained 
on the concept of trail distortion. With a vertically directed beam an 
echo can only be obtained by specular reflection if the meteor trail is 
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almost horizontal, or if the meteor emanates from a radiant close to the 
horizon. Under these circumstances the rate will be low (Prentice 1948). 
Hence a larger proportion of the echoes observed with a vertical beam are 
due to trail distortion. A simple arithmetical estimate, taking into account 
the relative collecting areas (Clegg 1948) shows that on this hypothesis, 
the percentage of range drifts observed with a narrow vertical beam should 
exceed that observed with a horizontal beam by a factor of the order of 20. 
This agrees well with the ratio obtained experimentally. 

As the radio wavelength is lowered the echo duration increases, and 
hence the occurrence of reflection at angles other than the normal to the 
initial trail should increase. This has been substantiated by the work of 
McKinley and Millman at 32-7 Mc./s., who find that trails occurring 


Fig. 3. 


visually within + 10° of a direction perpendicular to the line of sight from 
the station to the trail produce radio echoes after a delay of about 1-0 sec. 
whereas for angles of departure lying between 50° and 90° the delay 
increases to 10-0 sec. 


(iii) Turbulent Motion—Relative movement of short sections of the 
trail. 


The third process considered in § 2 concerned the relative movements of 
short sections of the trail. The explanation of the high speed amplitude . 
fluctuations in terms of this type of trail movement has been discussed by 
Greenhow (1950). In the presence of such turbulence the aspect sensitivity 
of the trail for radio reflections should slowly vanish, and appreciable 
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range spreading of long duration echoes should gradually occur, being more 
-marked on the lower radio frequencies. Evidence for this effect has been 
obtained on 72 Me./s. (Prentice, Lovell and Banwell 1947), while McKinley 
and Millman 1949, have found marked evidence of the effect on 32-7 Mc./s. 


§ 4. CoRRELATION OF METEOR TRAIL FLUCTUATIONS WITH DIRECT 
OBSERVATIONS OF IONOSPHERIC FADING. 


Studies of the fading of radio waves returned from the ionosphere 
suggest the presence of moving ionospheric irregularities (Ratcliffe 1948). 
These irregularities are considered to be caused by wind motion (Pawsey 
1935, Findlay 1949, Mitra 1949). These measurements have been made 
at frequencies ranging from 16 Ke./s. to 6 Mc./s., and at such frequencies 
meteor ionization from individual trails will remain an effective echo 
source for periods of many minutes. 


TABLE LI. 


Period of fading due to Relative Movement of Tonized Air at Heights 
above 80 km. 


Equivalent period 
Method Author Period of fading at 72 Mc./s. 
Short Long 


Tonospheric echoes 


Amplitude fading Ratcliffe 5/min. at 3 Mc./s. 0:7 sec. 
curves (1949) (approximate mean 
from spread of results) 
Phase fading Findlay 1/sec. at 2-4 Me./s. 0-03 sec. 
(1949) 
Phase and amplitude Bracewell 0-2/min. at 16 Ke./s. | 0-07 sec. 
fading (1948) 0-67/hr. at 16 Ke./s. 1:2 sec. 


Meteor trail echoes 


Amplitude flutter Greenhow 0-04 sec. at 72 Mc./s. 0-04 sec. 
(1950) 0-09 sec. at 36 Me./s. 0:04 sec. 
Longer period fading Greenhow 1:0 sec. at 36 Mc./s. 0-5 sec. 
(1950) 
Amplitude flutter McKinley and| 5-10/sec. at 32:7 Me./s. | 0-09— 
Millman 0-05 sec. 


(1949) 


It has already been suggested by several authors that sporadic E effects 
are in part due to meteor ionization (summary by Lovell 1948). Evidence 
will now be produced which indicates that some of the irregularities causing 
the fading of radio waves from the ionosphere may also be due directly to 
wind distortion of meteor trails. 

Since no measured wind velocities are yet available from meteor studies, 

‘the mean period of fading must be used as the point of comparison. In 
order to allow comparison of the results of various workers, it will be 
assumed that the speed of fading is proportional to the radio frequency 
employed. This assumption has already been made by Ratcliffe (1948), 
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who has found some evidence of its validity. Further experimental 
evidence supporting this assumption has now been obtained by obser- 
vation of the fluctuations occurring in the amplitude of radio waves 
reflected from meteor trails on 36 and 72 Mc./s. sumultaneously (Greenhow 
1950). 

Two distinct periods of ionospheric fading have been recognized on the 
very long wavelengths (Bracewell 1948) and although a large scaling up in 
frequency is required it is remarkable that the figures at 72 Mc./s. agree 
closely both with the short period flutter and with the longer period 
observed from meteor trails. A recent figure quoted by Ratcliffe at 
intermediate frequencies agrees with the longer period, and the figure 
given by Findlay (1949) is almost identical with the meteor flutter period. 

There appears, then, to be a close correspondence in the flutter observed 
from meteor trails and the ionospheric fading obtained from unspecified 
“centres ’’. Both show similar characteristics and have the same mean 
period. It therefore seems not unlikely that some of the ionospheric 
fading is, in fact, directly due to meteor ionization at E region heights.. 

In the absence of further results, no great significance can as yet be 
attached to the correspondence of the longer periods. 
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By A. G. Caruson *, J. E. Hooper and D, T. Kine, 
The H. H. Wills Physical Laboratory, University of Bristol f. 


[Received May 24, 1950.] 
[Plate XXIV.] 


SUMMARY. 

The spectrum of the y-radiation in the atmosphere at 70,000 ft. has 
been determined by observations on the scattering of pairs of fast electrons 
recorded in photographic emulsions exposed in high-flying balloons. The 
detailed form of the spectrum is consistent with the assumption that the 
y-rays originate by the decay of neutral mesons. It is found that the mass 
of the neutral mesons is 295-20 m,, and that they are created in nuclear 
explosions with an “energy spectrum” similar to that of the charged 
7-particles. 

A method is described for determining the lifetime of the neutral 
mesons and their frequency of occurrence compared with charged 
a-particles. The lifetime, 7, is less than 5107 sec. It may be 
possible to determine the lifetime, by observations of greater statistical 
weight, if it is longer than 2x10-“sec. The ratio of the number of 
neutral mesons to charged mesons, produced in nuclear explosions of 
great energy, is equal to 0-45+0-10. 


1, INTRODUCTION, 

ReEcENT experiments by Bjorklund, Crandall, Moyer and York (1950), 
strongly suggest that the high energy y-rays observed when matter is 
bombarded by protons with an energy greater than about 200 MeV. are 
of secondary origin; that they arise by the decay into two quanta of 
neutral mesons, of mass about 300 m, and lifetime less than 107" secs., 
which are the primary products of the nuclear collisions. Further, 
Bradt, Kaplon and Peters (1950), have made a detailed analysis of a 
disintegration produced in a photographic emulsion by an «-particle with 
an energy between 10” and 10¥ eV., which leads to the emission of fifty- 
seven shower particlés. In the “core” of the shower particles, most of 
which are mesons (Piccioni 1950, Fowler 1950), these authors observe the 
production of several pairs of electrons, and conclude that the charged 
mesons are accompanied by y-radiation. If this radiation is also created 
by the decay of neutral mesons—and if the latter particles are produced 
in the nuclear explosion with values of the kinetic energy similar to those 
of the charged mesons—it may be shown that the proper lifetime of the 
neutral mesons is less than 3 x 10~1 sees. 


* On leave from the University of Upsala. F 
+ Communicated by Professor C. F. Powell, F.R.S. 
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Finally, the existence of neutral mesons of mass about 280 m,, and their 
mode of decay into two y-rays, appears to have been put beyond doubt as 
a result of experiments by Panofsky, Aamodt and York (1950) on the 
capture of 7--particles in liquid hydrogen. This process leads to trans- 
mutations which may be represented by the following equations :— 


(i) Hi+a-1 > n°+7° ; a —> 2hy (hv~70 MeV.) ; 
and 
(ii) Hi+a-1+ n°+hv (hv~140 MeV.) ; 


7® representing a neutral meson. 

The observed degree of homogeneity of the y-rays from reaction (i), 
which is the more probable, indicates that the 7°-particles are emitted 
with small kinetic energy. This observation, and the known masses of 
the other particles involved in the transmutations, proves that the mass 
of the particle 7° is only a few electron masses less than that of the charged 
a-particles ; viz. m,.~280 m,. 

These observations appear to have an important bearing on experiments 
carried out in the last few years with Wilson chambers by Cocconi e¢ al. 
(1946), Fretter (1949), Chao (1949), Cocconi (1949), and Gregory et al. 
(1950). These authors have shown conclusively that the cascade showers 
of the soft component are commonly associated with penetrating showers 
of charged mesons. 


2. SCOPE OF THE PRESENT EXPERIMENTS. 


In Section I. of the present paper, observations are described which show 
that the detailed features of the spectrum of the y-radiation accompanying 
the penetrating showers are consistent with the assumption that the quanta 
are created by the decay of neutral mesons of mass 295+-20 m,; and that 
the distribution in energy of the neutral particles is closely similar to that 
of the charged mesons of the “ showers ”’. 

In Section IT., a method of determining the lifetime, 7, of the neutral 
mesons is described, and the results prove that 7,.<510-™ sec. There 
is some indication that 7,.—3x10- sec. ; and if this conclusion can be 
confirmed by further measurements, the observations would establish the 
independent existence of neutral mesons. Finally, the results show that 
the relative frequency of production of the neutral and charged mesons, 
AN (n°)|N (7*)=0-45-+-0-1. It appears, therefore, that on the average, 
one neutral meson is produced for every two charged z-particles in the 
explosive disintegrations which give rise to the penetrating showers. This 
result, taken in conjunction with observations by Rossi on the relative 
amounts of energy carried by the “ hard”’ and “ soft ’? components of 
the cosmic radiation in the atmosphere, shows that the latter arises largely 
as a result of nuclear interactions, but indirectly through the production of 
neutral mesons; and that at 70,000 ft., the fraction which need be 
attributed to other processes such as the decay of u-mesons, ‘‘ knock-on ”’ 
electrons, etc. is negligible.. 
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Section I. 
Determinations of the Mass and the Distribution in Energy of Neutral Mesons. 
3. THE y-Ray Spectrum at 70,000 FT. 


During the past year observations have been made in this laboratory on 
photographic plates exposed at an altitude of 70,000 ft. ; see Camerini et al. 
(1949). A number of these plates have been examined under the micro- 
scope with relatively high magnification, in order to detect pairs of electrons 
and other events involving only particles of charge e moving at relativistic 
velocities. Such events commonly escape observation when the plates are 
searched at lower magnification. By determining the energy of the pairs 
of electrons by the scattering method, in those favourable cases in which 
the tracks in the emulsion are of sufficient length, it has been possible to 
deduce the energy spectrum and direction of motion of the photons of the 
“ soft ’’? component at 70,000 ft. 


TaB_eE I. 
Energy region Cell length Probable error in energy 
(MeV.) (microns) (per cent) 
10-100 100 19 
100-200 200 28 
200-400 250 32 


The above errors are given on the assumption that the track under 
consideration is just 1000 microns long. 


The method employed in the determination of the energy of the electrons 
was similar to that described by Fowler (1950) and King (1950). Only 
those pairs were accepted for measurement of which the track of each 
component was longer than 1mm. In the case of tracks of length less 
than 2000 p, it is not possible to determine the energy of the corresponding 
particles unless it is less than 500 MeV. On the other hand, satisfactory 
estimates can be made up to 400 MeV. with tracks only 10004 long. 
Characteristic values of the probable errors of the measurements, for tracks 
1000 1 long, are shown in Table I. The energy of a pair is usually 
substantially greater than that of one of its components. We may, 
therefore, say that the total energy, E, of a pair may be measured with a 
probable error of less than 30 per cent, for values of E <600 MeV., if both 
tracks are longer than 1000; and for E<1500 MeV., if the tracks are 
longer than 2000 py. 

The results of measurements on pairs of very low energy require a 
small correction to allow for the influence of scattering on the escape of 
particles from the emulsion. Such corrections are, however, very small 
for E> 20 MeV., and they have not been applied in making the present 
measurements. The results are shown in fig. 1. 


302 
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In determining the orientation of the pairs of electrons with respect to 
the vertical line through the plates during the exposure, it was considered 
sufficient to determine only the inclination, 0, of the projection of the 
‘bisector ” of the pair. For pairs of electrons produced by radiation with 


Bigeel; 


NUMBER of PAIRS per 50 MeV. 


100 500 
MEV {000 
Distribution in energy of the pairs of electrons observed at 70,000 ft. The full 
line shows the expected photon spectrum arising by the decay of neutral 
mesons of mass 280 m, and with kinetic energies similar to those of charged 
a-particles emitted from stars. 


Angular distribution of the directions of motion of photons at 70,000 ft. 
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values of the quantum energy in the range under consideration, the 
bisector ” gives a measure of the direction of the parent y-ray which is 
subject to a probable error of less than 0-2°. 

The distribution in the values of 6, deduced from observations of 500 
pairs is shown in fig. 2. It will be seen that the directions of motion of the 
quanta are distributed nearly isotropically in the range of angles from 
0° to 90° with the direction vertically downwards, and that the intensity of 
any upward directed radiation is relatively very small. Charged mesons 
ejected in nuclear disintegrations in this exposure exhibit the same angular 
distribution ; Davies, Franzinetti, and Perkins (1950). 


4. CHARACTERISTICS OF THE RADIATION PRODUCED BY THE DECAY OF 
NEvutTRAL MEsons. 


Following the publication of the results of Bjorklund et al., it occurred 
to us that if the y-radiation in the atmosphere arises by the decay of 
neutral mesons—and if it is examined in conditions in which it has not been 
seriously modified by “ bremsstrahlung ”’ and other processes associated 


hy, 


with the development of the cascade showers of the soft component— its 
spectrum should display certain features of an internal consistency, and 
should allow a determination to be made of the mass of the postulated 
neutral meson. The method depends on the following considerations :— 

Suppose a neutral particle of mass m,), moving with velocity fc, trans- 
forms into two quanta, of energy hv, and hv, emitted in directions making 
angles 6, and 0, with the line of motion of the parent particle—see fig. 3. 
The equations for the conservation of energy and of momemtum may be 


written : 
h h 
hyy sin 0,= a sin Oo $ — cos 6,+ = COS O.=MPe/4/(1—B?) 5 
hy +-hvg=mye*/ (1 —f?)=Bmye*=e. 
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It is then easily shown that the energy of the quantum emitted at an angle 
6 is given by the well known formula for relativistic Doppler effect, 


~~ mc? 
~ 2B(1—f cos 6)’ 


Characteristic curves showing the variation of hy with 6, for different 
values of B and taking m,c? equal to 140 MeV., are given in fig. 4. It 
follows from (1) that for a given value of 8, the maximum and minimum 
values of hy are given by the relations 


hv (1) 


Bm,¢* 
2 


hy (max)= epee (1-+8); hy (min) = 


(1-6 )r 


Polar diagram showing the variation of energy with the direction of emission of 
the photons resulting from the decay of neutral mesons. Calculations for 
two values of the total energy of the neutral mesons. 


and that the greater the velocity of a homogeneous beam of neutral mesons, 
the greater the width of the spectrum of the emitted y-radiation. Assum- 
ing that the y-radiation is emitted isotropically in a coordinate system 
moving at a velocity equal to that of the mesons, the form of the spectrum 


of the y-rays can be calculated—see Appendix—and can be shown to be 
given by the radiation relation 


N(E) . dE=k . dE, 


where N(E).dE is the number of emitted quanta of energy between 
EK and E+dE. & is a constant for values of E lying between hv (max) 
and hy (min); and k=0 for values of E outside this interval. Character- 
istic spectra for different values of B, normalized to correspond to the same 
total number of quanta, are shown in fig. 5. 

In the case of an inhomogeneous beam of neutral mesons, dispersed in 
velocity, the resulting y-ray spectrum can be shown to have the features 
illustrated in fig. 6. There is a cusp in the spectral distribution if a finite 
number of neutral mesons are emitted with zero velocity. For any value 
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Fig. 5. 


ARBITRARY UNITS 


500 {000 MeV ISCO 


ENERGY DISTRIBUTION of PHOTONS from a 
MONOKINETIC BEAM of NEUTRAL MESONS. 


Distribution in energy of the photons produced by the decay of neutral mesons. 
Hig, 6. 
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DS 
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Number of photons per 20 MeV. 


) 100 200 300 400 500 600 700 800 
MeV. 


Spectrum of the y-radiation at 70,000 ft. The full line is that giving the best 
fit with the experimental results, and the mass of the neutral particle has 
been deduced from it. 
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of the intensity below the maximum, there are two corresponding values, 
of the quantum energy E, and E,, and these can be shown to be related to 
the rest-mass m, of the parent neutral particles by the relation 

4/ (Hi Be)==7g02- 9. - 
The proof of this theorem is included in the Appendix. 

It follows from the above result that if most of the y-radiation at 
70,000 ft. is produced by the decay of neutral mesons—and if it is not 
seriously modified by the processes which lead to the formation of the 
cascade showers—there should be an inner consistency in the form of the 
spectrum corresponding to the relation displayed in equation (2). It 
appeared to us that the observation of such a consistency would give very 
powerful support for the view that the y-radiation does in fact arise by a 
decay of neutral mesons and would allow an estimate to be made of their 
mass. 


TABrE Ui: 


Measurements of the Neutral Meson Mass from the Observed 
Energy Distribution of the Photons. 


Weights 
(No. BE uaa) €)/2=V(E, Es) 


Weighted 


75-5 
75:0 
75:0 
43:0 
79-0 
75:5 


Total 96 Weighted mean for ¢)/2—75-8 MeV. 
All energies are given in MeV. Mass of neutral particle mj, —=295+20m,. 


5, COMPARISON WITH EXPERIMENT ; THE MASS OF THE NEUTRAL MESON. 


Fig. 6 shows the distribution in energy of the y-radiation. It has been 
deduced from that of the pairs (see fig. 1) by assuming a variation with 
energy of the cross-section for pair production given by the quantum 
electro-dynamics (Heitler 1944). The full line in best accord with the 
distribution has been drawn, and intercepts made at various arbitrarily 
chosen intensities to determine values of E, and E,. The corresponding 
values of 4/(H,E,) are shown by circles. It will be seen from fig. 6 and 
Table II. that they show a remarkable degree of consistency, and indicate 
that the y-radiation arises by the decay of neutral mesons of mass 
295+20m,. The results may be subject to a small systematic error 
arising from a variation with energy of the scattering constant k, defined 
by the equation E=k/«, where « is the mean, observed angle of scattering 
for a cell length of 100 u. Such a variation may exist as a result of the 
effects of screening. 
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The stated errors in the estimated mass of the neutral meson corre- 
spond only to those due to statistical: fluctuations in the numbers of 
measured pairs of electrons, and when the result was first obtained it was 
not clear to what extent it was influenced by the development of the 
cascade showers. The recent observations of Panofsky et al. (1950) 
suggest, however, that the “ distortion ”’ in the form of the spectrum is 
small, and this appears to be associated with the following features of the 
experiments. 

Firstly, the plates were exposed near the main meson production layer 
in the atmosphere. At this altitude the probability that a given quantum 
will have produced a cascade, even in the case of quanta produced at the 


Fig. 7. 


NUMBER of MESONS 


100 MeV 1000 
ENERGY 


The distribution in energy of the neutral mesons emitted from nuclear disinte- 
grations. Spectrum of neutral mesons shown thus: x ; and of charged 
mesons, as observed by Camerini et al. (1950), shown thus: O. 


top of the atmosphere, is small. Secondly, if the cascade process starts in 
the plate assembly itself, it will commonly give rise to relatively narrow 
showers, made up of a number of particles moving nearly parallel to one 
another and close together. The development of such cascade showers 
has been observed in the present experiments, but they are rare, and any 
pairs associated with them have been excluded from the analysis. 
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For these reasons we regard our observations as giving very strong 
evidence for the view that a large part of the y-radiation in the atmosphere 
is due to the production and decay of neutral mesons, of mass ~295 +20m,, 
which may be identified with the particles called for by the results of 
experiments by Bjorklund et al. (1949), and of Panofsky ¢é al. (1950). 


6. DisTRIBUTION IN ENERGY OF THE NEUTRAL MESONS. 

From the observed form of the y-ray spectrum it is possible to deduce 
the distribution in energy of the parent neutral mesons and the result, 
determined by the methods described in the Appendix is shown in fig. 7. 
The corresponding distribution for the charged mesons, emerging from 
the nuclear explosions recorded in the same plates and determined by 
Camerini et al. (1950 a), is included in fig. 7, and it will be seen that the two 
‘spectra’ are similar in form. Alternatively, the spectrum of the y-rays 
which would be produced by the decay of neutral mesons, of mass 280 m, 
and with an energy distribution similar to that of the charged mesons, can 
be easily computed. Using the results of Camerini et al., for the charged 
mesons, the result thus obtained is shown by the smooth line in fig. 1. 
It will be seen that the calculated form of the spectrum is indistinguishable 
from that observed, and added weight is thus given to the basic assumption 
regarding the origin of the y-rays, and for the view that the parent neutral 
mesons are created with a distribution in energy closely similar to that of 
the charged particles. 


Section II. 


Observations to Establish the Independent Existence of Neutral Mesons, and 
to Measure their Lifetime and Frequency of Occurrence. 
7. A METHOD OF DETERMINING THE LIFETIME. 

The evidence for the existence of neutral mesons presented hitherto, 
whilst very strong, is circumstantial and depends on observing y-radiation 
of which the detailed characteristics are consistent with an assumed origin 
by the decay of neutral mesons. Although the mass of evidence which has 
now been accumulated is so large that the conclusion to which it leads can 
hardly be doubted, it is nevertheless important to obtain an unambiguous 
proof of the independent existence of the particles. Further, it has been 
pointed out by Fukuda and Miyamoto (1949) that the determination of the 
lifetime of the neutral meson is of great importance for distinguishing 
between different theoretical treatments of the nature of the particle. 
The following considerations indicate that such evidence could be obtained 
by an application of the photographic method if the lifetime of the particles 
is greater than about 2x 10-' sec. 

Suppose that a neutral meson emerges, with a total energy «= Bm,c?, 
from a nuclear explosion occurring in the emulsion; and that it decays 
into two photons, after traversing a distance 1. We have seen that the 
two y-rays will in general be emitted in directions inclined to that of the 
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parent particle. Suppose, further, that one of the photons produces a 
pair of electrons in the emulsion within a distance of 1 or 2 mm. of its point 
of creation. The probability of such an occurrence is about 3 per cent for 
each emitted neutral meson ; for the average length of path in the emulsion 
of a y-ray before pair production is 46mm. *, there are two photons 
produced by each neutral particle, and a considerable proportion of the 
events will escape observation because of the finite thickness of the 
emulsion. 

We have seen that for neutral mesons with kinetic energies in 
the range of values with which we are concerned, the pair of electrons will 
move in directions inclined to one another at an angle of less than 0-2°. 


Principle of the method for the determination of the lifetime of the neutral 
meson. 


parted to the nucleus involved in the form- 


the bisector of the angle between the tracks 
re of the direction of motion 


f the order of 0-2° (see fig. 8). 


Apart from the momentum im 
ation of the pair (Bethe 1934), 
of the two electrons therefore gives a measu 
of the parent y-ray which is subject to errors 0 ig. 8 
Tt follows that if the mean line of the pair is produced backwards, it will 
not pass precisely through the disintegrating nucleus in which the neutral 
meson originated, but at a distance ” fromit. The greater the value of the 
lifetime, the greater the average distance ] and the corresponding values 


ofr. 


* This value has been calculated by taking 
spectrum of the y-radiation, and the variation of th 


quantum energy. 


king account of the observed energy 
e conversion length with 
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The value of r for any particular event will depend upon the distance 
travelled by the neutral meson before it decays, and the direction of emis- 
sion of the photon relative to the line of motion of the parent particle. For 
a homogeneous beam of neutral mesons, moving with a given energy, 
the expected distribution in the values of r can be calculated by the methods 
outlined in the Appendix. From the results thus obtained, it is then 
possible to proceed to the general case in which the neutral mesons have a 
distribution in energy. Fig. 9 represents the results of such computations, 
assuming an energy distribution similar to that of the charged mesons 
deduced by Camerini et al. (1950 a), and a lifetime of the neutral mesons 
equal to 10-} sec. 


Giga. 


LIFETIME 10° sec 


SYIVd JO YSEWNAN 


100 200 p 
DISTANCE Fr 


The calculated distribution in the values of the quantity r+ for an assumed 
lifetime of the neutral meson of 10-1 sec., together with the calculated form 
of the distribution due to unrelated pairs. 


It is a particular feature of the method that in the case of mesons of 
high energy, B=2, the calculated distributions are almost independent of 
the value of B, and of the distribution in energy of the neutral particles. 
This is due to the fact that although the emitted radiation tends to be 
closely collimated about the line of motion of the parent mesons as B 
increases—a factor which leads to a reduction in the values of y—the effect 
is compensated by the relativistic extension of the time-scale of the moving 
particles. This tends to make the neutral mesons survive over a longer 
length of path in the emulsion. 


8. EXPERIMENTAL RESULTS. 


In order to apply the above method an extensive search has been made 
at high magnification of the region of the emulsion in the immediate 
neighbourhood of nuclear explosions occurring in the plates in which three 
or more “shower” particles are created (n,$3). All pairs of electrons 
observed in the volume of the emulsion defined in the manner illustrated in 
fig. 10 were recorded, and the corresponding values of r determined by 


observations on the direction of the projection of the mean line of the pair 
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in relation to the star, and observations of “dip”. The methods are 


discussed in more detail in a later paragraph with special reference to the 
errors of measurement. 


Fig. 10. 
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The area around a “shower ”’, which is scrutinized for “ pairs ”’. 

The observed distribution in the measured values of r is given in fig. 11, 
and displays a group of events corresponding to small values of r together 
with a pronounced general background. These features may be 
interpreted in the following manner:—We have seen that there is a 
random distribution of pairs of electrons in the plates due to the flux 
of y-radiation, and these will sometimes, by chance, occur in the volume 
of emulsion round a star which is subject to scrutiny. The frequency 
of occurrence of these chance juxta-positions, which we may refer to as 
“unrelated ”’ pairs, can be calculated from the observed numbers, per 
unit area of the plates, of the showers (7,3), and of the pairs of electrons. 


Fig. 11. 
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Observed distribution in the values of r for r<400 yu. The group of observations 
for small values of r, due to “ related ”’ pairs is clearly shown. 


The calculated variation with r of the frequency of occurrence of these 
unrelated pairs is shown by the full line in fig. 11; and it will be seen to 
be in good accord with the observations. Thus in the interval in r from 
40 p to 200m, the expected number of pairs is 22+5 and the observed 
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number 23. On the other hand, most of the fifteen pairs for which r<40 pp 
must be attributed to y-radiation associated with the neighbouring 
explosion. Thus the expected number of unrelated pairs for this interval 
in the value of r is only 0-2 compared with the 15 observed. The 
observations therefore prove that y-radiation is produced, by some 
mechanism, as a result of nuclear interactions leading to the production 
of shower particles of multiplicity n,>3. 

The observations on the related pairs prove that the lifetime of the 
neutral mesons is less than 5 «10-14 sec., and they appear to indicate that 
t,=3 10-4 sec. It has already been emphasized that such a result, 
apart from its interest in establishing the decay-constant for the neutral 
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Determination of the quantity r. 


m-particles, would also provide a direct proof for their independent 
existence. It is therefore important to examine what weight can be 
put on the observations; and, in particular, whether they can be 
interpreted, alternatively, as due to errors of measurement, the lifetime 
of the 7°-particles being much shorter than can be established by this 
method. 


9. ERRORS OF MEASUREMENT IN THE DETERMINATION OF THE QUANTITY ¢. 


The quantity r which is taken as a measure of the relationship of a 
pair to a particular nuclear explosion, is illustrated in fig. 12. It is the 
perpendicular distance from the centre of the star to the bisector of the 
lines of motion of the two electrons into which the y-ray transforms. 
To determine r, the perpendicular distance, x, from the star to the vertical 
plane containing the bisector of the electron-pair is first measured. By 
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“ vertical ”’ is meant here the direction normal to the plane of the emulsion 
which is taken as “ horizontal’’. «is represented by the line OC in fig. 12, 
and is determined by turning the plate until the bisector is parallel to 
one of the stage motions, Y, and coincident in the field of view with a 
line in the eye-piece graticule. The graticule is also provided with a 
scale perpendicular to the line. After moving the stage parallel to the 
Y axis, and refocusing, the star is eventually brought into the field of 
view. Its distance from the line can then be measured by means of the 
eye-piece scale. 

The above method of determining x is subject to errors which depend 
on the distance BC=d in fig. 12, and on the precision with which the 
bisector of the pair can be defined. Studies of scattering of fast particles 
(Davies et al. 1949) indicate that the alignment with a rectilinear track 
can be made with an error of less than 0-1°. In the case of the pairs 
produced by high energy photons, the two tracks are not appreciably 
scattered in the first 100 u of the trajectory, and the angle between them, 
5, is small; for 100 MeV. y-rays, 5~0-2°, and for 400MeV., 5~0-05°. 
For such energetic radiation, which constitutes most of the spectrum 
(see fig. 1) the definition of the bisector can be made with a precision of 
the order of 0-25° and the corresponding error in x can be written dv~ 500° 
dx and d both being measured in microns. 

The second quantity necessary for the determination of r is obtained 
from observations on the angle of “dip” of the bisector of the pair. 
In fig. 12, AC is the projection of the bisector on a plane, parallel to the 
surface of the emulsion, which contains the centre of the star, O, and f 
is the depth below this plane of the point of origin of the pair, Ke At 
a distance Ad along the bisector, at the point E, the depth is f+Af. It 
is then easily shown that 


2—CD?=(d. Af—f . Ad)?/(4d?+ 4f?), 
and that 
it a a 


In order to determine the probable errors resulting from the errors 
of measurement of the depth, numerous individual observations have 
been made on a number of events, and the standard deviations in the 
results deduced empirically. Similar observations have been made on 
protons of great energy which are observed to produce disintegrations, 
and which are known from scattering measurements to suffer an average 
change in direction along the trajectory of less than 0-01° per 100 p. 
Values of “7” can thus be deduced for tracks which are known to pass 
through the centre of the star, and of which the departure from 
rectilinearity is negligible. 

These observations indicate that the probable errors in the determination 
of the angle of dip are about 0-2° and that the uncertainty in the true 
values of 7 is contributed to almost equally by the two measurements, 
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ofwandofz. According to the particular characteristics of the individual 
events, one or other of the two measurements may be the more important. 
Taking both measurements into account, the probable error in the 
determination of the direction of the bisector is estimated to be 0-3". 
The important question remains of the precision with which the 
bisector gives a measure of the line of motion of the y-ray. In general 
the lines of motion of the two electrons are not equally inclined to that 
of the y-ray which produces them. Further, the nucleus involved in 
the creation of the pair receives an impulse of order of magnitude m,c. 
(Bethe 1934, Modesitt and Koch 1950). At the energies with which we 


Taste III. 


Individual ry measurements and their errors. 


Pair | Star type x dip= 4f/Ad 


0-050 -+0-013 2-4+0-5 
0-096 +0-010 9-5+0-7 
0-120+0-010 | 60-0-+8-0 


0-122+0-004 | 91-5+0-8 
1-011+0-020 | 132-9-+-0-4 
0-094-+40-004 | 31-8-+0-5 
0-172+0-003 | 113-6-+0-4 


MAI Ou Who 


0-158+0-002 | 157-5-+-0-3 
0-337+0-013 | 78-6--0-5 
0-042+0-016 | 27-3-+40-4 


0008+ -0008 | 17-4-+1-2 
0-320+0-006 | 145-0+1-1 


0-049-++0-003 | 81-6+0-8 | 29-:0+4-0 


0-170+0-004 | 94:6-+0-7 | 38-0+4-5 


All the errors are expressed in terms of the standard deviation. 
-a, d, f and r are all measured in microns. 


are concerned, the probable errors resulting from these effects can hardly 
exceed 0-2°, We thus obtain 0-4° as the final figure for the probable 
error in the determination of the line of motion of the photon. It follows 
that for errors in the measurement of the perpendicular distance from the 
star to this line, we may write 6r2d/150. 

Details of the observations on the fifteen related pairs are given in 
Table III. The limits of error given in the last column are those deduced 
from the standard deviations associated with the inconsistence in a large 
number of measurements on the same event. They do not, therefore, 
include contributions due to a deviation between the direction of the 
y-ray and the bisector of the pair. 
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A histogram showing the values of r on a more extended scale has been 
reproduced in fig. 13. Included in the figure are curves indicating 
the distributions that would be expected for standard deviations in the 
errors of measurement, arising from all sources, of 0-25 and 0°5°, the 
lifetime of the particles being assumed to be short compared with 
10-™ sec. It will be seen that there are only two or three pairs for 
which the values of r are substantially greater than those expected 
from a standard deviation of 0-5°. Reference to Table III. shows, 
however, that the values of 7 are many times greater, for these four 
events, than the standard deviations. It is very improbable that four 
such events should be found among a sample of fifteen related pairs, if 
the y-rays diverge, in effect, from the nuclear explosion; and we have 
seen that the number of expected unrelated pairs for r<40 pis only 0-2. 


Fig. 13. 
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Distribution in the observed values of r for related pairs. The curves (A) and 
(B) show the expected distributions due to probable errors of measurement 
in the determination of the line of motion of the y-ray of 0-25° and (O59 
respectively, the lifetime of the neutral mesons being assumed indefinitely 
short. 


The number of events observed is not great enough, however, to allow 
us to put any weight on the indications of a finite lifetime given by the 
experiments, and we only note that if the apparent features can be 
confirmed by further observations of greater statistical weight, it may be 
possible to establish the magnitude of the lifetime if it is equal to or 
greater than 2x 10-™ sec. The conclusion that the lifetime is less than 
510-14 sec. appears to be established. It will be remarkable if the 
lifetime of the neutral particles is found to be just sufficiently great 
to be determined by this method, the only one available with the resources 
at present at our disposal. 
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10. THe INTENSITY OF THE NEUTRAL MESONS. 


In order to determine the frequency of occurrence of the neutral 
mesons in comparison with the charged mesons, the total length, L, of 
the tracks of the shower particles in the scrutinized area round each star 
has been determined. Recent experiments indicate that about 80 per cent 
of these shower particles are 7-mesons, and the contribution of fast 
protons can thus be excluded. Suppose that, on the average, the emission 
of each charged z-particle is accompanied by ¢ neutral mesons : 


N (°)|N (w*)=¢. 


The length of path in the emulsion of the y-rays formed by the decay of 
the neutral mesons is then given, without serious error, by 


4 
=L.2.$=1-6L¢. 


The observation show that in this length of path, 15 pairs of electrons 
were created. Taking into account the small change with quantum energy 
of the ‘‘ conversion-length ”’, the average length of path of a y-ray in the 
emulsion before it produces a pair was calculated to be 4:6 cm., for the 
particular spectral distribution with which we are concerned (Heitler 


~ 1-6L 
1944). We may therefore write i a =4:6; with L=95 em. it follows 


that 
o=0:45-40-10. 


The above result suggests that in the production of showers of charged 
particles, for every two charged mesons, one neutral particle is created. 
This conclusion receives support from a recent study by Camerini, Fowler, 
Lock and Muirhead (1950 b) in this laboratory, of the energy balance 
in the nuclear processes leading to the production of showers. These 
authors have compared the energy of the primary particles with that of 
the mesons and nucleons ejected in the disintegration. They have thus 
been able to show that there is an apparent disappearance of energy 
which can be accounted for if neutral 7-particles are created. A similar 
result has been obtained by Green and Fretter (1950). 


11. ConcLUsIONs. 


We regard our experiments as proving that in the explosive disinte- 
gration of nuclei, neutral mesons are created in numbers equal to about 
half those of the charged z-particles. The mass of these neutral mesons 
is 295+20m,, and they disintegrate into two y-rays with a lifetime less 
than 5x10~-™ sec.; they may therefore be identified with the neutral 
a-mesons discovered in Berkeley. It is well known from the work of 
Steinberger (1949), and of Fukuda and Miyamoto (1949), that the mode 
of déxcay into two y-rays indicates that the neutral meson is of spin zero. 


\ 
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The distribution in energy of the neutral mesons is closely similar to 
that of the charged mesons, and they arise in nuclear processes of the 
same type. They may therefore be regarded as having properties similar 
to those of neutral “‘ heavy quanta” of Yukawa’s theory. 

The observed frequency of occurrence of the neutral mesons allows a 
simple explanation to be given of an observation by Rossi (1948) 
that the flows of energy represented by the “hard” and “soft” 
components of the cosmic radiation are of closely similar magnitude ; 
an observation which appeared as an empirical fact at the time but 
which, as Rossi suggested, might have an important, deeper significance. 
Our results confirm the suggestion that both the “hard” and “ soft ” 
components arise in the nuclear disintegrations produced by protons and 
heavier particles of great energy, which are accompanied by the creation 
of charged and neutral z-particles. It appears to be unnecessary to 
attribute any of the soft component to “ primary ”’ electrons. 

The lifetime of the 7°-particles is certainly less than 5 10-" see. ; 
and it will be difficult to determine it directly, and to establish the 
independent existence of the particles, with the methods at our disposal. 
A determination of the lifetime may, however, be possible by the present 
methods if it lies in the interval from 2 x 10—" sec. to 5x 10-" sec. 
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For completeness, we outline the main features of the calculations. 
We have seen, that the energy of a quantum emitted at an angle 6 with 
respect to the line of motion of a neutral meson of velocity v=Be, total 
energy Bm,c’=<«.is given by the equation 
i MC” 
~ 2B(1—f cos 6) ° 

5 2 


hv=E (1) 
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(a) Angular Distribution of the Radiation. 


The quanta produced by the spontaneous decay of a homogeneous. 
beam of neutral mesons may be assumed to be emitted isotropically in 
a coordinate system—the C-system—moving with the same velocity as. 
the particles. The motions in this system are then transformed to the 
laboratory system—the L-system. 

Suppose the angle of ejection of a quantum in the C-system, relative to 
the direction of motion of the meson, is ~«. The components of the 
velocity, c, of the quantum, parallel and perpendicular to the direction 
of motion of the particle, are c. cos « and c. sin «, respectively (fig. 14). 
Let uw, and u, be the corresponding components in the L-system. By 
the Einstein theorem for the addition of velocities we obtain 


where B=1/,4/(1—8?). It follows that 


U sin a 
tan @= -2 = 


u,  B(cos a+)” 
Writing sin €=/(1—f) and cos =f, it follows that : 


tan 6— Sie sin e ep taehas ay 
cos a-++cos € ern © ) 
From the identity sec’ @=tan* 6+-1, we then obtain the useful relations 


cos cos 
cos 6= ot os 


igecos sé cosa fc 
and 
: sin € sin « 
ge a 
TEP Goats Cone aa ee es 


and from (2), (3) and (4) 


da (cos €+cos«)?. sec? 


d§ sin &(1fcos&.cosa)’ © * °° «+ @) 
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Let the intensity of the radiation in the L-system, measured as the number 
of quanta per unit solid angle, be I(#); and in the C-system, I(«). We 
can then write 


I(@) . sin @ . dO=I(«) . sin« .da=K . sina. da. 


Putting K=1 for convenience, we obtain 


sina da 
a sin 0° ao’: 
which reduces to 
DOP COS:0)" oo whs: 6 ss eomep ee hO) 


The form of the distribution, for two values of B, is shown as a polar 
diagram in fig. 15. 


Fig. 15. 


Polar diagram showing the distribution in angle of the intensity of the y-rays 
produced by the decay of homogeneous beams of neutral mesons of different 
energies. 


(b) Distribution in energy of the gamma-rays. 
We may write N(E) .dE=31(@).sin@.d@. From (1) and (6) we obtain 
dé _—-.2(1—f cosé)? 
dE ~—S Be(1— 8? )sin 8” 
and, thence, 
IGG ee Pe==CONSUAT te. °... . lee Wea) 
The distribution extends from 
Reeerieea\e2- “to Hj (1 A)e2e 7 (8) 


(c) Properties of the y-ray spectrum. 

We now consider an inhomogeneous beam of neutral mesons distributed 
in energy, the number in the interval from « to «+de being F(e) . de. 
The y-ray intensity due to the neutral mesons in this energy interval is 
given, according to (7) by the equation 


GT AN (9) 
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between the limits 
E,=}{e—+/(e—e)} and E,=${e+ V(e—eG)}, - - (10) 


where ¢, is the rest-energy of the neutral meson. At a y-ray energy, E, 
the intensity is obtained by integrating equation (9) between the limits 
«=E+<2/4E and infinity ; for every neutral meson with energy above « 
contributes to the y-ray intensity at the energy E. This holds both for 
energies E, <e,/2 and E,><«,/2, because, when solving (10) for «, we obtain 
in both cases 


sek pedi. eee US a ee 
Thus the y-ray spectrum is given by 
aoa He) de 
N(E)= | —>—s- -- 12 
@=) Ved “ 
By differentiating (12) we obtain 
F(E+¢?/4E)=E ax (13) 


a relation which allows the distribution in energy of the mesons to be 
derived from the observed spectrum of the y-radiation. At the quantum 
energy E=e,/2, the intensity of the y-radiation has a maximum value. 
It is clear from (13) that the slope, dN/dE, at this point is zero only when 
F(e))=0; 7. ¢., when there are no mesons of zero velocity in the beam. 


(d) The Rest Energy of the Neutral Meson. 


It follows from (12), which holds both for E,<e,/2 and E,>e,/2, that 
the intensity of the y-radiation is the same for E, and E, when according 
to (11) 

Ej +¢)/4E,=E,+ 5/4, 
which gives 
e?=4E,E, 

or 

gr Wire hig)e cee os) Ae Ah yf. Poe) 
This relation, which permits the derivation of the rest energy, ¢, from the 
observed y-ray spectrum, holds generally, being a.result of the Doppler 
effect. It is therefore independent of the form of the “ spectrum ” of 
the neutral mesons. 


(e) Derivation of the r Distributions. 


Consider a meson emerging from a star (fig. 8), travelling a distance | 
before decaying. Assume that it then emits one of its y-rays in the direction 
making an angle 6 with the direction of motion of the neutral meson, 
so that the orthogonal distance from the star to the line of motion of the 
y-ray isr. It foHows that 


v=), Ain 0% 
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and, as we have N(@) . d0=1(6) . sin @ . d0, we may write 


do 
dra 
Now /=cfBt, where t is the time interval between the instants of creation 
and the decay of the neutral meson in the C-system. We assume that 


the meson decays according to the usual laws of radioactive decay, and 
then find that the number decaying between times ¢ and ¢+dt is 


N(r)=N(6) 1(0) . tan 6/1. 


N(¢) dt= eS (—t/79) dt, 


om) 


T, being the proper mean lifetime of the meson. Hence, 


N(0) dl= ol exp (—l/c8B7p) . dl; 


BBro 
and since 
l=r/sin 8, 
N(J) dl= Nokes pall as exp (—1/cBBry sin @) dé ; 
cBBz, tan @ sin 6 " : 
and 


N(r)= N,/c8B7 - [ "1(6) exp (—r/¢BBr» sin 8) d8. 
/ 0 


We have performed this integration graphically for several values of the 
energy, and found that the calculated distributions in the values of 7 
are nearly exponential for values greater than 10. For values BeS2 
the curves are very nearly identical. 

The final distribution is obtained by adding the curves for different 
values of B, each appropriately weighted according to the energy spectrum 
of the neutral mesons. 


(f) Background Pairs. 


The probability of finding an unrelated pair of which the bisector 
passes the star at a distance lying between r and r+dr may be calculated 
from the known angular distribution and frequency of occurrence of such 
pairs in the plates. 

Our observations show that y-rays are isotropically distributed in 
the lower hemisphere (fig. 2). As a result of the geometry of the scanned 
area (fig. 10), two-thirds of the unrelated pairs point towards the star. 

If N(r) dr is the distribution in r formed by these pairs, we may write 


ns dé 
Ni) =4 | (0) .sin 6. =. N(l) dl 
where / is the distance between the star and the origin of the pair, @ is the 


angle between the bisector and the line joining the origin of the pair 
with the star, and N(J) di is the number of pairs of which the origin lies 
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in the spherical shell of radii] and J+dl. Thus N(1) dl=p4ql* dl where p 
is the observed number of pairs per unit volume. The angular distribution 
1(@) is constant. R is the effective radius of the So area. Using 
the relation r=/ . sin 0 we obtain 


/ 


R 
N(r)=constant x | 1. tan 6 d0 


=constant x r/(R?—??). 


For r<600 py, we have N(r)~const .R.7, a relation which is true to within 
10 per cent. The constant is then deduced from the observed density 
of “ background ” pairs and the known value of R. 
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A large shower of mesons (6+ 4p) showing the origin of a related pair very close 
to the star. The pair has been mar ked by an arrow in the microphotograph. 


Observer: Miss E. James. To face page 724 
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LXIV. CORRESPONDENCE 


Gap Measurement as a Method of Analysing Cosmic Ray Stars 
in Photographic Emulsions. 


By P. E. Hopason*. 


[Received May 31, 1950.] 


OwtnG to the recent increase in sensitivity of photographic emulsions 
for nuclear research, tracks of particles of moderate energy are clogged, 
that is, the grains coalesce to form continuous blocks. Since grain 
counting is a well-known method of identifying and estimating the energy 
of such particles, it is desirable to extend the method to such tracks. 

While in regions of moderate clogging some simple convention, such as 
one micron of solid track is equivalent to three grains, is satisfactory, 
this is highly artificial where the clogging is severe. In such regions it is 
natural to measure the total gap length per unit length, which is defined 
as the gap density. As grain counting is easier and more rapid than gap 
measurement, it is desirable to use it to the highest possible grain densities, 
which in practice are found to be in the region of 1-5 grains per micron, 
which corresponds to a gap density of about 0-1. Beyond this region 
gap measurement must be used. 

In order to relate the gap density to the energy loss, gap measurements 
versus residual range were made on the tracks of a number of particles 
of different mass and charge ending in the emulsion of an Iford G5 plate 
exposed on the Jungfraujoch. Figs. 1 and 2 show the discrimination 
obtained by means of graphs of total gap length against residual range. 
Tf the track is longer than 400 microns, gap measurement gives almost 
certain identification, except between protons, deuterons and tritons, 
and between 7- and y-mesons. The two 7t-mesons are those recently 
reported (Harding 1950). 

The energy losses at points along the tracks were calculated from the 
known range-energy curves (Bradner et al. 1950), and hence a graph of 
gap density versus energy loss was constructed. ‘T’hese results are plotted 
on fig. 3, together with the standard deviations of measurements along 
130 micron samples of track. In practice more than 130 microns of track 
are often available, and so the accuracy is higher than indicated by the 
figure. It is found that, above energy losses of 9 KeV./y, the gap density 
and energy loss satisfy the empirical formula 

ele =0-53 (aR) = ) 
dR dR 
Pepe i Se 
* Communicated by Sir George Thomson, F.R.S, 
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_ where dH/dR is the energy loss in KeV./u and dG/dR the gap density. 


The exact form of this expression depends on the type of emulsion and 
development, and so each batch of plates must be calibrated. 
At lower grain densities, where the clogging is small, the gap density 


is given by 


30 20 ISMeV 


HA PARTICLE ENERGY 


Gap density. 


. x : 
PROTONS ssi 
ALPHA- feel ea 


PARTICLES 


y 
Energy loss in KeV. per micron. 


Relation between gap density and energy loss with standard deviations 
of measurements along 130u samples of track. 


where dN/dR is the grain density in grains per micron empirically 
related to dE/dR, and A~0-6u is the mean grain diameter. This curve 
is included in fig. 3 and agrees with the experimental points for high gap 
densities. It should be noted that the gap density is a pure number, and 
so, unlike the grain density, is independent of the size of the eye-piece 


scale divisions. 
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The procedure to be followed in order to identify and measure the 
“energy of those particles from cosmic ray stars in electron sensitive 
emulsions producing clogged tracks depends on whether the particle 
under investigation ends in the emulsion or not. 

If the track ends in the emulsion, a gap measurement is made and the 
particle identified by comparison with figs. 1 and 2. Visual examination 
of the track for scattering and 8-rays (Sérensen 1949) confirms the 
identification. Use of range-energy curves gives the energy of emission 
of the particles from the star. 

If the track does not end in the emulsion, the gap density identifies 
the particle, provided the track is sufficiently long, except perhaps for 
some uncertainty between a slow meson and a fast proton, and between 
a slow proton and a fast alpha-particle. Inspection of scattering in the 
first case and $-rays in the second, resolves this difficulty. Using fig. 3, 
the energy loss can be found from the gap density and hence the energy 
is known from the energy-energy loss curves. 

The above analysis can only be carried out if the track makes a small 
angle with the plane of the emulsion. The correction factor which must 
be applied to the gap densities of dipping tracks may be derived from a 
simple geometrical model. It may be shown that a gap of length g 
appears of length (g+A) cos 6’—A when viewed from a direction making 
an angle 6’ with the normal to the track. Assuming that the distribution 
of gap lengths may be replaced by a mean gap length g,, and that the 
lengths of blocks of grains remains constant during shrinkage, it may 
further be shown that 


dG cosé ; A(1—cos 6’)\ /dG\’ ; cos 6 
dR cos 6’ gi gr) * ~ cos 6’ 


where tan 6/tan 0’=shrinkage factor=2:7, 


AGN! : 
(an) apparent gap density, 


Jm—=Mean apparent gap length. 


For 6<10° this correction is small, while for 6>20°, gap measurement 
is unreliable. ‘The correction factor is large for tracks of low gap density, 
since the mean gap length decreases with gap density. 

This method, although more laborious, is better than simply counting 
the number of gaps because it gives better discrimination between 
particles and also it is less sensitive to the angle of dip of the track. 
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A Note on the Use of Resistance Thermometers for Measurement of Rapidly 
Changing Temperatures. 


By C. B. Daisy, D. H. Fenper and A. J. Woopatt. 
Military College of Science, Physics Branch *. 


* Communicated by the Authors. 
[Received June 2, 1950.] 


In a recent issue of the Philosophical Magazine (1950, 41, 468-477), Leah, 
Rounthwaite and Bradley describe the use of resistance thermometers in 
studying the high temperatures produced in explosions. These investi- 
gators sought to follow a temperature which rises more than 1000° C. in 
20 or 30 milliseconds by using a platinum alloy wire of length | inch and 
diameter 0-0005 inch, and to eliminate radiation errors by a calibration 
experiment designed to disentangle natural convection from radiation. 

In investigations we are making, we have followed the temperature- 
changes of a long, fine platinum wire suddenly heated and then allowed to 
cool in an enclosure under conditions in which conduction through the gas 
is the dominant method of transference of heat. In such conditions the 
cooling curve of the wire is roughly exponential and experiment sub- 
stantiates our theory, which assigns to the cooling process a time-constant 
depending on the thermal conductivity of the gas, on the thermal capacity 
per unit length of the wire and on the geometry of the system, including the 
dimensions of the containing vessel. The plain wire used by these investi- 
gators would display in such conditions a time-constant of the order of 
20 millisec., and the quartz-coated wire would behave appreciably differ- 
ently. The relatively stout leads would have a time-constant of the order 
of half a minute, but this will be considerably reduced if natural convection 
is allowed to come into play. 

If this system of wires and leads is exposed to a sudden change in 
temperature, since convection currents take time to form, the initial heat 
transfer will be solely by radiation and conduction through the gas and 
along wires. Conduction to the leads is not negligible. Indeed a calcu- 
lation shows that in the extreme case when the surface of the wire loses 
heat by conduction alone, the average temperature of the wire only rises 
by 93 per cent of the rise at the centre, followed by a subsequent slow drift 
as the leads warm up. Since the end effects will be very different in the 
main and the calibration experiments, an error will be thereby introduced, 
‘and the effect will moreover be appreciably different for the quartz-coated 
wire. 

Further, the attack on the radiation correction in the explosion experi- 
ments relies on a subsidiary experiment on natural convection. As 
mentioned above, natural convection can have little bearing on the critical 
first stage of the explosion experiments and the method of obtaining a 
radiation correction is therefore invalid. 
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Lastly, but probably of most importance, the true shape of the temper- 
ature-time curve must be badly distorted by the lag in response of the wire. 
The fact that the plain and the quartz-covered wires have different 
time-constants will ensure that the distortions are not the same in the two 
cases and this factor alone might conceivably account for the differences 
observed in the experiments. Certainly the effect of this factor must be 
very carefully considered before any reliance can be placed on deductions 
from the curves presented. 

It was the conviction that the response of resistance thermometers and 
thermocouples to sudden changes warranted more attention which 
prompted us to undertake our investigations some time ago, and we hope 
shortly to be in a position to make a more adequate communication on the 
subject. 


May 22nd, 1950. 


LXV. Notices of New Books and Periodicals received. 


Relations entre les phénomenes solaires et géophysiques. Méthodes de calcul dans 
des problémes de mécanique. Oinétique et mécanisme des réactions d’inflamma- 
tion et de combustion en phase gazeuse. (Centre National de la Recherche 
Scientifique : Paris.) 


THESE are records of three of the international colloquia which have been 
organized in recent years with the help of the French Scientific Research Centre. 

The first, held at Lyons in 1947, dealt with solar and geophysical phenomena. 
A main problem in this field is to establish the complex relations between the 
incident solar radiations and effects observable at the earth’s surface. Some 
evidence is obtainable by balloon and, more recently, by rocket flights, but 
most of the experimental material presented here is of radio observation on 
the ionosphere and optical investigation of such phenomena as the aurora 
and the luminosity of the night sky. On the theoretical side, there is a paper 
by S. K. Mitra on unsolved problems in the field and articles by H. S. W. Massey 
and Marcel Nicolet on the photochemistry of the high atmosphere. 

The second volume is divided into two parts referring to meetings at 
Marseilles and Paris in April 1948. The first part contains an exposition by 
D.N. de G. Allen of Southwell’s relaxation method, and an account by T. Vogel 
of the escalator method of Morris and Head for the determination of normal 
vibrations. The second is devoted to problems of fluid flow, including in 
particular articles by J. M. Burgers on turbulent flow, and by L. Malavard on 
the use of electrical analogues. 

The third volume records the colloquium on the initiation and propagation 
of combustion in gases which was held in Paris at the end of April 1948. 
It is concerned for the most part with the kinetics and detailed molecular 
mechanism of the reactions involved. It would give a wrong impression 
to cite individual articles. Considerable advance has been made in this subject 
during and since the war, and this colloquium contains much interesting 
and valuable information. 

The benefit of colloquia such as these is greatest for those fortunate enough 
to take part, but these printed records will be of interest to many not so 
experienced in such special field as to have contributed thereto. G. W. 


Notices of New Books and Periodicals received. plod 


The Chemical Elements and their Compounds. By N. V. Stpawick. Two vols. 
[Pp. xxxii+1703.] (Oxford: Clarendon Press, 1950.) Price 70s. 


‘Tus work stresses the basic principles underlying the reactions and properties 
of the chemical elements and compounds. A simple and sensible layout 
based on the classification of the periodic table, and the emphasis placed on the 
relations between elements of a group and the groups themselves, provide a 
clear and satisfying picture of chemistry as a coherent whole. Perhaps the 
most striking feature of the discourse is the rational treatment accorded to the 
organic and metal-organic compounds, which are carefully fitted into the 
pattern instead of being relegated to more highly specialized works. Technical 
details are reduced to a minimum, with the result that the fundamental structure 
of chemistry is thrown into relief. Clear exposition combines with an extensive 
bibliography to make these two volumes a necessity for the chemist’s library. 
Dara: 


Non-Linear Vibrations in Mechanical and Electrical Systems. By J. J. STOKER. 
[Pp. xix+273.] (Interscience Publishers, Inc.) Price 40s. 


‘Tue author is concerned with the free and forced oscillations in two ordinary 
non-linear differential equations, and although there is no claim to give a 
survey of the work of the past twenty years, he does take account of much 
recent work and includes an up-to-date bibliography. 

After a discussion of the stability of the first order homogeneous differential 
equation near a singularity, he considers the forced oscillation in Duffing’s 
Equation, in which a small non-linear element occurs in the restoring force. 
Sections are included on jump phenomena, that is, discontinuities in the 
solution as the period of the forcing term is varied and on the appearance of 
subharmonic oscillations of a permanent type. 

The other differential equation is Van der Pol’s, in which the damping term 
is non-linear, so that the damping is negative for small velocities and positive 
for large velocities. Self-excited oscillations are considered first of all and 
then the stability of the periodic solutions arising from a small forcing term. 
A chapter is included on Mathieu’s equation and its application to the stability 
of the solution of Duffing’s equation. 

The presentation up to this point is elementary and within the reach of 
non-mathematicians, but to conclude the book there are six appendices proving 
some results of a more abstruse nature. K. 


Les spectres de rayons X et la structure électronique de la matiere. By YvTETE 
Caucnors. (Paris: Centre National de la Recherche Scientifique.) 


Tyars short book contains an admirable survey of a subject to which the school 
of Mlle. Cauchois in Paris has made important contributions. The first part 
gives a rapid conspectus of the field and discusses especially the influence of 
chemical binding on the position of the discontinuities in X-ray absorption 
spectra, while the remaining and greater part of the book deals with the 
interpretation of the fine structure of these spectra in terms of the permitted 
electronic energy levels, especially for crystalline absorbers. Emission spectra 
are also fully treated, though in this field there is a much smaller mass of 
reliable experimental work. ae 
This subject is of the highest importance as giving substantial information 
con the behaviour of electrons in molecules and condensed matter, and it is an 
especial merit of the book to draw attention to the decidedly unsatisfactory 
state of theoretical work on these questions. Many of the unsolved problems 
raised are of considerable interest. This book should be read by all interested 
in the electronic structure of matter. G. W. 


he 
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Dynamics of real Fluids. By E. G. Ricnarpson. [Pp. 143.] (London : 

Edward Arnold, 1950.) Price 21s. 

TERE has for some time been a need for an up-to-date text-book on the physics. 
of real fluids, avoiding both the mathematical complexity of Lamb’s classical 
treatise on ideal fluids, and the engineering details of, say, Dodge and Thomson. 
One hoped that Dr. Richardson’s book would meet this need, and a glance at 
the chapter headings seemed encouraging. A short introduction on_ classical 
hydrodynamics is followed by a long chapter on “ Fluids of Small Viscosity ” 
which deals largely with boundary layer theory and turbulent wakes. There: 
follow chapters on compressible fluids, fluids with free surfaces, fluids with 
anomalous viscosities, and problems involving thermal effects. 

As a text-book for students the book seems to the present reviewer to have. 
serious faults. As a summary of the present state of knowledge on the topics 
treated, the book is probably unique: but to the non-specialist it is very 
difficult reading. The author’s very familiarity with his subject-matter is. 
probably a contributory cause—his researches figure in almost every chapter. 
The treatment attempted is largely non-mathematical, but equations are 
sufficiently numerous to necessitate an attempt to understand them, and. 
insufficiently detailed to assist the understanding. The exposition seems too 
rather obscure in places, and the number of mistakes and misprints is such that. 
more than twenty were noticed at a single reading. 

In short the book is excellent in conception, but the execution leaves much to: 


be desired. N. TE: 


An Introduction to Heat Transfer. By M. FisHENDEN and O, A. SauNnDERsS. 

[Pp. 205.] (Oxford: Clarendon Press, 1950.) Price 15s. 

Tuts book is intended primarily for engineering students but others will also. 
find it a useful summary of much information of practical value. It attempts 
to explain the fundamental physical principles of the three mechanisms of 
heat transfer, conduction, convection and radiation, and at the same time to 
provide quantitative formule and numerical data in a form suitable for 
practical applications. In most cases the basic formule are quoted without 
proof, but are followed by an adequate comparison with available experimental 
data. Hach chapter contains a number of worked examples to illustrate the 
use of the formule. 

The first third of the book deals with radiation and conduction in two 
chapters of about equal length, whilst the remainder is devoted to problems of 
natural and forced convection. The chapter on thermal conduction contains 
an account of Schmidt’s graphical method of solving problems of unsteady 
heat flow, and emphasizes the importance of numerical and analogue methods 
generally. The treatment of thermal convection is based on formule which 
satisfy the principle of dynamical similarity, and is preceded by an elementary 
discussion of the application of dimensional analysis to convection problems. 

The present reviewer feels that it is a valid criticism that this book attempts 
to cover too large a field in a limited space. In consequence, the treatment of 
some important aspects of the heat transfer problem is rather brief ; for example, 
forced convection to a gas flowing at high speed is covered in a single paragraph. 

E. P. Hicks. 


ERRATA. 


In the article by K. C. Westfold, of the June issue, there should be added. 
to the References : 


Smurp, 8. F. and Wesrrop, K. C., 1949, Phil. Mag., 40, 831-48. 


[The Editors do not hold themselves responsible for the views 
expressed by their correspondents. | 


